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Abstract 

The last decade has witnessed great experimental progress that has led to mea- 
surements of near threshold cross sections — polarized as well as unpolarized — of 
high accuracy for various inelastic nucleon-nucleon collision channels. These data, 
naturally, pose challenges to theorists to develop methods by which they can be 
understood and explained in commensurate detail. 

In this work we review the status of the present theoretical understanding of 
this class of reactions with special emphasis on modeMndependent methods. We 
discuss in detail not only the many observables involved in the reactions, but also the 
physical questions that can be addressed by studying them in the various reaction 
channels. The special advantages of nucleon-nucleon induced reactions are stressed. 
Foremost among these is the use of the initial and final states as a spin/isospin filter. 
This opens, for example, a window into the spin dependence of the hyperon-nucleon 
interaction and the dynamics of the light scalar mesons. 
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1 Introduction 



1.1 Strong Interactions at low and medium Energies 

Quantum Chromo Dynamics (QCD), as the well-accepted theory of the Strong 
Interaction, unfolds its impressive predictive power in particular at high ener- 
gies. However, with decreasing momentum transfer the mathematical structure 
of QCD becomes increasingly complicated because the perturbative expansion 
in the coupling constant no longer converges. The fact that particles carrying 
strong charge (color) have never been observed in an isolated state is just one 
prominent example of the non-perturbative character of QCD at low energy 
Although a large amount of new data is available from measurements with 
electromagnetic probes (e.g., from MAMI at Mainz, ELS A at Bonn, and JLAB 
at Newport News), there is still much to be learned about the physics with 
hadronic probes at intermediate energies, comprising the investigation of pro- 
duction, decay, and interaction of hadrons. An important class of experiments 
in this context is meson production in nucleon-nucleon, nucleon-nucleus, and 
nucleus-nucleus coUisions close to the production thresholds. Recently, two 
reviews on this subject [1,2] were published, however both had their main 
emphasis on the experimental aspects. In this report we will focus on recent 
theoretical developments and insights. In addition in neither of these reviews 
was the potential of using polarization — which will be the main emphasis of 
this work — discussed in detail. 

It is interesting to ask a priori what physics questions can be addressed with 
the production of various mesons near threshold in pp, pd, and dd collisions 
with stored and extracted polarized beams. We therefore present here a brief 
list, naturally influenced by the personal taste of the author, that is not aimed 
at completeness, but more at giving the reader a flavor of the potential of 
meson production reaction in nucleonic collisions for gaining insight into the 
strong interaction physics at intermediate energies. 

• Final state interactions: In practice it is difficult to prepare secondary beams 
of unstable particles with sufficient accuracy and intensity that allow to 
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study experimentally the scattering of those particles off nucleons. Here 
production reactions are an attractive alternative. One prominent example 
in this context is the hyperon-nucleon interaction that can be studied in 
the reaction pp — > pYK, where Y denotes a hyperon (A, E etc.) and K 
denotes a kaon. Prom the invariant mass distributions of the YN system, 
information about the on-shell YN interaction can be extracted. We shall 
return to this point in sec. 2.2. 

• Baryon resonances in a nuclear environment: By design, in nucleon-nucleon 
and nucleon-nucleus collisions we study systems with baryon number larger 
than 1. This allows to study particular resonances in the presence of other 
baryons and excited through various exchanged particles. One example is 
the A^*(1535), which is clearly visible as a bump in any 1] production cross 
section on a single nucleon. A lot is known about this resonance already, 
however, many new questions can be answered by a detailed study of the 
reaction NN NNr]. In the closc-to-thrcshold regime, kinematics and 
the conservation of total angular momentum constrains the initial state to 
only a few partial waves. This, in combination with the dominance of the 
A^*(1535) in the 7] production mechanism, allows for a detailed study of 
the NN — > N*N transition potential. This is relevant to understanding 
the behavior of the N* in any nuclear environment as well as might reveal 
information on the structure of the resonance. For a detailed discussion we 
refer to sec. 7.3. 

In hadronic reactions there is additionally a large number of excitation 
mechanisms for intermediate resonances possible. Besides pseudoscalars and 
vector mesons, which are more cleanly studied in pion and photon induced 
reactions, scalar excitations are possible. It was shown recently that, for 
instance, the Roper resonance is rather easily excited by a scalar source [3], 
and thus hadronic reactions might be the ideal place to study this contro- 
versial resonance [4]. 

• Charge symmetry breaking (CSB): Having available several possible initial 
isospin states {pp, pn, dd, etc.), which all have several possible spin states, 
allows experiments that make the study of symmetry breaking easily ac- 
cessible. One example is the impact of charge symmetry breaking on the 
reaction pn — > d7r°: it leads to a forward-backward asymmetry in the dif- 
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ferential cross section (c.f. section 5 and Refs. [5,6]). It is important to note 
that in case of the pion production the leading charge symmetry breaking 
operators are hnked to the up-down mass difference. 

A forward-backward asymmetry in the reaction pn d{7r^r])(^s-~wave) 
should allow extraction of the CSB /o — oq mixing matrix element — a quan- 
tity beheved to give insights into the structure of the light scalar mesons. 

Once the leading symmetry breaking mechanism is identified, one can use 
the specific signals caused by CSB to extract information on symmetry con- 
serving matrix elements. For details on the idea as well as an application — 
extracting information on the existence of mesonic bound states — we refer 
the reader to Ref . [7] . 
• Effective field theory in large momentum transfer reactions: Pion production 
in nucleon-nucleon collisions is still amenable to treatment within chiral 
perturbation theory — if the expansion is adapted to the large momentum 
transfer typical of those reactions. This research is still in progress, however 
once completed it will allow not only to study systematically the CSB pion 
production mentioned above, but also to pin down the size of three body 
forces relevant for a quantitative understanding of pd scattering, and to 
include the dispersive corrections to nd scattering in a chiral perturbation 
theory analysis. The latter is necessary for an accurate extraction of the 
isoscalar ttN scattering length from deuteron reactions, which are otherwise 
difficult to access. This issue will be discussed in chapter 6.3. 

The wealth of information comes with the drawback that, apart from rare 
cases, it is difficult to extract a particular piece of information from the data. 
For example, resonances and final state interactions modify in a coherent su- 
perposition the invariant mass plot. Fortunately, polarization can act as a spin 
filter and different contributions to the interaction can be singled out because 
they show up in the angular distributions of different spin combinations. The 
various polarization observables will be discussed in detail in sec. 4. 

One characteristic feature of meson production in NN coUisions at the first 
glance looks hke a disadvantage with respect to a possible theoretical analy- 
sis: the large momentum transfer makes it difficult to reliably construct the 
production operator. However, as we will discuss in detail, due to this the pro- 
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a) 



b) 



Fig. 1. Comarison of two different reactions that can be used to extract parameters of 
a strong final state interaction. Shown is a typical low momentum transfer reaction 
(diagram a), e.g., K~d — > jAn) and a large momentum transfer reaction (diagram 
b), e.g., pp^pK+A). 

duction operator is largely independent of the relative energy of a particular 
particle pair in the final state ^ if we stay in the regime of small invariant 
masses. Because of this, dispersion relations can be used to extract low energy 
scattering parameters of the final state interaction and at the same time give a 
model-independent error estimate. In contrast, many reactions dominated by 
one-body currents that are characterized by small momentum transfers (c.f. 
Fig. 1). This makes a quantitatively controlled, model-independent analysis 
difficult. For example, it was shown in Refs. [8,9] that value for the neutron- 
neutron (AA^) scattering length extracted from 7r~d — >■ ^nn {K~d — > 7An), 
where the initial state is in an atomic bound state, is sensitive to the short- 
range behavior of the baryon-baryon interaction used ^ . 

1.2 Theoretical approaches to the reaction NN — > B1B2X 

In this chapter we will briefly describe the theoretical developments for the 
reactions NN — > B1B2X close to the threshold, where the Bi denote the two 

^ If there are resonances near by this statement no longer true. See the discussion 
in sec. 4. 

^ This implies that meson exchange currents should contribute at the same order 
of magnitude (c.f. discussion in sec. 6.4). 
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outgoing baryons (either NN or, e.g., in case of the production of a strange 
meson YN, where Y denotes a hyperon). 

The development of theoretical models especially for the reaction NN — > 
NNx has a long history. For a review of earlier works we refer the reader to 
the book by H. Garcilazo and T. Mizutani [10]. Most of the models can be 
put into one of two classes: namely those that use the distorted wave Born 
approximation, where a production operator that is constructed within some 
perturbative scheme is convoluted with nucleonic wavef unctions, and truly 
non-perturbative approaches, where integral equations are solved for the full 
{NN,NNx) coupled-channel problem. 

We start this section with a brief review of the second class. The obvious 
advantage of this kind of approach is that no truncation scheme is required 
with respect to multiple rescattering. This is a precondition to truly preserve 
three-body unitarity [11-13]. However, as we will see, this approach is tech- 
nically rather involved and there is significant progress still to be made that 
will allow for quantitative predictions, especially for reactions with unbound 
few-nucleon systems in the initial and final state. 

One obvious complication is to preserve the requirements of chiral symmetry, 
which are very important — at least for the reaction NN — > A^A^tt close to the 
threshold, as will be discussed in the following sections, since in its non linear 
realization it is necessary to consistently take into account one-, two- and 
more pion vertices as soon as loops are considered. Naturally, this can only 
be done within a truncation scheme that is consistent with chiral symmetry. 
As a possible solution, in Ref. [14] it was suggested to change to a linear 
representation, as given, e.g., by the linear sigma model. To our knowledge 
this area is not yet very well developed. 

All the calculations so far carried out for pion^ production within the full 
{NN,NNx) coupled-channel system, where performed within time ordered 
perturbation theory. Within this scheme, one way to impose unitarity is it to 

^ To be concrete in this section we will talk about pions only. They are, in fact, the 
only mesons that have so far been studied with integral equation approaches. 
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Fig. 2. Some possible contributions to the production operator for neutral pion pro- 
duction in pp collisions. Solid lines are nucleons, dashed lines are pions. 



demand, that the number of pions in any given intermediate state should not 
to exceed a given number. This, in combination with that there is no particle 
number conservation, implies the fact, that the tiNN coupling constant gen- 
erated by the equations at some intermediate time depends on the number 
of pions in flight at that time. The same is true for the renormalization of 
the nucleon pole. Put to the extreme: in the nuclear matter limit the dressed 
pion nucleon coupling vanishes [15] (see also the discussion in Ref. [16]). An- 
other problem inherently present in time-ordered formalisms was pointed out 
by Jennings [17]: namely that it is impossible simultaneously to fulfill the 
Pauli-Principle and the unitarity condition within these integral equations. A 
consequence of this is a wrong sign for T20 for elastic 7rd scattering at backward 
angles [18]. Both of these problems can be avoided within a covariant formu- 
lation of the integral equations as is given in Refs. [19,20]. Unfortunately, not 
only are those equations rather involved, they also need covariant T matrices 
for both subsystems. Although these exist (e.g., Ref. [21] for the nN system 
and Ref. [22] for the NN system), they were not constructed consistently. 
Also the three dimensional convolution approach of Refs. [23,16] does not suf- 
fer from the above mentioned problems, but also within this scheme so far no 
results for the {NN,NNx) were published. 

On the other hand there are the approaches that treat the production operator 
perturbatively (typically diagrams of the type shown in Fig. 2 are considered) , 
while including the nucleon-nucleon or, more generally, the baryon-baryon 
interaction in the final and initial state non-perturbatively. Therefore in this 
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formalism three-body unitarity can be achieved only perturbatively, however, 
this formalism at least fulfills the minimal requirement: to take fully into 
account the non-perturbative nature of the NN interaction. This is indeed 
necessary, as was confirmed experimentally, since a proper inclusion of the 
NN final state interaction is required to describe the energy dependence of 
the reaction pp — > ppX [1,2]. A production operator that is constructed per- 
turbatively allows in addition to treat the interactions between the subsystems 
without approximation and consistent with chiral symmetry (c.f. discussion 
in sec. 6.3). A priori, it is unclear according to what rules the production op- 
erator should be constructed and only a comparison with experiment can tell 
if the approach is appropriate or not. There is, however, one observation in 
favor of a pcrturbativc treatment of the production operator: in the case of 
pion production, effective field theory methods can be applied to show that 
as long as there is a meson exchange current in leading order it dominates the 
low energy amplitude. Loops undergo strong cancellations and thus might well 
be neglected. If, however, there is no meson exchange current at leading order 
(as is the case for the reaction pp ppir^) the situation is significantly more 
involved and loops might well be significant ^ . This issue will be discussed in 
detail in sec. 6. In this introductory chapter we only wish to stress, that to our 
present knowledge the case of 7r° seems to be an exception and, in addition to 
the reaction pp pn7r~^ all heavy meson production cross sections close to the 
threshold are indeed dominated by meson exchange currents ^ . We must add, 
however, that the of heavier mesons can also give a significant contribution. 

An approach that lies somewhere between those that solve the integral equa- 
tions for the whole NN NNtt system and those that treat the production 
operator perturbatively is that of Refs. [22,24,25]. These works document the 
attempt to extend straight-forwardly what is known about the phenomenol- 
ogy of nucleon-nucleon scattering (see for example [26]) to energies above the 
pion threshold. The three-body singularities stemming from both the energy 
dependent pion exchange as well as the nucleon and Delta self energies were 



^ Note, however, even then a perturbative treatment of the production operator is 
still justified. 

^ Here this phrase is to be understood to include resonance excitations as well. 
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taken into account. Via the optical theorem the pion production cross sections 
can then be calculated in a straight forward way. It turned out, however, that 
the models could not describe the close to threshold data. The reasons for 
this will become clear in chapter 6. In sec. 6.5 we will discuss how the recent 
progress in our understanding of pion production shows what is needed in 
order to improve, for example, the work of Ref . [24] . 

The construction of theoretical models for the production of mesons heavier 
than the pion started only recently, after the accelerator COSY came into 
operation. Most of these models are one meson exchange models. Those ap- 
proaches will be described to some extend in section 7.2. 

1.3 Specific aspects of hadronic meson production close to the threshold 

In the near-threshold regime the available phase space changes very quickly. 
Thus, especially when comparing different reactions, an appropriate measure 
of the energy with respect to that particular threshold is required. In pion 
production traditionally the variable f], defined as the maximum pion momen- 
tum in units of the pion mass, is used to measure the excess above threshold. 
For all heavier mesons the so-called excess energy Q, defined as 

is normally used. In Appendix A we give the explicit formulas that relate 
(5 to 77. It is useful to understand the physical meaning of the two different 
quantities: Q directly gives the energy available for the final state. The inter- 
pretation of 77 is somewhat more involved. In a non-relativistic, semiclassical 
picture the maximum angular momentum allowed can be estimated via the 
relation l^^x — Rq', where q' denotes the typical momentum of the correspond- 
ing particle and i? is a measure of the range of forces. If we identify R with 
the Compton wavelength of the meson of mass nix, we find Imax — q' /i^x — V- 
This interpretation was given for example in Refs. [27,28]. 

When trying to compare the cross sections of reactions with different final 
states one has to choose carefully the variable that is used for the energy. In 
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Fig. 3. Sketch of the production reaction showing the initial state interaction amongst 
the two nucleons (ISI), the final state interaction of all outgoing particles (FSI) as 
well as the production operator A as defined in the text. The left diagram indicates 
the complexity possible, whereas the right diagram shows the first and potentially 
leading term for the FSI only. 

Ref. [29] it is shown that the relative strength of the total cross sections for 
PP ~^ PP^^-i PP ~^ PPV PP ~^ PPV' is strikingly different when comparing 
them at equal 77 or at equal Q. Since the dominant final state interaction in all 
of these reactions is the pp interaction, it appears most appropriate to compare 
the cross sections at equal Q, for then at any given excess energy the impact 
of the final state interaction is equal for all reactions, which is not the case 
for equal values of rj. In Ref. [30] a different energy variable was suggested, 
namely the phase space volume. 

In order to produce a meson in nucleon-nucleon collisions the kinetic energy 
of the initial particles needs to be sufficiently large to put the outgoing meson 
on its mass shell. To produce a meson of mass nix an initial momentum larger 
than the threshold value 

|p-;(*'^'--)r = m.M^ + ^ (1) 

is necessary. As long as we stay in a regime close to the production threshold 
the momenta of all particles in the final state are small and therefore pf^^^^^ 
also sets the scale for the typical momentum transfer. In a non-relativistic 
picture a large momentum transfer translates to a small reaction volume, 
defined by a size parameter R ~ l/pi- Thus, already for pion production it! is 
as small as ~ 0.5 fm and is getting even smaller for heavier mesons. 
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Therefore the two nucleons in the initial state have to approach each other 
very closely before the production of the meson can happen. It then should 
come as no surprise that it is important for quantitative predictions to un- 
derstand the elastic as well as inelastic NN interaction in the initial state. 
However, since for all reactions we will be looking at the initial energy is sig- 
nificantly larger than the excess energy Q, the initial state interaction should 
at most mildly influence the energy dependence. The energy dependence of 
the production operator should also be weak, for it should be controlled by 
the typical momentum transfer, which is significantly larger than the typical 
outgoing momenta. 

On the other hand, in the near-threshold regime all particles in the final state 
have low relative momenta and thus undergo potentially strong final state 
interactions that can induce strong energy dependences. The different parts 
of the matrix element are illustrated in Fig. 3. 

In the next subsections we will briefly sketch some properties of meson produc- 
tion in nucleon-nucleon collisions in rather general terms with the emphasis 
on the gross features. Details, as well as selected results for various reactions, 
will be presented in the subsequent sections. 

1.4 Remarks on the production operator 

As was mentioned above one of the characteristics of meson production in 
nucleon-nucleon collisions is the large momentum transfer. This leads to a 
large momentum mismatch for any one-body operator that might contribute 
to the production reaction. 

This reasoning becomes most transparent when looking at it in real space, as 
shown in Fig. 4. In case of the one-body operator (panel a) the evaluation of 
the matrix element involves the convolution of a rapidly oscillating function 
with a mildly oscillating function, which generally will lead to significant can- 
cellations. If there is a meson exchange present, as sketched in panel b, the 
cancellations will not be as efficient. Thus one should expect that if there is 
a meson exchange current possible at leading order, it should dominate the 
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a) b) 

Fig. 4. Illustration of the momentum mismatch. The picture shows why one should 
a priori expect meson exchange currents to be very important close to the threshold. 

production process. This picture was confirmed in explicit model calculations 
(see e.g. Ref. [31]). 

As mentioned above one should in any case expect the production operator 
to be controlled by the large momentum transfer. Thus variations of the in- 
dividual amplitudes with respect to the total energy should be suppressed 
by (p'/pY, where p' (p) is the relative momentum of an outgoing to particle 
pair (of the initial nucleons). Therefore it should be neghgible close to the 
production threshold. This will be used for model-independent analyses of 
the production of heavy mesons and was checked within the meson exchange 
picture (see discussion in Ref. [32] and sec. 2.2). 

1.5 The role of different isospin channels 

One big advantage of meson production in nucleon-nucleon collisions as well 
as nucleon-nucleus and nucleus-nucleus collisions is, that besides the spin of 
the particles there is another internal degree of freedom to be manipulated: 
the isospin. 

Nucleons are isospin-(l/2) particles in an SU (2) doublet (isospin t — 1/2) with 
ts — +1/2 and ts — —1/2 for the proton (p) and neutron (n), respectively. 
Thus a two nucleon pair can be either in an isotriplet state (total isospin T — 1 
with the three possible projections of the total isospin — +1 for a pp state, 
— for a pn state and T'3 = — 1 for a nn state) or in an isosinglet state 



11 



(T = with T3 = for a pn state). 



Let us first concentrate on the isospin conserving situation. Even then various 
transitions are possible, depending on whether an isovector or an isoscalar 
particle is produced. Let us denote the allowed transition amplitudes by ^TiT^ , 
where 7i(Ty) denote the total isospin of the initial (final) NN system [27]. 
Then, in the case of the production of an isovector, the amplitudes An, Aiq, 
and Aqi are possible. They could be extracted individually from pp — > ppx'^ oc 
l^iiP) PP ~^ pnx'^ (X \Aii + Aiq\^, and pn — > ppx~ oc l/2|74ii + ^oiP; where 
the factor of 1/2 in the latter case stems from the isospin factor of the initial 
state (we have \pn) — l/-\/2(|10) + 100)) in which the isospin states are labeled 
by both the total isospin as well as their projection). On the other hand, for 
the production of an isoscalar there are two transitions possible, namely Aqo 
and All, where pp — > ppx oc and pn — > pnx oc l/2|Aii + AqoP- In this 

case the measurement of two different reaction channels allows extraction of 
all the production amplitudes. 

The deuteron is an isoscalar and thus acts as an isospin filter. Accordingly, one 
finds for the production of any isovector particle x, A{pp — > dx'^) — 2A{pn — > 
dx^) , as long as isospin is conserved. The consequences of isospin breaking will 
be discussed in sec. 5. 

The relative strength of the different transition amplitudes proved to provide 
significant information about the production operator. As an example let us 
look at the production of an isoscalar particle. If the production operator is 
dominated by an isovector exchange the corresponding isospin structure of the 
exchange current is Oiy — (71-72), where Tj denotes the isospin operator of 
nucleon i, which is clearly distinguished from the structure corresponding to 
an isoscalar exchange Ois ~ 1. Thus one finds ^ 

^ The easiest way to see this is to express Oj^ in terms of the Casimir operators of 
the underlying group: 

Tl-T2 = 4ti ■ t2 = 2(f 2 - £2 - tl) , 

with < >= T{T + 1) and < >= 3/4. 
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{TT^\0i^\TT:,) = 2T{T + l)-^ 



1 for T = 1 , 



< 



-3 for T = . 



On the other hand 



(TT3|0,,|TT3) 



1. 



For two given NN states a different total isospin imphes that also either the 
total spin or the angular momentum are different (c.f. next section). Therefore, 
different isospin states can not interfere in the total cross section and one may 
expect 



For X = 7] a recent measurement at Uppsala [33] gave for this ratio a value of 
6.5, clearly indicating a dominance of isovector exchanges. It should be stressed 
that the estimates given here only hold if effects from initial and final state 
interactions, the spin dependence of the production operator, as well as all 
other dynamical effects are completely neglected. However, even in the most 
detailed studies, the bulk of the ratio stems from the isospin factors [34-36]. 

1.6 Selection Rules for NN — > NNx 

In this section we will present the selection rules relevant for nucleon-nucleon 
induced meson production. Those rules are based on the symmetries of the 
strong interaction that imply the conservation of parity, total angular momen- 
tum and isospin. 

A two-nucleon system has to obey the Pauli Principle, which implies 
(^_^L+s+T _ ^_-|^^ where L, S and T denote the angular momentum, total 
spin and total isospin of the two nucleon system respectively. In the case of 
a two proton system, for example, where T—1, L + S needs to be even. Con- 
sequently, for T = 1 all even angular momenta, and therefore all even parity 
states, are spin singlet {S — 0) states and consequently have J — L. On the 
other hand, for T = states it is the odd parity states that are spin singlets. 



crtotjpn pnx) 
CTtotiPP ppx) 



5 if isovector exchanges dominate , 



1 if isoscalar exchanges dominate . 
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In addition, for a reaction of the type NN — > NNx we find from parity 
conservation 

(-)^ = nA-r-"'' , (2) 

where tTx denotes the intrinsic parity of particle x and L, L' and I' denote, 
respectively, the angular momentum of the incoming two nucleon system, of 
the outgoing two-nucleon system and of particle x with respect to the outgoing 
two-nucleon system. 

We may now combine the two criteria to write 

(_)(A5+AT)^^^(_y' ^ (3) 

where AS (AT) denotes the change in total (iso)spin when going from the 
initial to the final NN system. 

As an example let us look at the reaction pn — > d{T:r])s-wave — a reaction 
that should provide valuable information on the light scalar resonances (note: 
the nr] channel is the dominant decay channel of the scalar-isovector meson 
ao(980), c.f. chapter 7.5). Based on Eq. (9), the near-threshold regime is dom- 
inated by the lowest partial waves which in this case would be 5" = 1 and 
T" = (deuteron), together with a relative s wave between the two meson 
system and the deuteron (/' = 0). In addition, the (T^rii)s-wave system is a 
scalar-isovector system. Therefore we have to use Tr^: = +,/' = and AT = 1 
in Eq. (3), yielding AS* = 0. On the other hand the initial state needs to be of 
even parity, but even parity T = 1 states have 5* = 0. Thus, for the produc- 
tion of a scalar-isovector state the s-wave deuteron with respect to the two 
meson system is prohibited and, as long as isospin is conserved, this reaction 
has to have at least one p-wave in the final state. On the other hand, for the 
production of an isoscalar scalar particle, the s-wave final state is allowed. 

In a straight-forward way the corresponding selection rules for systems with 
nuclei in the initial and/or final state can be easily derived along the same 
lines. Bose symmetry, for example, demands a dd system to be symmetric 
under exchange of the two deuterons, forcing L + S* to be even. Therefore the 
lowest partial waves contributing to dd — > arj are ^Pq —>■ s and ^Di — > p. 

Another example of selection rules at work in a system other than NN — > 
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NNx. 
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The lowest partial waves for the production of a scalar particle in NN — ^ NNx. 
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a) b) c) d) 

Fig. 5. Diagrammatic presentation of the distorted wave Bom approximation. 
NNx will be given in sec. 4.3.1. 

2 The final state interaction 

The role of final state interactions in production reactions has been discussed 
in the literature for a long time and at various levels of sophistication. Ex- 
tensive discussions can be found in Refs. [37,38]. For more recent discussions 
focusing on the production of mesons in nucleon-nucleon collisions we refer to 
Refs. [39,40,32]. In this section we will approach the question of the role of the 
final state interactions from two approaches: first we will present a heuristic 
approach that should make the physics clear, but has only limited quanti- 
tative applicability; then, using dispersion integrals, quantitative expressions 
will be derived, but these are less transparent. This investigation will show 
under what circumstances scattering parameters can be extracted from high 
momentum transfer production reactions, and will address questions about 
the accuracy of this extraction. 

2.1 Heuristic approach 

As early as 1952 Watson pointed out the circumstances under which one would 
expect the final state interaction to strongly modify the energy dependence 
of the total production cross-section NN — > NNx [41]. His rather convincing 
argument went as follows: assume in contrast to the production reaction the 
pion absorption on a free two nucleon pair. Due to time reversal invariance the 
matrix elements for the two reactions are equal. Based on the observation that 
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the absorption takes place in a small volume, it is obvious that the reaction 
is a lot more probable if, in a first step, the two nucleons come together 
rather closely due to their attractive interaction and that this significantly 
more compact object then absorbs the pion. In this picture the infiuence of 
the two-nucleon system on the energy dependence of the cross section is very 
natural. Based on this argument, one should expect ^ 

fnixri 

(77) oc / d'q'p'dQ,,\^Uq')^NN{p')\' , (4) 
^0 

where q' {p') denotes the momentum of the outgoing meson (the relative mo- 
mentum of the outgoing two nucleon pair) and d?q' = dq'dflgi; energy conser- 
vation implies p'"^ = M]\j{E — Ug/), where ujqi denotes the energy of the meson. 
In Eq. (4) is the A^A^ wavefunction at the energy, 2E{p'), of the final 

NN subsystem, where E{p') = and denotes the wavefunction for the 
outgoing meson. If we now neglect the impact of the meson-nucleon interac- 
tion on the energy dependence, the meson wave function is given by a plane 
wave. For its component in the l-th partial wave we get [42] 

^l{q')^ji{q'R)c^q'' , (5) 

where ji{q'R) denotes the Ith Bessel function. The term q'^ is the first term in 
an expansion of q'R, where R represents the range of forces. Therefore, Eq. (5) 
should appropriately describe the q' dependence of the corresponding partial 
wave amplitude as long as (q'R) ^ 1. In addition, when the final state inter- 
action is very strong and shows a strong energy dependence (as, for example, 
the nucleon-nucleon interaction in the ^Sq for small relative energies), we may 
replace the NN wave function by the on-shell T matrix, which can be written 
in the L-th partial wave 

T\E,O = T,^e^'^^^^'^'^^^0i, (6) 

where dL^Ep/) denotes the corresponding phase shift as deduced from elastic 
NN scattering data. In practice this means that we neglect the term with 
the non-interacting outgoing two nucleon pair (Fig. 5a) und c)) compared to 
that with the pair interacting (Fig. 5b) and d)). In case of the NN final state 

^ Here all final state particles are treated non-relativistically for simplicity. 
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interaction this is indeed justified. In general, however, the plane wave piece 
is not negligible. Fortunately in this case dispersion integrals may be used to 
fix the relative strength of diagrams a and h of Fig. 5. This will be discussed 
in detail in the next section. 

At low energies the T matrix may be replaced by the first term in the effective 
range expansion ^ . 

p'^^^'coi5L{E,,)^- , (7) 

where ql denotes the scattering length of the A'^A^ interaction in the partial 
wave characterized by L. Thus we can write 

Since Eq. (5) only holds for (q'R) ^ 1, in Eq. (8) the a^^-term has to be 
kept in the denominator only if ai, ^ In case of the A" A" interaction, 

where R ~ this condition holds for the -S-waves only. When data for the 
reaction pp — > ppn^ close to threshold became available [43,44] Eq. 8 indeed 
turned out to give the correct energy dependence of the total cross-section 
[45]. Eq. (8) contains two very important messages. First of all that a strong 
final state interaction will significantly change the invariant mass distribution 
of the corresponding two particle subsystem — the relevant scale parameter is 
[aip'). This observation was used in Ref. [46] to extract information on the 
AA^ interaction. We come back to this in section 2. Secondly that due to the 
centrifugal barrier higher partial waves are suppressed in the close to threshold 
regime. As a rough estimate we can take non relativistic kinematics to derive 



(7 



\f^{NN)xJ 



for all those partial waves where there is no strong FSI. Here we used the 
relation between Q and rj as given in Eq. (A. 4) and /inn = M^/^ and fJ,(NN)x — 
rrix/ (1 + mx/ (2Mjv)) denote the reduced masses of the A" A" or more generally 



^ Unfortunately, in the literature there are two different conventions for the sign of 
the scattering length. The convention used here is that of Goldberger and Watson 
[37]; the usual convention for baryon-baryon interactions has a minus sign on the 
right hand side. 
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the baryon-baryon system and of the particle x with respect to the NN system 
respectively. 

Eq. (9) contains two pieces of information: first of all it shows that, for a 
given value of Q, heavier systems are more likely to be in a higher partial. 
The reason is that their maximum momentum allowed by energy conservation 
is larger than that of light particles. Secondly it should give a reasonable 
estimate of the energy dependence of the partial cross sections for all those 
partial waves where the NN system is in a partial wave higher than the S- 
wave. In case of the pion production this finding was confirmed recently in Ref. 
[47] for the reaction pp — > pp7r° as is demonstrated in Fig. 6, where the spin 
cross sections ^ctq and ^(Ji are shown as a function of r]. We use the notation 
of Ref. [47] as "^^^^(Tm, where S (m) denotes the total spin (spin projection) 
of the initial state. For the definition we refer to chapter 4.2. It is easy to 
show that the lowest partial waves in the final state that contribute to 
i^Ui) are Ps {Pp and -Ds), where capital letters denote the relative angular 
momentum of the two nucleons and the small letters that of the pion with 
respect to the two-nucleon system. Here we used that pion rf-waves should be 
suppressed compared to NN D-waves (c.f. Eq. (9)). Accordingly one should 
expect the Ps states to have an energy dependence given by rf whereas the 
Pp states should follow an rf dependence (c.f. Eq. (9)). The corresponding fits 
are shown in the Figure as a solid and a dashed line, respectively. The good 
agreement of these simple forms for the energy dependence of the different 
partial waves supports the conjecture that neither the initial state interaction 
nor the production operator and not even the FSI, for partial waves higher 
than jS-waves, induce a significant energy dependence beyond the centrifugal 
barrier. 

As we have seen, the energy dependence factorizes from the total amplitude. 
The situation is somewhat more complicated, however, when it comes to quan- 
titative predictions for the cross section: in the presence of very strong final 
state interactions it is not possible to separate model independently produc- 
tion operator and final state interaction. The reason for this is that what sets 
the scale for the size of the matrix element is the convolution integral of the 
production operator with the final state half off-shell NN T-matrix. This con- 
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Fig. 6. Energy dependence of two different spin cross sections for the reaction 
PP ~^ ppTT^- The data for ^ao (^cri), shown as solid (opaque) dots, are compared 
to a fit of the form y = Arf with r = 6 (t = %). The data are from table V of Ref. 
[47] with an additional error of 10 % on the total cross section. 

volution consists of a unitarity cut contribution in which the intermediate two 
nucleons are on their mass shell, and a principal value piece that is sensitive 
to the off-shell behavior of both the NN T-matrix as well as the production 
operator. This already shows that it is impossible to separate the production 
operator from the final state interaction in a model independent way. Since 
off-shell effects are not observable, a consistent construction of production 
operator and final state interaction is required to get meaningful results for 
the observables. 

In order to keep the formulas simple so far we have not discussed at all the 
effect of the Coulomb interaction, although its effect especially on near thresh- 
old meson production with two protons in the final state is known to be quite 
strong [45]. To account for the Coulomb interaction in the formulas given 
above one should use the Coulomb modified effective range expansion instead 
of Eq. (7). The corresponding formulas are e.g. given in Ref. [48]. 



21 



2.2 Dispersion theoretical approach 



The fact that meson production in nucleon-nucleon colhsions is a high mo- 
mentum transfer reaction allows for a rigorous treatment of strong final state 
interactions within a dispersion theoretical approach ^ . As we will see, in those 
cases where the effective range of a particular final state interaction is of the 
order of the scattering length, the modification of the invariant mass spectrum 
induced by the corresponding final state interaction is no longer proportional 
to the elastic scattering situation. However, the corrections that arise can be 
accounted for systematically. In addition to things that are are well known and 
presented in various text books [37,38], we derive an integral representation 
for the scattering length in terms of an observable [32]. In addition one gets 
at the same time an integral representation of the uncertainty of the method. 
The latter can be estimated using scale arguments or by doing a model calcu- 
lation. For example, in the case of the hyperon-nucleon interaction the scat- 
tering lengths as extracted from the available data on elastic scattering are of 
the order of a few fermi, but with an uncertainty of a few fermi — especially 
for the spin singlet scattering length (c.f. discussion in section 7.4). Using the 
kinematics for associated strangeness production to estimate of the theoretical 
uncertainty, we find a value of 0.3 fm. Thus it is feasible from the theoretical 
point of view to improve significantly our knowledge of the hyperon-nucleon 
scattering parameters through an analysis of the reaction pp — > pKA. 

To be definite, we will discuss the reaction pp pKA, however — as should 
be obvious — the complete discussion also applies to any subsystem with a 
sufficiently strong interaction. Let us assume the Ap system is in a single 
partial wave (this can be easily extracted from polarized experiments as will 
be discussed in sec. 4.3.1) and the third particle — here the kaon — is produced 
with a definite momentum transfer t = (pi —px)'^, where pi denotes the beam 
momentum. The full production amplitude may then be written as 



^ Actually, the formalism applies to all large momentum transfer reactions. In addi- 
tion to meson production in NN collisions, one may study e.g. ^yd — ^ BiB2+meson. 
This reaction is also discussed in Refs. [49-54]. 
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Fig. 7. Illustration of the sources of a discontinuity in the amplitude A from two 
baryon intermediate states. The dashed horizontal lines indicate the presence of a 
two baryon unitarity cut (Diagrams with initial state interaction are not shown 
explicitly (c.f. Fig. 5)). 

A{s,t,m)^- — — —dm +-/ — ^-dm , (10) 

TT J-oo m' ^ — m'^ TT Jml m' — m'^ 

where m^ = is the lefthand singularity closest to the physical region, 
mo = Mjv + Ma corresponds to the first right hand cut, and 

D{s, t, m^) = t, m^ + iO) - t, m^ - iO)). (11) 

Note that, in contrast to the representation of the amplitude in case of one 
open channel only as in NN — * tttt, D is not simply given by the imaginary 
part of A, since the initial state interaction also induces a phase in A (c.f. 
discussion in the next section). Thus D has to be written explicitly as the 
discontinuity of A along the appropriate two-particle cut. 

The second integral in Eq. (10) gets contributions from the various possible 
final state interactions, namely KK, NK and AA^, where the latter is the 
strongest. We may thus neglect the former two for a moment to get, for m^ > 
ml, 

D{s, t, m^) = A{s, t, m^)e-'^ sin b, (12) 

where (5 is the elastic kN phase shift. To see where this relation comes from we 
refer to Fig. 7: the appearance of the various unitarity cuts due to the baryon- 
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baryon final state interaction leads to a discontinuity in the amplitude along 
the corresponding branch cut. We thus may read almost directly off the figure 
that, indeed, the discontinuity in the amplitude is given by the production 
amplitude times the two-baryon phase space density times the (complex con- 
jugate of the) on-shell AN T-matrix^° . Eq. (12) follows by simply using the 
definition of the elastic phase shift: K,T{m'^) — e"^ sin S, where k, — p'tt/jl de- 
notes the phase space density here expressed in terms of the reduced mass of 
the AA^ system /i = {MnM^)/{Mn + Ma). 

We will discuss in sec. 7.4 how to control the possible infiuence of the K 
interactions. The solution of (10) in the physical region can then be written 
as (see [38] and references therein) 



A{s, t, m') = e"(^' + '^h{s, t, m') , (13) 
where 



TT J-oo m' ^ — 
and, in the absence of bound states. 



, , I S(m'^) , ,2 ,,,, 
u{z)^- , ' dm'^. 15 

The dependence of $ is dominated by the dependence of the production 
operator. The momentum transfer in the production operator, however, is 
controlled by the initial momentum and therefore one should expect the 
dependence of $ to be weak as long as the corresponding relative momentum 
of the A-nucleon pair is small. Thus, in a large momentum transfer reaction. 

The appearance of the complex conjugation here follows from the direct evaluation 
of the discontinuity of the T-matrix by employing the unitarity of the S matrix: 
from S^S = 1 it follows, that 

disc(r) := ^.{T{m^ + iO) - T{m^ - iO)) = -K\T{m^)\^ . 



24 



the dependence of the production amphtude A is governed by the elastic 
scattering phase shifts of the dominant two particle reaction in the final state! 
The relation between the phase shifts and the dependence of the amphtude 
as can be extracted from an invariant mass distribution is fixed by analyticity 
and unitarity! 

So far we are in line with the reasoning of the previous section, however the 
exponential factor in Eq. (13) is in general not simply the elastic scattering 
amphtude. For illustration we investigate the form of A{m'^) for a final state 
interaction that is fully described by the first two terms in the effective range 
expansion (c.f. Eq. (7)), p'ctg(5(m^)) = l/a+ ( 1/2) rp'^, where p' is the relative 
momentum of the final state particles under consideration in their center of 
mass system. Then A can be given in closed form as [37] 



where a — (l/r)(l + yl + 2r/a). Note: The case of an attractive interaction 
without a bound state in this convention corresponds to a positive value of a 
and a positive value of r. In the limit a — > oo, as is almost realized in NN 
scattering, the energy dependence of A{m'^) is given by 1/(1 — iap') as long 
as p' <^ 1/r. This, however, exactly agrees with the energy dependence for 
NN on-shell scattering. This prediction was experimentally confirmed by the 
near-threshold measurements of the reaction pp — > pp7r° [43,44]. However, for 
interactions where the effective range is of the order of the scattering length, 
the numerator of Eq. (16) plays a role and thus the full production amplitude 
is no longer given by the on-shell elastic scattering times terms whose energy 
dependence is independent of the scattering parameters. 

It is often argued that since the full production amplitude is given by a term 
with a plane wave final state as well as one with the strong interactions, that 
an appropriate parameterization is given by a Watson term (proportional to 
the elastic scattering of the outgoing particles) plus a constant term, and their 
relative strength should be taken as a free parameter. The considerations in 
the previous paragraphs show, however, that this is not the case: Eq. (13) (and 
thus also the special form given in Eq. (16)) describes the dependence of 



A{s,t,m'^) 



(p'2 + Q;2)r/2 



$(s,i,m') , 



(16) 



l/a+(r/2)p'2_^p- 
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Fig. 8. The integration contours in the complex m' ^ plane to be used in the original 
treatment by Geshkenbein [55,56] (a) and in the derivation of Eq. (18) (b). The 
thick lines indicate the branch cut singularities. 

the full invariant mass spectrum. 

One more comment is in order: In the previous section we stressed that, al- 
though the shape of the invariant mass spectra as well as the energy depen- 
dence of the total cross section are governed by the on-shell interactions in 
the subsystems, the overall normalization is not. This statement, based on 
effective field theory arguments, is illustrated in Appendix D. On the other 
hand, Eq. (13) gives a closed expression for the invariant mass spectrum for 
arbitrary values of m^. On the first glance this look like a contradiction. How- 
ever, it should be stressed that Eq. (13) does in general not allow one to relate 
the asymptotic form of the full production amplitude with information on the 
scattering parameters of the strongly interacting subsystems in the final state 
to the amplitude in the close to threshold regime. Besides the trivial obser- 
vation that Eq. (13) holds for each amplitude individually and that far away 
from the threshold there should be many partial waves contributing, there 
is no reason to believe that the function $ is constant over a wide energy 
range. Therefore, we want to emphasize that the dependence of Eq. (13) is 
controlled by the scattering phase shifts of the most strongly interacting sub- 
system only in a very limited range of invariant masses. In fact, the leading 
singularity that contributes to an dependence of $ is that of the t-channel 
meson exchange in the production operator. The rn? dependence that origi- 
nates from this singularity can be estimated through a Taylor expansion of 
the momentum transfer and thus turns out to be governed by {p' /pY, where p' 
denotes the relative momentum of the final state particle pair of interest. For 
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p we may use that of the initial state particles at the production threshold, all 
to be taken in the over all center of mass system. As a consequence, for values 
of that do not significantly deviate from ml (or equivalently p' <^ p) the 
assumption of $ being constant is justified. However $ can show a significant 
dependence over a large range of invariant masses. 

Let us now return to our main goal: namely to derive from (13) a formula 
that allows extraction of the elastic scattering phase shifts from an invariant 
mass spectrum. In this course we will derive a dispersion relation in terms 
of the function |A(m^)p. Here we will use a method similar to the one used 
in [55,56]. However, in contrast to the formulas derived in these references, 
we will present a integral representation for the elastic scattering phase shifts 
from production data that involves a finite integration range only. 

For this we first observe that Eq. (15) holds for purely elastic scattering only 
and thus is of very limited practical use. On the other hand a significant 
contribution stems from large values of m! ^ which depends only weakly on 
in the near-threshold region, and therefore can be absorbed into the function 
Let 



exp 



1 r 

TT Jm 



5(m'2) 



m' ^ — rn? — iO 



dm 



I 2 



(17) 



where ^{ml^^^.m^) = $(m2)$„2 (m^), with 



*m2 

'"'max 



[m^) = exp 



-/ 

TT Jn, 



5{m'^) 



I 2 



iO 



-dm! 



I 2 



(18) 



The quantity m^^^ is to be chosen by physical arguments in such a way 
that both $(m^), and $m2^^^(m^) vary slowly on the interval (^0''"^max)- 
Obviously, it needs to be sufficiently large that the structure in the amplitude 
we are interested in can be resolved . On the other hand, it should be as 
small as possible, since this will keep the infiuence of the inelastic channels 



One immediately observes that the integral given in Eq. (22) goes to zero for 
i^max i^o- At the same time the theoretical uncertainty Sa"^""^'', defined in Eq. 



(24), tends to the value of the scattering length. 
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small. In order to estimate the minimal value of rrimax we see from Eq. (16) 
that the values of p' that enter in the integral given in Eq. (22) should be 
at least large enough that the scattering length term plays a significant role. 
Thus we require p'^ax ~ ^/O'typ- In case of the pn interaction (the spin triplet 
pn scattering length is 5.4 fm) this would correspond to a value of p'^ax 
only 10 MeV corresponding to a value for e^ax — ''^max ~ "^o — the maximum 
excess energy that should occur within the integral in Eq. (17) (as well as Eq. 
(22) below) — as small as 2 MeV. On the other hand, if we want to study the 
hyperon-nucleon interaction, we want to 'measure' scattering lengths of the 
order of a few fermi we choose atyp = 1 fm, leading to a value of 40 MeV for 
^max- Note that 40 MeV is still significantly smaller than the thresholds for 
the closest inelastic channels (that for the KYjN final state is at 75 MeV, that 
for the TiKKN channel at 140 MeV). 

The integral in (17) contains an unphysical singularity of the type 
log {ra^rnax " '"^^), which is canceled by the one in ^(m^^^., m^), but this does 
not affect the region near threshold. 

Notice next that the function 



has no singularities in the complex plane except the cut from to ni^^^ 
(c.f. Fig. 8b) and its value below the cut equals the negative of the complex 
conjugate from above the cut. Hence, 



log {A{m 



1 




(19) 



5(m2) 





dm! 



I 2 



(20) 



It is an important point to stress that 
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'max 



m' ^ — mo(m' ^ — m 



1 



2^1 \ ml 




dm' 



I 2 







(21) 



as long as is in the interval between m^ and m^^^,. Therefore, if the func- 
tion $ only weakly depends on m^, as it should in large momentum transfer 
reactions, it can well be dropped in the above equation. In addition, up to 
kinematical factors, \A\'^ agrees with the cross section for the production of 
a Ap pair of invariant mass m. We can therefore replace it in Eq. (20) by 
the cross section where, because of Eq. (21), all constant prefactors can be 
dropped. Thus we get 



This is the desired formula: namely, the scattering length is expressed in terms 
of an observable. Note that Eq. (22) is applicable only if it is just a single YN 
partial wave that contributes to the cross section as- In sec. 4.3.1 we will 
discuss how the use of polarization in the initial state can be used to project 
out a particular spin state in the final state. 

Up to the neglect of the kaon-baryon interactions, Eq. (20) is exact. Thus, 
^g(th) — ^YiQ theoretical uncertainty of the scattering length extracted using Eq. 
(22) — is given by the integral 



Since log |$(m^„^, m^)^ = log |$(m2)p + log |$^2^^Jm2)|2 we may write 
§Q{th) _ §(j{ihc) _|_ ^^TUmax^ whcrc thc former, determined by $(m^), is con- 
trolled by the left hand cuts, and the latter, determined by $^2^^^(m^), by 
the large-energy behavior of the AA^" scattering phase shifts. The closest left 




(22) 





(23) 
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hand singularity is that introduced by the production operator, which is gov- 
erned by the momentum transfer. Up to an irrelevant overall constant, we may 
therefore estimate the variation of $ ~ 1 + 5{p' /pf, where we assume 5 to be 
of the order of 1. Evaluation of the integral (23) then gives 

Sa^"^'^^ ~ Sip'^Jp') ~ 0.05 fm , 

where we used p'^^x ^ — '^f^^max with emax ~ 30 MeV and the threshold value 
for p ~ 900 MeV. On the other hand, using the definition of $^2^^^(m^) given 
in Eq. (18) one easily derives 



\Sa" 



8{y)dy 



(1 + ^2) (3/2) 



< ^^\5max\ , (24) 
T^Pmax 



where — {m — rrimax) I ^max- Thus, in order to estimate 5a"^"""^ we need to 
make an assumption about the maximum value of the elastic h.N phase for 
> Tn^^^. Recall that we implicitly assume the inelastic channel not to play 
a role. This assumption was confirmed in Ref. [32] within a model calculation. 
The denominator in the integral appearing in Eq. (24) strongly suppresses 
large vahies of y. Since for none of the existing AA^ models does 5max exceed 
0.4 rad, we estimate 

^^m^a. _ 0.2 fm . 

When using the phase shifts as given by the models directly in the integral, 
the value for SaJ^""^'' is significantly smaller, since for all models the phase 
changes sign at energies above rn^^^. Combining the two error estimates, we 
conclude 

5a < 0.3 fm . 

In the considerations of the theoretical uncertainty of the parameters extracted 
so far we did not talk about the possible effect of the kaon-baryon interaction. 
Unfortunately it can not be quantified a priori; however, it can be controlled 
experimentally. As will be discussed in sec. 4, a significant meson-baryon 
interaction, due, for example, to the presence of a resonance, will show up 
as a band in the Dalitz plot. If this band overlaps with the region of final 
state interactions we are interested in, interference effects might heavily distort 
the signal [57]. By choosing a different beam energy, the FSI region and the 
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resonance band move away from each other. Thus there should be an energy 
regime where the FSI can be studied undistorted. It is therefore necessary for 
a controlled extraction of FSI parameters to do a Dahtz plot analysis at the 
same time. An additional possible cross-check of the influence of a meson- 
baryon interaction would be to do the full analysis at two different beam 
energies separated in energy by at least the typical hadronic width (about 
100-200 MeV). If the parameters extracted agree with each other, there was 
no substantial influence from other subsystems. 

The possible gain in experimental accuracy for the extraction of the AA^ scat- 
tering lengths is illustrated and discussed in detail in sec. 7.4. 

To our knowledge so far the corresponding formalism including the Coulomb 
interaction is not yet derived. 

3 The initial state interaction 

As mentioned earlier, the collision energy of the two nucleons in the initial 
state is rather large, especially when we study the production of a heavy me- 
son. However, as long as the excess energy is small, only very few partial waves 
contribute in the final state and thus the conservation of total angular momen- 
tum as well as parity and the Pauli principle allow for only a small number 
of partial waves in the initial state. In addition, only small total angular mo- 
menta are relevant. In such a situation the standard methods developed by 
Glauber, which include the initial state interaction via an exponential sup- 
pression factor [58] , cannot be apphed. The role of the initial state interaction 
relevant to meson production reactions close to the threshold was discussed 
for the first time in Ref. [40]. 

It should be clear a priori that for quantitative predictions the ISI has to play 
an important role: to allow a production reaction to proceed the nucleons in 
the initial state have to approach each other very closely — actually signifi- 
cantly closer that the range of the NN interaction. Thus, especially for the 
production off heavy mesons, a large number of elastic and inelastic NN reac- 
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tions can happen before the two-nucleon pair comes close enough together to 
allow for the meson production. Thus the effective initial current gets reduced 
significantly. 

Technically the inclusion of the initial state interaction in a calculation for a 
meson production process requires a convolution of the production operator 
with the half off-shell T-matrix for the scattering of the incoming particles 
(c.f. Fig. 5c) and d)). Contrary to the final state interaction, for the high 
initial energies the energy dependence of the A'"A'" interaction is rather weak. 
One might therefore expect that the real part of the convolution integral is 
small and that the effect of the initial state interaction is dominated by the 
two nucleon unitarity cut [40]. This contribution, however, can be expressed 
completely in terms of on-shell NN scattering parameters: 



where p denotes the relative momentum of the two nucleons in the initial state 
with the total energy E, and Sl {rji) denote the phase shift (inelasticity) in 
the relevant partial wave L. Each partial wave amplitude should be multiplied 
by Xl, defined in Eq. (25), in order to account for the dominant piece of the 
ISI. This method was used e.g. in Ref. [59]. Using typical values for phase 
shifts and inelasticities at NN energies that correspond to the thresholds of 
the production of heavier mesons, Eq. (25) leads to a reduction factor of the 
order of 3, thus clearly indicating that a consideration of the ISI is required 
for quantitative predictions. One should keep in mind, however, that Eq. (25) 
can only give a rough estimate of the effect of the initial state interaction 
and should whenever possible be replaced by a full calculation. This issue is 
discussed in detail in Ref. [34] . 

Unfortunately, the apphcability of Eq. (25) is hmited to energies where scat- 
tering parameters for the low NN partial waves are available. Due to the 
intensive program of the EDDA collaboration for elastic proton-proton scat- 
tering [60] and the subsequent partial wave analysis documented in the SAID 
database [61], scattering parameters in the isospin-one channel are available 




(25) 
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up to energies that correspond to the production threshold. The situation 
is a lot worse in the isoscalar channel. Here pn scattering data would be re- 
quired. At present, those are available only up to the 77 production threshold. 
It is very fortunate that there is a proposal in preparation to measure spin 
observables in the pn system in this unexplored energy range at COSY [62] . 

4 Observables 

In this section the various observables experimentally accessible are discussed. 
After some general remarks about the three body kinematics in the final state, 
in the next subsection we will discuss unpolarized observables. In the following 
subsection we will then focus on polarized observables. 

Even if x denotes just a single meson, the reaction NN B1B2X is subject 
to a five dimensional phase: three particles in the final state introduce 3 x 
3 = 9 degrees of freedom, but the four-momentum conservation reduces this 
number to 5. As we will restrict ourselves to the near-threshold regime, the 
final state can be treated non-rclativistically. The natural coordinate system is 
therefore given by the Jacobi-coordinates in the overall center of mass system 
(c.f. Ref. [42]), where first the relative momentum of one pair of particles is 
constructed and then the momentum of the third particle is calculated as the 
relative momentum of the third particle with respect to the two body system. 
Obviously there are three equivalent sets of variables possible. In the center 
of mass this choice reads 



where we labeled the three final state particles as ijk; qk' — Pk holds in the 
over all center of mass system only. For simplicity in the following we will 
drop the subscripts when confusion is excluded. For reactions of the type 
NN — > B1B2X it is common to work with the relative momentum of the two- 
nucleon system and to treat the particle x separately. From the theoretical 



Pij' 



Mipj - MjPi 

Mi + M.j 





(26) 



Mi + Mj + Mfe 
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Fig. 9. Illustration of the choice of variables in the over all center-of-mass system. 

point of view this choice is most convenient, for one is working already with the 
relative momentum of the dominant final state interaction and the assignment 
of partial waves as used in sec. 1.6 is straightforward. Note that for any given 
relative energy of the outgoing two nucleon system e = p'^/Mjv the modulus 
of the meson momentum is fixed by energy conservation; we thus may 
characterize the phase space by the 5 tuple 

e = {e,Qp,,Q,,} , (27) 

where VLk = {cos{9k), (pk) denotes the angular part of vector k. The coordinate 
system is illustrated in Fig. 9. The center of mass nucleon momentum in the 
initial state will be denoted by p. Throughout this report we choose the beam 
axis along the z axis. Explicit expressions for the vectors appearing are given 
in Appendix B.l. 

Due to the high dimensionality of the phase space it is not possible to present 
the full complexity of the data in a single plot. At the end of sec. 4.2 we will 
discuss a possible way of presenting the data through integrations subject to 
particular constraints. A different choice would be simply to present highly 
differential obscrvables as is done for bremsstrahlung , or at least to use the 
two dimensional representation of the Dalitz plot to show some correlations. 

It should be noted that in the early bremsstrahlung experiments the method of 
integration was not possible, because the detectors used had very limited angular 
acceptance. 
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Fig. 10. Sketch of a Dalitz plot for a reaction with a three-body final state (particles 
are labeled as 1,2 and 3). Regions of possibly strong final state interactions in the 
{ij) system are labeled by F(^ij^; a resonance in the {12)-system would show up like 
the band labeled as R. The region of a possible interference of the two is labeled as 
I. 



4-1 Unpolarized observables 

As long as the initial state is unpolarized, the system has azimuthal rotation 
symmetry, reducing the number of degrees of freedom from 5 to 4. 

In the case of a three-particle decay (a 1 ^ 3 reaction) the physics does 
not depend on the initial direction. The same is true for reactions where the 
initial state does not define a direction, like in the experiment series for pp 
annihilation at rest carried out by the Crystal Barrel collaboration at LEAR 
(see Ref. [63] and references therein). Therefore the number of degrees of 
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freedom is further reduced from four to two . Those are best displayed in 
the so-called Dahtz plot (see [64] and references therein) that shows a two 
dimensional representation in the plane of the various invariant masses m^^ = 
{pi+PkY- In reactions with two particles in the initial state and three particles 
in the final state (2 — > 3 reactions), however, the initial momentum defines an 
axis and therefore a single fully differential plot is no longer possible. Especially 
in the example given below it will become clear how the appearance of the 
additional momentum changes drastically the situation. 

To account for the higher complexity of the 2 — > 3 reactions, if enough statis- 
tics were collected in a particular experiment, one can present differential 
Dalitz plots — a new plot for each different orientation of the reaction plane 
[65]. However, in what follows we will call Dalitz plot the representation of a 
full data set/calculation in the plane of invariant masses, ignoring the initial 
direction. Obviously this means throwing out some correlations, for in general 
any integration reduces the amount of information in an observable. We will 
briefly review the properties of a Dalitz plot, closely following Ref. [64]. In 
Fig. 10 a schematic picture of it is shown. In this plot also different regimes 
are specified: those of potentially strong final state interaction for the subsys- 
tem (ij) are labeled F(^ijy They should occur when particle i and j move along 
very closely. On the other hand, resonances will show up as bands in the Dalitz 
plot. Labeled as R in the figure, the effect of a resonance in the (12) system is 
shown. The total area of the Dalitz plot scales with the phase space volume. 
Therefore, as the excess energy decreases, resonance and final state interaction 
signals or the regions of different final state interactions might start to overlap, 
leading to interference phenomena (c.f. the hatched area / in Fig. 10). It was 
demonstrated recently [57] that those can rather strongly distort resonance 
and final state interaction signals. At least for known final state properties, 
these patterns might help to better pin down resonance parameters [57]. On 
the other hand, as the excess energy increases the different structures move 
away from each other and one should then be able to study them individually. 



This is quite obvious, since one can look at a three particle decay as the crossed 
channel reaction of two body scattering which is well known to be characterized by 
two variables in the unpolarized case. 
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This observation is of great relevance if one wants to extract parameters of a 
particular final state interaction from a production reaction (c.f. discussion in 
section 2). 

In case of particle decays of spinless or unpolarized particles, the Dalitz plot 
not only contains the information about the occurrence of a resonance, but 
also its quantum numbers can be extracted by projecting the events in the 
resonance band (labeled as R in Fig. 10) on the appropriate axis (for the 
example of the figure this is the 23 axis). In the case of 2 — > 3 reactions, 
however, this projection is not necessarily conclusive. To explain this statement 
we have to have a closer look at the angles of the system. First of all there is 
a set of angles, the so called helicity angles, that can be constructed from the 
final momenta only. One example is 



cos(^^,,0 = , (28) 

\p \\q \ 

were p' and q' where defined in Eq. (26). These angles that can be extracted 
from the Dalitz plot directly. In addition there are those angles that are related 
to the initial momentum — the Jackson angles. One example is 



^°«(^.'.) = §JJ§\ ■ (29) 

It should be stressed that it is not in the distribution of the helicity angles but 
that of the Jackson angles that the subsystems reflect their quantum num- 
bers [64] . Therefore a pure Dalitz plot analysis is insufficient for production 
reactions and the distributions for the Jackson angles have to be studied as 
well. This will be illustrated in an example in the following subsection. 

Note that in the presence of spin there are even additional axes in the problem. 
This will be discussed in detail in sec. 4.2. 

^^As is stressed in Ref. [64], only under special conditions, namely for peripheral 
production as it occurs in high energy experiments, the information in the helicity 
angles and the Jackson angles agrees. Close to threshold, however, meson production 
is not at all peripheral. 
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Fig. 11. Angular distributions for the reaction pp dK^K^ measured at Q = 
MeV [66]. The solid line shows the result of the overall fit including both KK 
s-waves as well as p-waves. To obtain the dashed (dotted) line, the parameters 
for the p-wave (s-wave) were set to zero (see text). The small error shows the sta- 
tistical uncertainty only, whereas the large ones contain both the systematic as well 
as the statistical uncertainty added linearly. 

4.1.1 Example: analysis of pp — > dK^K^ close to threshold 

Recently, a first measurement of tlie reaction pp dK^K^ close to the pro- 
duction threshold was reported [66] at an excess energy Q = 46 MeV. The 
data, as well as the corresponding theoretical analysis, based on the assump- 
tion that only the lowest partial waves contribute, will now be used to illustrate 
the statements of the previous section. It will become clear, especially, that 
the information encoded in the distributions of the Jackson angles and the 
helicity angles is rather different. 

As can be seen in table 2, a final state that contains s-waves only is not 
allowed in this reaction. For later convenience we work with the relative 
momentum of the kaon system {pkk' = p'] c.f. Eqs. (26)) and the deuteron 
momentum with respect to this system {q^' = q'). Given our assumptions, 
that only the lowest partial waves contribute, the amplitudes that contribute 
to the production reaction are either linear in q' or linear in p'. In Ref. [67] 

In sec. 7.5 it will be argued, tfiat the reaction pp dK^K^ can be used to study 
scalar resonance aQ'(980). Thus we are especially interested in the partial waves of 
the kaon system, that should show a strong final state interaction. 
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Table 3 

Results for the C parameters from a fit to the experimental data. The parameters 
are given in units of Cq . 









c{ 


C2 + ^Ca 


OzbO.l 


1 zb 0.03 


1.26 ± 0.08 


-0.6 ±0.1 


-0.36 ±0.17 



the full amplitude was constructed (c.f. also discussion in sec. 7.5), however 
it should be clear that the spin averaged square of the matrix element can be 

written as 

\M\' = CW" + Cgp" + C({q' ■ pf + Ct{p' ■ pf 

+C2{p' ■ q') + Cs{p' ■ p){q' ■ p) , (30) 

since all terms in the amplitude are either linear in p' or linear in q'. Here 
p — p/\p\ denotes the beam direction. Since the two protons in the initial state 
are identical, any observable has to be symmetric under the transformation 
p — —p. This is why p appears in even powers only. 

Figs. 11 and 12 show the data as well as a fit based on Eq. (30). The parameters 
extracted are given in table 3. The first two panels of Fig. 11 contain the 
distributions of the angles {p' ■p)/p' and {q' ■p)/q'- The last panel contains the 
distribution of the helicity angle {p' ■ q')/{q'p')- The solid line corresponds to a 
complete fit to the data including both p-waves in the K^K~^ system as well as 
those in the d{K^K~^) system. For the long dashed line the K^K^ p- waves were 
set to zero, whereas for the dotted hue the d{K^K^) p-wmes (corresponding 
to KK s-waves) were set to zero. Thus, the first two panels truly reflect the 
partial wave content of the particular subsystems individually, whereas the 
hehcity angle (which can also be extracted from the Dahtz plot) shows a flat 
distribution only if both subsystems are in a p-wave simultaneously. Therefore, 
the helicity angle can well be isotropic although one of the subsystems is in a 
high partial wave. Note that this statement is true even if all particles were 
spinless. The only thing that would change is that Cq and Cq would vanish 
(c.f. subsec. 4.3.2). In Fig. 12 two Dalitz plot projections (invariant mass 
spectra) are shown. The flrst one {da/drnxx) is needed to disentangle Cq 
and Cq . The second one does not give any additional information. 
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M^^ [GeV] M^^ [GeV] 

Fig. 12. Various mass distributions for the reaction pp dK^K~^ (Line code as 
in Fig. 11). The small error bars show the statistical uncertainty only, whereas the 
large ones contain both the systematic as well as the statistical uncertainty ( c.f. Ref. 
[66]). 

What is now the information contained in the two-dimensional Dahtz plot? 
Since it does not contain any information about the initial state, the parts of 
the squared amplitude that can be extracted from the Dalitz plot are easily 
derived from Eq. (30) by integrating over the beam direction [64], giving 

/ dn,\M\^ = (c^ + ic?') q" + {c^,' + icf') p" 

+ (Cs + ^^4) ip' ■ q') . (31) 

Therefore from the Dalitz plot one can extract the total KK s-wave strength 
{Ci + (l/3)Cf ), the total KK p-wave strength (C^' + (l/3)Cf'), as well as 
the strength of the interference of the two (C3 + (1/3)C74). Note that in this 
particular example, all the coefficients given in Eq. (31) (and even the Cq 
and Ci individually) can as well be extracted from the angular distributions 
given in Fig. 11 and the KK invariant mass distribution (left panel of Fig. 
12) directly; the Dalitz plot here does not provide any additional information. 

Obviously, as we move further away from the threshold, the complexity of the 
amplitude increases and the Dalitz plot contains information not revealed in 
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the projections. To summarize, in order to allow for a complete analysis of 
the production data, in addition to the Dalitz plot the angular distributions 
of the final momenta on the beam momentum need to be analyzed as well. 
The latter distributions are the ones that give most direct access to the partial 
wave content of the subsystems. 

4.2 Spin dependent observables 

Polarization observables for 2 ^ 2 reactions are discussed in great detail in 
Ref. [68]. In our case, however, we have one more particle in the final state 
and therefore there are more degrees of freedom available. Here we will not 
only derive the expressions for the observables in terms of spherical tensors 
but also relate these to the partial wave amplitudes of the production matrix 
elements. In this section we closely follow Refs. [47,69]. 

In terms of the so called Cartesian polarization observables, the spin- 
dependent cross section can be written as 



where (Tq{^) is the unpolarized differential cross section, the labels i,j and 
k can be either x, y or and P^, Pt and Pf denote the polarization vector 
of beam, target and the first one of the final state particles, respectively. All 
kinematic variables are collected in ^, defined in Eq. (27). The observables 
shown explicitly in Eq. (32) include the beam analysing powers Ajo, the cor- 
responding quantities for the final state polarization Dq^, the spin correlation 
coefficients A^j, and the spin transfer coefficients D^j. In this context it is im- 
portant to note that baryons that decay weakly, as the hyperons do, have a 
self analyzing decay. In other words, the angular pattern of the decay particles 
depends on the polarization of the hyperon, therefore, the hyperon polariza- 



(j(e, A, Pu Pf) = ^o(0 1 + E((^f)^Ao(0 + {Pf)iDoim 



i 



+ Y^muPthMo + (n).(P/),A,(0) 




(32) 



ijk 
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tion in the final state can be measured without an additional polarimeter (see 
e.g. Ref. [70]). All those observables that can be defined by just exchanging 

— * — * 

Pf) and Pt, such as the target analyzing power Aq^, are not shown explicitly. 
Prom Eq. (32) it follows that for example 

where the af"^ (<^j*^) the Pauli matrices acting in the spin space of beam 
and target, respectively. The production matrix element is denoted by Al. In 
addition, Sb (st) denote the total spin of the beam (target) particles. 

It is straightforward to relate the polarization observables to the partial wave 
amplitudes that can be easily extracted from any model. Por this purpose it 
is convenient to use spherical tensors defined through 

n-kzqzMqi _ \ 71 r (/i)t (/2)tA/< (6) (<) ka]\ (oa\ 

■^kiqi,k2q2 - (^2sf, + l){2St + 1) ^ "^"^ ''^ •'^^^kiqjk2q2'^^^ J ' ^"^^^ 

where the T^q denote the spherical representation of the spin matrices 



TlO = (^z , n±l = T;^(cra; ± ic^y) , ^00 = 1- (35) 

To relate the observables to the spherical tensors, the easiest method is to use 
the definitions of Eqs. (35) inside the various Eqs. (34). Por the observables 
for which the final polarization remains undetected, the relations between the 
various T and the corresponding observables arc shown in Table 4. In Table 5 a 
few of the observables that contain the final state polarization are listed. Triple 
polarization observables arc not listed explicitly, but it is straightforward to 
derive also the relevant expressions for these, such as 

^xx,x + ^yy,x ~ {^xy,y ~ ^yx,y) = Iie{Tiii_i) . 

After some algebra given explicitly in Appendix C, one finds 

m^) = ^E^|.apK^. (36) 

LIX 

where p = {/cigi, k2q2: hqa, k^q^} and Q = qi + q2 - qs - qi- All the angular 
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Table 4 

Relations between spherical tensors and some observables that do not contain the 
final state polarization following Ref. [69]. To simplify notation, the indices specify- 
ing the final state polarization are dropped. The symbol * indicates a possible set of 
independent observables. Note: For pp induced reactions more observables become 
equivalent, as described in the text. Those are marked by (*). The linear combi- 
nations of spin correlation observables appearing in the table are defined in Eqs. 
(39). 43 
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Table 5 

Relations between spherical tensors and some observables that do contain the final 
state polarization. The symbol * indicates a possible set of independent observables. 
The linear combinations of spin correlation observables appearing in the table are 
defined in Eqs. (40). 
dependence is contained in 



Air 



and 



(37) 



(38) 



^kA-EC|g''M"(M'^)t. 

a, a 

Here a and a are multi-indices for all the quantum numbers necessary to 
characterize a particular partial wave matrix element and C denotes a coupling 
coefficient that can be expressed in terms of Clebsch-Gordan coefficients. Its 
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explicit form is given in Eq. (C.9). Prom Eq. (36) one can derive the angular 
dependences of all observables. 

In what follows it is convenient to define the following quantities: 

— -^xx ~l~ -^yy ) — -^xx -^yy ) and — A^-y Ay^. (39) 
and, analogously, 

-Ds = + Dyy , Da = Dxx - Dyy , and Ds = D^y - Dy^ . (40) 

Using the conservation of parity and the explicit expression for the C coeffi- 
cient given in Eq. (C.9) 

(ja,a,p _ ^_-j(fci+fe2+fe3+fc4)^a^a,p {^^) 

Since the parameter C is real, the analyzing powers are proportional to the 
imaginary part of M^M?, whereas the differential cross section as well as the 
spin correlation parameters depend on the real part of M.aM.*a- Thus, it is 
either the real part or the imaginary part of B that contributes to the angular 
structure. As we will see in the following subsection, this observation allows 
for a straight-forward identification of the possible azimuthal dependences 
of each observable. Another obvious consequence of Eq. (41) is, that those 
observables for which ki is odd have to be small when only a single partial 
wave dominates. Thus, at the threshold, analyzing powers will vanish. 

The structure of Eq. (36) is general — no assumption regarding the number 
of contributing partial waves was necessary. However, if we want to make 
statements about the expected angular dependence of observables, the number 
of partial waves needs to be restricted. For example, if we allow for at most 
p-waves for the A'^A'^ system as well as the particle x with respect to the A^A^ 
system, then the largest value of L and / that can occur is 2, which strongly 
limits the possible ^-dependences that can occur in the angular function B 
defined in Eq. (37). 

If the spin of the particles is not detected there is no interference between 
different spin states in the final state. This leads to a severe selection rule 
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in case of the reaction pp — > ppX: since the final state is necessarily in a 
T — 1 state the Pauli principle demands that different spin be accompanied 
by different angular momentum. Therefore, for the a reaction with a pp final 
state the partial waves can be grouped into two sets, namely {Ss, Sd, Ds, Dp} 
and {Pp, Ps}, where only the members of the individual sets interfere with 
each other. 

4-2.1 Equivalent observables 

All the angular dependence of the observables is contained in the function B 
defined in Eq. (37). In this subsection we will discuss some properties of B 
and relate these to properties of particular observables. 

The functional form of B enables one immediately to read off the allowed 
azimuthal dependences for each observable as well as to identify equivalent 
observables. To sec this we rewrite Eq. (37) as 

^S,a(«^^) =E/L[,A,Q,n(^P',^.')exp{i[(g-ri)0,,+ri0,,]} , (42) 

n 

which directly translates into the following ^-dependences for the spherical 
tensors (c.f. Eq. (36)): 

AT 

%{pA)^ ^p,n(ep/,eg/)exp{i[(Q-n)0p/ + n0g/]} . (43) 

n=-N 

Note that N is given by the highest partial waves that contribute to the 
reaction considered: 

-Lmax <{Q-N) < Lmax and N < Imax , (44) 

where L^ax {lmax) is given by twice the maximum baryon-baryon (meson) an- 
gular momentum. These limits are inferred by the Clebsch-Gordan coefficient 
appearing in the definition of B in Eq. (37). 

Eq. (43) directly relates the real and the imaginary parts of the spherical 
tensors: 

lm{Tp{ep>, <Pp, + 7r/(2g), 0,, + 7r/{2Q))) = Re{T,{9p>, (j^p,, 9^,, </),0) • (45) 
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Thus, two observables are equivalent if they are given by the real and imagi- 
nary parts of the same spherical tensor with Q 7^ .In table 4 the relations 
of the various observables to the spherical tensors are given. Thus, using Eq. 
(45) we can identify the following set of pairwise equivalent observables: 

^2/0 = -^xO -I -^Qy = -^Ox i -^xx ~ -^yy = -^xy + -^yx i (^6) 

and analogously for observables for which the final state polarization is mea- 
sured as well. Notice that there is no connection between A^x + Ayy and 
A^y — Ayx, for these have Q = and therefore there is no transformation, 
such as the one given in Eq. (45) , that relates real and imaginary parts of the 
spherical tensors. 

For identical particles in the initial state, as in pp induced reactions, all observ- 
ables should be equivalent under the exchange of beam and target. This further 
reduces the number of independent observables, for now the beam analyzing 
powers are equivalent to the target analyzing powers and A^z is equivalent to 
Azx- In tables 4 and 5 a possible set of independent observables is marked by a 
*. Those of these that are not independent for identical particles in the initial 
state are labeled as (*). 

From the discussion in the previous section it follows (c.f. Eq. (41)), that all 
those observables with an even (odd) value of ki + k2 + + k^ lead to a 
real (imaginary) value for gfp^„, defined in Eq. (43). As a consequence, for all 
coefficients appearing in the expansion of the observables, the 0-dependence 
is fixed (c.f. table (4)); for example, the terms that contribute to ctq, Cq^zz, 
and (TqAy, behave as cos{n{(f)qi — (pp')). 

As was stressed above, the phase space for the production reactions is of 
high dimension. To allow for a proper presentation of the data as well as 
of calculations, one cither needs high dimensional plots (sec discussion in the 
previous section) or the dimensionality needs to be reduced to one dimensional 
quantities , while, however, still preserving the full complexity of the data. As 

^^For the experimental side this is the far more demanding procedure, for the 
angular dependence of efficiency as well as acceptance needs to be known very 
well over the full angular range for those variables that are integrated in order not 
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Fig. 13. Some polarization observables reported in Ref. [47] for the reaction 
pp ppn^ at Tj = 0.83 as a function of the pion angle compared to predictions 
of the model of Ref. [71]. The solid lines show the results for the full model whereas 
contributions from the Delta where omitted for the dashed lines. 



can be seen in Eq. (43), each polarization observable is described by 2N + 1 
functions gp,n{Gp',Gq'), where the number of relevant terms is given by the 
number of partial waves. In order to allow disentanglement of these functions, 
in Ref. [47] it was proposed to integrate each observable over both azimuthal 
angles under a particular constraint, 

$ = m([)p + n(t)q = c . (47) 

This integration projects on those terms that depend on $ or do not show any 
azimuthal dependence at all [47]. To further reduce the dimensionality of the 
data, either the relative proton angle or the meson angle can be integrated to 
leave one with a large number of observables that depend on one parameter 
only. Those are then labeled as Afj{6k), where k is either p or q. In Figs. 13 

to introduce false interferences. 
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Fig. 14. Some polarization observables reported in Ref. [47] for the reaction 
pp pp-K^ at T] = 0.83 as a function of cos{9p) compared to predictions of the 
model of Ref. [71]. The solid lines show the results for the full model whereas con- 
tributions from the Delta were omitted for the dashed lines. 



and 14 some observables reported in Ref. [47] are shown for the energy with 
highest statistics, namely t] = 0.83. In the figures the data is compared to 
the model predictions of Ref. [71]. The solid lines are the results of the full 
model whereas for the dashed lines the contribution from the Delta isobar was 
switched off. In section 6 this model will be discussed in more detail. 

In the case of Ref. [47] the complete set of polarization observables for the 
reaction pp —>■ ppn^ is given. Since the particles in the initial state are identical 
there are 7 independent observables, all functions of 5 independent parameters 
(c.f. table 4). In the analysis of the data presented in the same reference it 
was assumed that only partial waves up to p-waves in both the A^A^ as well 
as the {NN)'n' system were relevant. Thus the various integrations described 
in the previous paragraph lead to 32 independent integrated observables that 
depend only on a single parameter. On the other hand, if the assumption 
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about the maximum angular momenta holds, only 12 partial waves have to be 
considered in the partial wave analysis. Since the amplitudes are complex and 
two phases are not observable (c.f. discussion at the end of section 4.2), a total 
of 22 degrees of freedom needs to be fixed from the data. Thus a complete 
partial wave decomposition for the reaction pp ppir^ seems feasible. 

4-3 General structure of the amplitudes 

In this section we give the recipe for constructing the most general transition 
amplitude for reactions of the type A'^A'^ B1B2X, where we focus on spin 

1/2 baryons in the final state. A generalization to other reactions is straight- 
forward. For further applications we refer to a recent review [72]. 

For simplicity let us restrict ourselves to those reactions in which there is only 
one meson produced. The system is then characterized by three vectors, 

p, q' , and p' , 

denoting the relative momentum of the two nucleons in the initial state, the 
meson momentum, and the relative momentum of the two nucleons in the final 
state, respectively — in the over all center of mass system. In addition, as long 
as X denotes a scalar or pseudoscalar meson, we find 6 axial vectors, namely 
those that can be constructed from the above: 

i{p X p') , i{p X q') , and i{p' x q') ; 

and those that contain the final or initial spin of the two nucleon system 

S , S' , and i{SxS') , 

where S — Xi'^y^X2 and S' — xl^'^yixl)'^ ■ Here Xi denotes the Pauli spinors 
for the incoming (1,2) and outgoing (3,4) nucleons and a denotes the usual 
Pauli spin matrices. If x is a vector particle, an additional axial vector, namely 
the polarization vector of the vector meson e * occurs. In addition, if instead of 
a two nucleon state in the continuum a deuteron occurs in the final state, its 
polarization direction will be characterized by the same e*. Since the energy 
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available for the final state is small (we focus on the close to threshold regime) , 
we restrict ourselves to a non relativistic treatment of the outgoing particles. 
This largely simplifies the formalism since a common quantization axis can be 
used for the whole system. 

In order to construct the most general transition amplitude that satisfies parity 
conservation, we have to combine the vectors and axial vectors given above 
so that the final expression form a scalar or pseudoscalar for reactions where 
the produced meson has positive or negative intrinsic parity, respectively. The 
most general form of the transition matrix element may be written as 

M = H{IX') + iQ ■ {SX') + iA ■ {S'X) + {Si Sj ')Bij , (48) 

where X = {x2'^yXi) ^^(^ ^' = (x4'^y(X3 )^) • addition, the amplitudes have 
to satisfy the Pauli Principle as well as invariance under time reversal. This 
imposes constraints on the terms that are allowed to appear in the various 
coefficients. 

The 9 amplitudes that contribute to B may be further decomposed according 
to the total spin to which S and S ' may be coupled 

Bij = lf5ij + hlcijk + h\j , 

where the superscripts indicate if the combined spin of the initial and the final 
state are coupled to (s), 1 {y) or 2 (i), where h\j is to be a symmetric, trace 
free tensor. 

Once the amplitudes are identified the evaluation of the various observables 
is straightforward. In this case the polarization comes in through 

Xixl = ^(l + ^i-<?) , 

— * 

where Pi denotes the polarization direction of particle i. Using the formulas 
given in Appendix B.2 one easily derives (summation over equal indices is 
implied) : 
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4^0 = \H\^ + + \A\^ + \Bmn? , (49) 
AAoiao = +ie,,fe (Q*Qfc + 5*5^^) + 2Im (S^^l^ - Q*//) , (50) 

AD^,ao = -ie,,k {A*Ak + B*^Bik) - 2Im (5*Qz - A*if) , (51) 

+ 2Re {eiM^H - A*^Bi^) - Q\Qj - B*^B^^) , (52) 
4AjO-o = 2Re (Q*Aj + eumQlBmj 

+^jmlBilA*^ + -eUmtjnkBinBmk + B*jH^ , (53) 

4A,-,feo(Xo = Im(5i,(2iJ*Afc - 2Q;fi;fe - e,,g,(A;A^ - 5;„5;^)) 
+ '^^iij{Q*iAk — enmkA*^Bi„i + HBi^) 

+2{Q*Bjk + QjBik) 

iB*mBjn + B*M) (54) 



For illustration we also give here the explicit expressions for As and defined 
in the previous section: 



(ToAj: = - {-\H\'' + |g,|2 - I A|2 + 

CTqAa = - - (QIQx - Q*yQy + Bl^B^rn " B*^By^^ 



As was stressed in the previous section, the method of spherical tensors is very 
well suited to identifying equivalent observables. This is significantly more 
difficult in the amplitude approach. However the amplitude method becomes 
extremely powerful if — due to physical arguments or appropriate kinematical 
cuts — one of the subsystems can be assumed in an s-wave, for then the num- 
ber of available vectors is reduced significantly and rather general arguments 
become possible (c.f. sec. 4.3.1). 

Since any amplitude can be made successively more complex by multiplying 
it by an arbitrary scalar, in most of the cases an ordering scheme is demanded 
in order to make the approach useful. In the near-threshold regime this is 
given by the power of final momenta occurring — in analogy with the partial 
wave expansion. Actually, the amplitude approach presented here and the 
partial wave expansion presented in the previous subsection are completely 
equivalent. However, in the near-threshold regime the amplitude method is 
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etot(deg) etot(deg) 

Fig. 15. Sensitivity of the analyzing power as well as the differential cross section 
for the reaction pn — > pp-K~ to the sign of the ^Pq Sqs amplitude ai. The lines 
corresponds to the model of Ref. [73]: the solid line is the model prediction whereas 
for the dashed line the sign of ai was reversed. The experimental data are from Ref. 
[74] and [75] at TLab = 353 MeV (rj = 0.65). 

more transparent. As one goes away from the threshold the number of partial 
waves contributing as well as the number of the corresponding terms in the 
amplitude expansion increases rapidly. As a consequence the construction of 
the most general transition amplitude is rather involved. The partial wave 
expansion, on the other hand, can be easily extended to an arbitrary number 
of partial waves. 

As follows directly from Eq. (48), in the general case the matrix element Ai 
is described by 16 complex valued scalar functions. One can show, e.g. by ex- 
plicit construction, that for general kinematics of the reaction NN NNx 
all 16 amplitudes are independent. A possible choice is given explicitly in 
Table 6 for the reaction pp pp-l-(pseudoscalar). However, for particu- 
lar reaction channels or an appropriate kinematical situation their number 
sometimes reduces drastically. For example in coUinear kinematics (where 
p , q' , and p' are all parallel) the number of amplitudes that fully de- 
scribes the reactions pp — * ]9]9+(pseudoscalar) is equal to 3 (this special case 
is discussed in Ref. [76]). This can be directly read off Table 6, for under 
coUinear conditions all cross products vanish and all structures of one group 
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Table 6 

List of a possible set of the independent amplitude structures that contribute to 
the reaction pp pp+ (pseudo scalar). In the last column shows the lowest partial 
waves for the final state that contribute to the given amplitude (using the notation 
of section 1.6). To keep the expressions simple we omitted to give the symmetric, 
trace free expressions for the terms listed in the last two lines. 



Amplitudes 


Structures 


Lowest pw 


H 


{pxp')q'{p-p') 


Ds 


Q 


P 


p'ip-p') 


q'ip- q') 


Ss, Ds, Sd 


A 


P' 


p{p-p') 


q'ip'-q') 


Ps, Pd 


Bij 


eijkPkip' ■ q') 


(-ijkp'kip- q') 


eijkq',{p' ■ P) 


Pp 




Sij{pxp')q' 


ip X p')iq'j 


{P X q%p'j 






PiPjip X p')q' 


{p X p')iPj{p- q') 


{p X q')iPj{p-p') 





that are given by the different vectors of the system collapse to one structure 



Another interesting example is that of elastic pp scattering. The Pauli Principle 
demands (c.f. section 1.6) that odd (even) parity states are in a spin triplet 
(singlet). Time reversal invariance requires the amplitude to be invariant under 
the interchange of the final and the initial state. When parity conservation is 
considered in addition, one also finds that the total spin is conserved. Thus, 
only H and Bij will contribute to elastic pp scattering. In addition, from 
the two vectors available in the system (p and p') one can construct only 4 
structures that contribute to the latter, namely 5ij{p ■ p')bi, eijk{p x p')kb2, 
{PiP'j+PjPi-{Vmj){p-p')h, and {PiPj+p'jP',-{2/3)p\j){p-p')h.'' Thus, 
pp scattering is completely characterized by 5 scalar functions. 

As a further example and to illustrate how the formalism simplifies in the vicin- 
ity of the production threshold, we will now discuss in detail the production of 
pions in NN collisions. Throughout this report, however, the formalism will be 

Note, this choice of structures in not unique; we could also have used {p x p')i{p x 
p')j as a replacement of any other structure in 6* . 
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applied to various reactions. In our example there are three reaction channels 
experimentally accessible, namely pp — > pp7r°, pp — > pmr'^, and pn — > ppTr~, 
that can be expressed in terms of the three independent transition amplitudes 
ATiTf (as long as we assume isospin to be conserved), where Ti(Tf) denote 
the total isospin of the initial (final) NN system [27] (c.f. section 1.5). As in 
section 1.6, we will restrict ourselves to those final states that contain at most 
one p-wave. We may then write for the amplitudes of An, 



H 



11 



0, 



A^^ = a2p' + a3 
for the amplitudes of ^lo. 



(p • p')p - -p' 



(55) 



?io 



~*10 



{q' ■p)p- ^q' 



for the amplitudes of ^oi, 



(56) 



■■ ciq" + C2 



pip ■ q") - ^q' 



PiPr, 



(57) 



where p denotes the initial NN momentum, normalized to 1, and p' and 
q' denote the final nucleon and pion relative momentum, respectively. The 
coefficients given are directly proportional to the corresponding partial wave 
amplitudes as listed in Table 1; e.g., 03 is proportional to the transition ampli- 
tude — >^ P2S. When constructing amplitude structures for higher partial 
waves care has to be taken not to list dependent structures. In order to remove 
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dependent structures the reduction formula Eq. (B.7) proved useful. In addi- 
tion, one should take care that the number of coefficients appearing exactly 
matches the number of partial waves allowed. For example, the partial waves 
that contribute to are ^Si -^^ Pis, ^Di -^^ Pis and P2S. 

The large number of zeros appearing in the above list of amplitude structures 
reflects the strong selection rules discussed in section 1.6. 

As an example we will calculate the beam analyzing power and the differential 
cross section for the reaction pn pp7i~ . These observables were measured 
at TRIUMF [74,75] and later at PSI [77,78]— here, however, with a polarized 
neutron beam). In accordance with the TRIUMF experiment, where the rela- 
tive A^A^ energy in the final state was restricted to at most 7 MeV, we assume 
that the outgoing NN system is in a relative S-w&Ye. This largely simphfies 
the expressions. We then find 



(7o = ||ai|^ ^?'R-e (al (^ci + ^C2^^ cos(6') , 
aoAy = (^^q'Hm{clc2) sm{20) 

-^q'lm (at (ci - ^2)) sin(^)) cos(0) , (58) 

where we used the definitions (g' x p)y = —q' sin{6) cos{(f)) and (g' ■ p) = 
q' cos{6) (c.f. Appendix B.l). Thus, the forward-backward asymmetry in ao 
as well as the shift of the zero in Ay directly measure the relative phase of 
the ^Pq Sqs transition in An (ai) in the pion p-wave transitions of Aqi 
(ci and C2)), as was first pointed out in Ref. [79]. This issue will be discussed 
below (c.f. sec. 6). Note: As before we neglected here pion (i-waves, since they 
are kinematically suppressed close to the threshold (c.f. Eq. (9)). 



4.3.1 Example I: Polarization observables for a baryon pair in the ^Sq final 
state 

As an example of the efficiency of the amplitude method, in this subsection we 
will present an analysis of the angular pattern of some polarization observables 
for the reaction pp — > B1B2X under the constraint that the outgoing two 
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baryon state {B1B2) is in the ^Sq partial wave and x is a pseudoscalar, as is 
relevant for the reaction pp — > pKA. The dependence of the observables on the 
meson emission angle is largely constrained under these circumstances, as was 
stressed in Ref. [80]. In contrast to the discussion in sec. 1.6, here we will not 
assume Bi and B2 to be identical particles. The results of this subsection will 
show how to disentangle in a model independent way the two spin components 
of the hyperon-nucleon interaction [32] . 

The analysis starts with identifying the tensors that are to be considered in 
the matrix element of Eq. (48) . For this we go through a chain of arguments 
similar to those in sec. 1.6. Given that we restrict ourselves to a spin-zero final 

— * 

state, only H and Q can be non-zero. In addition, the quantum numbers of 
the final state are fixed by l^;, the angular momentum of the pseudoscalar with 
respect to the two baryon system, since 



J^l, and TTtoi = (-)(^-+^) (59) 

for the total angular momentum and the parity, respectively. Conservation of 
parity and angular momentum therefore gives 



(-)^ = {-y^+^ = {-y+^ s^i, (60) 

and consequently we get H — 0. In addition, for odd values of Ix we see from 
the former equation that L must be even. In pp systems, however, even values 
of L correspond to -S" = states, not allowed in our case. Therefore must 
be even. We may thus make the following ansatz: 



Q^ap + Pq'{q'-^ , (61) 

where a and (3 are even functions of p, p' and {p ■ p') . All other coefficients 
appearing in Eq. (48) vanish. This has serious consequences for the angular 
dependences of the various observables. For example, the expression for the 
analyzing power collapses to 

Note, here we could as well refer to table 6 to come to the same conclusion as in 
this section, for the state is allowed for the pp system. However, the argument 
given is quite general and instructive. 
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^ijk 



{Q*Qk) = -lm{P*a){q'-p){q' x p}i ■ 



(62) 



Therefore, independently of the partial wave of the pseudoscalar emitted, for 
a two-baryon pair in the ^5*0 state the analyzing power Agy vanishes if the 
pseudoscalar is emitted either in the xy-plane or in the 2;x-plane. In Ref. [32] 
this observation was used to disentangle the different spin states of the AN 
interaction. 

4-3.2 Example 11: Amplitude analysis for pp — > dK'^K^ close to threshold 

In sec. 4.1.1 we discussed in some detail the data of Rcf. [66] for the reaction 
pp dK^K^ based on rather general arguments on the cross section level. In 
this subsection we will present the corresponding production amplitude based 
on the amplitude method presented above. This study will allow us at the 
same time to extract information on the relative importance of the Oq and the 
A(1405) in the reaction dynamics. 

As in sec. 4.1.1, we assume that either the KK system or the deuteron with 
respect to the KK system is in a p-wave, whereas the other system is in an 
s-wave, calling for an amplitude linear in p' or g', respectively We use the 
same notation as in sec. 4.1.1. Therefore the final state has odd parity and 
thus also the amplitude needs to be odd in the initial momentum p. An odd 
parity isovector A'^A'^ state has to be S* = 1 and thus has to be linear in 
defined in section 4.3. In addition, the deuteron in the final state demands 
that each term is linear in the deuteron polarization vector e. We therefore 
get for the full transition amplitude, in slight variation to Eq. (48) due to the 
presence of the deuteron in the final state. 



where e appears as complex conjugate, since the deuteron is in the final state, 
and 



Bij^ aspPiQj' + bsp{p-q')6ij + CspQi'pj + dspPiPj{q' ■ p) 

+ apsPiPj' + bps{p-p')5ij + cpsPi'pj + dpsPiPj{q' ■ p) , (64) 



M = BijSit-j 



(63) 
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p p p p 

a) b) 

Fig. 16. Illustration of the coordinate system used in the analysis for the reaction 
pp^dK°K+. 

where capital letters in the amplitude label indicate the partial wave of the 
KK system and small letters that of the deuteron with respect to the KK 
system. 

Once the individual terms in the amplitude are identified, it is straightforward 
to express the Cj defined in sec. 4.1.1 in terms of them. We find, for example. 



Cl'^l{\as,\'+\cs,\') , 



2 

Cf = \bsp\'^ + - \bsp + dsp]"^ + Re a*spCsp + {asp + csp)*{bsp + dsp) 

^2 = 0'*SpO'Ps + C*SpCps ■ 

A fit to the experimental data revealed that, within the experimental uncer- 
tainty, Cq is compatible with zero. Thus, given the previous formulas, both 
aps and cps have to vanish individually. 

Based on the strongly populated KK s-waves in Ref. [66] it was argued that 
the reaction pp — > dK^K'^ is governed by the production of the Oq . In Ref. 
[81], however, it was argued that the strong Kd interaction caused by the 
proximity of the A(1405) resonance should play an important role as well. This 
FSI should enhance the Kd s-wave. We now want to calculate the contribution 
of this partial wave relative to the Kd p-wave based on the amplitudes given 
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in table 3. This will illustrate a further strength of the amplitude method, for 
within this scheme changing the coordinate system is trivial. The coordinate 
system suited to study resonances in the KK system and in the Kd system 
are illustrated in the left and right panels of Fig. 16. All we need to do now is 

— * — * 

express the vectors that appear in Eq. (30) in terms of Q' and P'. We find 

q'^P'-aQ' and p' ^ ^{{2 - a)Q' + P') , 

where a = mrf/(mrf + m^). Obviously, the squared amplitude expressed in the 
new coordinates reveals the same structure as Eq. (30): 

\M\^ = B^'q'' + B^'P" + B^'iQ' ■ pf + b({P' ■ pf 

+B2{P' ■ Q') + Bs{P'-P){Q'-P) , (65) 

where the coefficients appearing can be expressed in terms of the C coefficients 
of Eq. (30) so that, for example. 



Q, _ i2-af , , a{2-a) 1 
-Do — ^ <-^0 +0i -O2 , 



Bf = ^^^C( + a^Cr - ^i^iCs . 



(66) 



With these expressions at hand it is easy to verify, that KK s-waves contribute 
to 83 % to the total cross section, whereas Kd s-waves contribute to 54 % 
only [82]. Here we used the total s-wave strength for KK, (Cg + (l/3)Ci ), 
and the total s-wave strength for Kd, {Bq + (1/3)5^ ), as a measure of the 
strength of the partial waves. 

As we do not see a significant population of the Kd s-wave, it appears that 
the A(1405) does not play an essential role in the reaction dynamics of pp — > 
dK^K'^ close to threshold in contrast to the Oq . 



4-4 Spin Cross Sections 



As early as 1963, Bilenky and Ryndin showed [83], that from the spin cor- 
relation coefficients that can be extracted from measurements with polarized 
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Table 7 

List of the lowest partial waves in the final state that contribute to the individual 
spin cross sections. Capital letters denote the baryon-baryon partial waves whereas 
small letters that of the meson with respect to the baryon-baryon system. 





possible final state 

Six 


s for pp bb'x 

Plx 




'S,p 
^Sos, ^Sod, ^Sid 
^Sis, ^Sod, ^Sid 


^PjS, ^PlS 
^PjP, ^Pip 
^Pjp, ^Pip 



beam and target, the cross section can be separated into pieces that stem from 
different initial spin states. Their results were recently re-derived [84] and the 
formalism was generalized to the differential level in Ref . [85] . With these so- 
called spin cross sections it can easily be demonstrated how the use of spin 
observables enable one to filter out particular aspects of a reaction. We begin 
this subsection with a derivation of the spin cross sections using the amplitude 
method of the previous subsection and then use the reaction ^ — > pmr'^ as 
an illustrative example. 

Since we have the amplitude decomposition of the individual observables given 
in Eqs. (49)-(54), one easily finds 

(7o(l - - Ayy - A,,) = \H\'' + =: Vq , (67) 
(7o(l + Axx + Ayy - A,,) = |g,|2 + =: Vo , (68) 

ao(l + A,,) = ^ (|g,|2 + \Qy\^ + \Bxn\' + \Byn\^) =: Vi , (69) 

where the assignment of the various spin cross sections ^'^^'^^^cfms: with S (Ms) 
the total spin (projection of the total spin on the beam axis) of the initial state 
can be easily confirmed from the definition of the amplitudes in Eq. (48). 

As was shown in section 1.6, in case of two nucleon initial or final states 
restrictive selection rules apply. For example, for pp final states the isospin 
of the final NN system is 1 and therefore states with even (odd) angular 
momentum have total spin (1). From Table 7 it thus follows, that the pp S- 
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Fig. 17. Demonstration of the selectivity of the spin cross sections. Shown are the 
spectra of the reaction pp pnir^ as a function of e for an excess energy of 25 
Me V. In each panel the solid line shows the full result for the corresponding cross 
section, the dot-dashed, long-dashed, dashed, and dotted lines show the Sp, Ss, Sd, 
and Pp contribution respectively. The curves are from the model of Ref. [71 ]. 

wave in connection with a meson s-wave contributes only to ^(Jq. This example 
of how spin observables can be used to filter out particular final states was 
used previously in sec. 2 (c.f. discussion to Fig. 6). 

For illustrative purposes we show in Fig. 17 the spectra for the various cross 
sections for the reaction pp —>■ pmr^, as a function of the relative energy of the 
nucleon pair in the final state. The curves correspond to the model of Ref. [71] 
that very well describes the available data in the vr"*" production channel. The 
model is described in detail in sec. 6. In the upper left panel the unpolarized 
cross section is shown. It is dominated by the Sp final state (the dominant 
transition is Sip), and from this spectrum alone it would be a hard 

task to extract information on final states other than the ^Si NN state. This 
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Fig. 18. The lowest meson production thresholds for single meson production in 
proton-proton collisions together with the corresponding energy ranges of the modem 
cooler synchrotrons. 

can be clearly seen by the similarity of the shapes of the dot-dashed line and 
the solid line. The spin cross sections, however, allow separation of the spin 
singlet from the spin triplet initial states. Naturally ^ctq (lower left panel of Fig. 
17) is now saturated by the Sp final state, but in ^(Tq and ^ai other structures 
appear: the former is now dominated by the transition ^Pq — >^ SqS and the 
latter by ^Pi Sis. 



4-5 Status of Experiment 



In this presentation we will be rather brief on details about current as well 
as planned experiments, as this subject was already covered in recent reviews 
[1,2]. Here we only wish to give a brief list of observables and reactions that 
are measured already or are planned to be measured in nucleon-nucleon and 
nucleon-nucleus induced reactions. 

In the case of pion production, measurements with vector- and tensor- 
polarized deuteron and vector-polarized proton targets and polarized proton 
beams have been carried out at lUCF [134]. 

Because of good An detection of photons and charged particles, CELSIUS, at 
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Table 8 

List of observables measured for various NN — > NNx channels for excess energies 



uptoQ = 40 MeV. 
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least in the near future, is well equipped for studies involving rj mesons. For 
the production of heavier mesons, COSY, due to its higher beam energies and 
intensities, but due most importantly to the possibility to use polarization, is 



64 



in a position to dominate the field during the years to come. The energy range 
of the different cooler synchrotrons is illustrated in Fig. 18. 

In table 8 a list is given for the various A^A^ induced production reactions mea- 
sured in recent years at SACLAY, TRIUMF, PSI, COSY, lUCF, and CELSIUS 
in the near-threshold regime {Q < 40 MeV) . The corresponding references are 
listed as well. In the subsections to come we will discuss some examples of the 
kind of physics that can be studied with the various observables in the many 
reaction channels. 

5 Symmetries and their violation 

As was already mentioned in several places in this article, symmetries strongly 
restrict the allowed pattern for various observables. This leads to observable 
consequences. Naturally, it is therefore also straight forward to investigate 
the breaking of these symmetries by looking at a violation of these symmetry 
predictions. 

Probably the most prominent example of an experiment for a storage ring 
is that proposed by the TRI collaboration to be performed at the COSY 
accelerator. The goal is to do a null experiment in order to put an upper limit 
on the strength of T-odd P-even interactions via measuring Ay^^z in polarized 
proton-deuteron scattering, which should vanish if time reversal invariance 
holds [135]. For details we refer to Ref. [136]. 

5.1 Investigation of charge symmetry breaking (CSB) 

If the masses of the up and down quark were equal, the QCD Lagrangian were 
invariant under the exchange of the two quark fiavors. In reality these masses 
are not equal and also the presence of electro magnetic effects leads to small 
but measurable charge symmetry breaking effects. Note, the mass difference 
of a few MeV [137] is small compared to the typical hadronic scale of 1 GeV. 
Quantifying CSB effects therefore allows to extract information on the light 
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quark mass differences from hadronic observables. 

As should be clear from the previous paragraph, CSB is closely hnked to the 
isospin symmetry. However, does isospin symmetry demand an independence 
of the interaction under an arbitrary rotation in isospin space, a system is 
charge symmetric, if the interaction does not change under a 180 degree ro- 
tation in isospin space. Therefore isospin symmetry or charge independence 
is the stronger symmetry than charge symmetry. For an introduction into the 
subject we refer to Ref. [5]. 

The advantage of reactions with only nucleons or nuclei in the initial state is, 
that one can prepare initial states with well-defined isospin. This is in contrast 
to photon induced reactions, since a photon has both isoscalar and isovector 
components. Therefore, in the case of meson photo- or electro-production, 
CSB signals can only be observed as deviations from some expected signal. 
(As an example of this reasoning see Ref. [138].) In the case of hadron induced 
reactions on the other hand, experiments can be prepared that give a non- 
vanishing result only in the presence of CSB. This makes the unambiguous 
identification of the effect significantly easier. 

One complication that occurs if a CSB effect is to be extracted from the 
comparison of two cross sections with different charges in the final state is 
that of the proper choice of energy variable. Obviously, the quantity that 
changes most quickly close to the threshold is the phase space, so that it 
appears natural to compare two reaction channels that have the same phase 
space volume. However, due to the differences in the particle masses, this calls 
for different initial energies. In the case of a resonant production mechanism 
this might lead to effects of the same order as the effect of interest. In Ref. 
[139] this is discussed in detail for the reactions pp — > dn'^ and pn — > dn^. 

In this report we will concentrate on the implications of CSB on observables 
in NN and dd coUisions. For details of the mechanisms of CSB we refer to 
Ref. [5]. In the corresponding class of experiments the deuteron as well as 
the alpha particle play an exceptional role since as isoscalars they can act as 
isospin filters. 
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The most transparent example of a CSB reaction is 

dd ^ ax , 

where x is some arbitrary isovector. The reaction with x = tt° was recently 
measured at lUCF for the first time close to the threshold [140]. The initial 
state as well as the a are pure isoscalars. Thus the final state as an isovector 
can not be reached as long as isospin is conserved. 

In addition pn reactions can be used for producing clean signals of CSB. 
More generally, whenever a pn pair has a well-defined isospin they behave 
as identical particles. Specifically, the differential cross section needs to be 
forward-backward symmetric because nothing should change if beam and tar- 
get are interchanged. Any deviation from this symmetry is an unambiguous 
signal of CSB [5] . An experiment performed at TRIUMF recently claimed for 
the first time a non-vanishing forward-backward asymmetry in pn — > d7r° [6] . 
Note that not every forward-backward asymmetry in pn reactions stems from 
isospin violation. A counter example was given at the end of section 4.3, where 
the differential cross section for pn — > pp7r~ is discussed in detail. There, in 
contrast to the previous example, T — and T — 1 initial states interfere. 

In case of pion production in nucleon-nucleon collisions it is possible to define a 
convergent effective field theory (see section 6). Within this theory it is possible 
to relate effects of CSB in these reactions directly to the up-down quark mass 
difference [5,141]. Preliminary studies show, that the relative importance of 
different CSB mechanisms in the reaction pn — > d7r° and dd Q;7r° are very 
different [142] , and thus it should be possible to extract valuable information 
on the leading CSB operators from a combined analysis of the two reactions. 

In the arguments given all that was used was the isovector character of the 
meson produced. Thus, the same experimental signals will be seen also in 
pn d{TT^ri) [143] and dd ain^rj) [144]. In Ref. [67] also the analysing 
power was identified as a useful quantity for the extraction of the — /o 
mixing matrix element (see also discussion in Ref. [145]). The {Tr'^r))s-wave is 
interesting especially close to the KK threshold, since it should give insight 
into the nature of the light scalar mesons ao(980) and /o(980) (c.f. sec. 7.5). 
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a) b) 

Fig. 19. Illustration of different sources of charge symmetry breaking: diagram a) 
shows CSB in the production operator through tt — rj mixing and diagram b ) shows 
CSB in the propagation of the scalars. Thin solid lines denote nucleons, thick solid 
lines scalar mesons and dashed ones pseudoscalar mesons. The X indicates the 
occurrence of a CSB matrix element. 

Note that the ni] channel is the dominant decay channel of the ao, which is 
an isovector-scalar particle. It should be stressed that the charge symmetry 
breaking signal in case of the scalar mesons is significantly easier to interpret 
in comparison to the case of pion production. The reason is that the two 
scalar resonances of interest overlap and therefore the effect of CSB as it 
occurs in the propagation of the scalar mesons is enhanced compared to mixing 
in the production operator [82]. To make this statement more quantitative 
we compare the impact of /o — ao mixing in the propagation of the scalar 
mesons (Fig. 19b) to that of tt — mixing in the production operator (Fig. 
19a). We regard the latter as a typical CSB effect and thus as a reasonable 
order-of-magnitude estimate for CSB in the production operator. Observe 
that the relevant dimensionless quantity for this comparison is the mixing 
matrix element times a propagator (c.f. Fig. 19). In the production operator 
the momentum transfer — at least close to the production threshold — is given 
by t — —Mnttir, where ttir denotes the invariant mass of the meson system 
produced (or equivalently the mass of the resonance) and denotes the 
nucleon mass. Thus, the appearance of the rj propagator introduces a factor of 
about 1/t into the amplitude, since t ^ m^. On the other hand, the resonance 
propagator is given by l/{mRrR), as long as we concentrate on invariant 
masses of the outgoing meson system close to the resonance position. Here Tr 
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denotes the width of the scalar resonance. Thus we find using = 50 MeV 
[146] that the CSB in the production operator is kinematically suppressed by 
a factor of more than Vr/Mj^ ~ 1/20 as compared to CSB in the propagation 
of the scalars. In addition the mixing matrix element is enhanced in the case 
of /o — ao mixing (c.f. sec. 7.5) and therefore it should be possible to extract 
the /o — ao mixing matrix element from NN and dd induced reactions. 

6 The reaction NN NNn 

The production of pions in nucleon-nucleon collisions has a rather special role. 
First of all, it is the lowest hadronic inelasticity for the nucleon-nucleon inter- 
action and thus an important test of our understanding of the phenomenology 
of the NN interaction. Secondly, since pions are the Goldstone bosons of chi- 
ral symmetry, it is possible to study this reaction using chiral perturbation 
theory. This provides the opportunity to improve the phenomenological ap- 
proaches via matching to the chiral expansion as well as to constrain the chiral 
contact terms via resonance saturation. Last but not least, a large number of 
(un)polarized data is available (c.f. table 8). 

After a brief history, we continue this chapter with a discussion of a particular 
phenomenological model for pion production near the threshold, followed by 
a presentation of recent results from chiral perturbation theory. 

6.1 Some History 

In section 1.2 it was argued, that for the near-threshold regime the distorted 
wave born approximation is appropriate, and here we will concentrate on those 
models that work within this scheme . 

Pioneering work on pion production was done by Woodruf [147] as well as 
by Koltun and Reitan [31] in the 1960s. The diagrams included are shown in 

In the next section a further argument in favor of the distorted wave born ap- 
proximation will be given. 
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Figs. 2a and 2b where, in these early approaches, the ttN — > ttN transition 
amphtudes (denoted by T in diagram 2b) were parameterized by the scattering 
lengths. 

When the first data on the reaction pp — > pp7r° close to the threshold were 
pubhshed [43,44] , it came as a big surprise that the model of Koltun and Reitan 
[31] underestimated the data by a factor of 5-10. This is in vast contrast to the 
reaction pp pnir'^ reported in Ref. [91], where the discrepancy was less than 
a factor of 2. On the other hand, it was shown that the energy dependence of 
the total cross section can be understood from that of the NN FSI once the 
Coulomb interaction is properly included [45] (c.f. sec. 2). 

Niskanen investigated whether the inclusion of the Delta isobar, as well as 
keeping the rather strong on-shell energy dependence of the nN interaction, 
could help to improve the theoretical results for neutral pion production [148]. 
Although these improvements lead to some enhancement, the cross section was 
still missed by more than a factor of 3. 

The first publication that reported a quantitative understanding of the pp — > 
ppn^ data was that by Lee and Riska [149] and later confirmed by Horowitz et 
al. [150] , where it was demonstrated that short range mechanisms as depicted 
in Fig. 2c, can give a sizable contribution. However, shortly after this discovery 
Hernandez and Oset demonstrated, using various parameterizations for the 
ttN — > ttN transition amplitude and qualitatively reproducing earlier work 
by Hachenberg and Pirner [151] that the strong off-shell dependence of that 
amplitude can also be sufficient to remove the discrepancy between the Koltun 
and Reitan model and the data. Gedalin et al. came to the same conclusion 
within a relativistic one boson exchange model [152]. In Ref. [153] the nN 
amplitude needed as input for the evaluation of diagram 2b was extracted from 
a microscopic model. Also there a significant although smaller contribution 
from the pion rescattering was found. Thus, in this model still some additional 
short range mechanism is needed. 

In the succeeding years many theoretical works presented calculations for the 
pp — > pp7r° cross section. In Refs. [154,155] covariant one boson exchange mod- 
els were used in combination with an approximate treatment of the nucleon- 
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nucleon interaction. Both models turned out to be dominated by heavy meson 
exchanges and thus give further support to the picture proposed in Ref. [149]. 
However, in Ref. [156] the way that the anti nucleons were treated in Refs. 
[149,150,154] was heavily criticized: the authors argued that the anti nucleon 
contributions get significantly suppressed once they are included non pertur- 
batively. It is interesting to note, that also in Bremsstrahlung the contribution 
from anti-nucleons in a non-perturbative treatment is significantly reduced 
compared to a perturbative inclusion [157,158]. Additional short range con- 
tributions were also suggested, namely the p — uj meson exchange current 
[159], resonance contributions [160] and loops that contain resonances [161], 
all those, however, turned out to be smaller compared to the heavy meson 
exchanges and the off-shell pion rescattering, respectively. 

At that time the hope was that chiral perturbation theory might resolve the 
true ratio of rescattering and short range contributions. It came as a big sur- 
prise, however, that the first results for the reaction pp — > pp7r°[162,163] found 
a rescattering contribution that interfered destructively with the direct contri- 
bution (diagram a in Fig. 2), making the discrepancy with the data even more 
severe. In addition, the same isoscalar rescattering amphtude also worsened 
the discrepancy in the %^ channel [164,165]. Some authors interpreted this 
finding as a proof for the failure of chiral perturbation theory in these large 
momentum transfer reactions [155,166]. Only recently was it demonstrated, 
that it is possible to appropriately modify the chiral expansion in order to 
make it capable of analysing meson production in nucleon-nucleon collisions. 
We will report on those studies in sec. 6.3 that in the future will certainly 
prove useful for improving the phenomenological approaches (c.f. subsec. 6.5). 

Before we close this section a few remarks on the off-shell t^N amphtude are 

necessary. For this purpose we write the relevant piece of the ttA^ interaction 
T-matrix in the following form 

T^N = ^ ATtTT • 'kN + t^-^ ATtr • (tt X tt) A^ , (70) 

where t^'^^ {t^~^) denote the isoscalar (isovector) component. Note that it is 
only the former that can contribute to the reaction pp — > pp7r° , for the 

This is only true if we do not include the A isobar explicitly, as will be discussed 
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isospin structure of the latter changes the total isospin of the two nucleon 
system. 

As long as we neglect the distortions due to the final and initial state interac- 
tions, what is relevant for the discussion in this paragraph is the half off-shell 
t:N amplitude. We thus may write t = t{s, k^), where s denotes the invariant 
energy of the ttN system and is the square of the four momentum of the 
incoming pion. At the threshold for elastic ttN scattering (s = (m^r + M^Y, 
k'^ — rri^) we may write 

(so, ml) := Y (l + ^) (ai + 203) = (-0.05 ±0.01) m;^ , 
t^-)[s,,ml):^^{l + ^){a,-a,) = (1.32 ± 0.02) m;^ (71) 

where the 02/ denote the scattering lengths in the corresponding isospin chan- 
nels /. The corresponding values were extracted from data on -R^d atoms in 
Ref. [167]. Note that the dominance of the isovector interaction is a conse- 
quence of the chiral symmetry: the leading isoscalar rescattering is suppressed 
by a factor m.,^/Mj^ compared to the leading iso-vector contribution — the so 
called Weinberg-Tomozawa term [168,169]. However, it is still remarkable, 
that also the higher order chiral corrections are small leaving a value consis- 
tent with zero for the isoscalar scattering length. A detailed study showed, 
that this smallness is a consequence of a very efficient cancellation of sev- 
eral individually large terms that are accompanied with different kinematical 
factors [170]. To be concrete: to order O(p^) one finds 

= A [-2mlc, + g;A;o [c., + , (72) 

where k and q' denote the four momentum of the initial and final pion respec- 
tively. The values for the various q are given in table 9. For on-shell scattering 
at the threshold {q' — k — (m^,0)) one gets t^'^\sQ,ml) — —0.2Am~^ using 
the values of Ref. [171] . Please note, that the linear combination of the q ap- 
pearing above turns out to be an order of magnitude smaller than the individ- 

in the next section. 

^■"^Note, this value is inconsistent with the empirical value given in Eq. (71). To 
come to a consistent value one has to go to one loop order as discussed in Ref. [172]. 
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ual values. As a consequence, the on-shell isoscalar amplitude shows a rather 
strong energy dependence above threshold. It should not then come as a sur- 
prise, that the transition amplitude corresponding to Eq. (70), when evaluated 
in the kinematics relevant for pion production in NN collisions, within chiral 
perturbation theory turns out to be rather large numerically [162,163]. For the 
non-covariant expression given in Eq. (72) this translates into q' — {111^^, 0) and 
k = (m7r/2, k), leading to t'^'^\{m^ + M^f, -MnttIt,) = 0.5 This is why 
the first calculations using chiral perturbation theory found a big effect from 
pion rescattering — unfortunately increasing the discrepancy with the data. In 
Ref. [173] the tree level chiral perturbation theory calculations where repeated 
using a different prescription for the energy of the exchange pion . The au- 
thors found agreement with the data, but with a sign of the full amplitude 
opposite to the one of the direct term. 

One year earlier it was shown, that within phenomenological approaches the 
isoscalar transition amplitude evaluated in off-shell kinematics is also signif- 
icantly different from its on-shell value. For example, in the Jiilich meson 
exchange model it is the contribution from the iterated p i-channel exchange 
and the a exchange that are individually large but basically cancel in threshold 
kinematics in the isoscalar channel [153]. This cancellation gets weaker away 
from the threshold point. This rescattering contribution, however, turned out 
to interfere constructively with the direct term. 

Since chiral perturbation theory as the effective field theory for low energy 
strong interactions is believed to be the appropriate tool to study pion reac- 
tions close to threshold, it seemed at this stage as if there were a severe prob- 
lem with the phenomenology. However, as was shown in section 4.3, there are 
observables that are sensitive to the sign of the s-wave pp ppn^ amplitude — 
relative to a p-wave amplitude that is believed to be under control — and the 
experimental results [74,75,77,78] agree with the sign as given by the phe- 
nomenological model (c.f. Fig. 15) . 

Does this mean that chiral perturbation theory is wrong or not applicable? No. 
As we will discuss in the subsequent sections, it was demonstrated recently 

This prescription was later criticized in Ref. [174]. 
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that the chiral counting scheme needs to be modified in the case of large 
momentum transfer reactions. No complete calculation has been carried out 
up to now, but intermediate results look promising for a consistent picture to 
emerge in the years to come. The insights gained so far from the effective field 
theory studies call also for a modification of the phenomenological treatment. 
This will be discussed in detail in section 6.5. 

6.2 Phenomenological approaches 

As stressed in the previous section, the number of phenomenological models 
for pion production is large. For dcfinitcncss in this section we will focus on one 
particular model, namely that presented in Refs. [71,73,153], mainly because 
it incorporates most of the mechanisms proposed in the literature for pion 
production in nucleon-nucleon collisions, its ingredients are consistent with the 
data on vrA^ scattering, and it is the only model so far whose results have been 
compared to the polarization data recently measured at lUCF [47,94,108]. 

The model is the first attempt to treat consistently the A^A^ as well as the ttN 
interaction for meson production reactions close to the threshold: both were 
taken from microscopic models (described in Refs. [175] and [176] for the A^A^ 
and the ttN interaction, respectively). These were constructed from the same 
effective Lagrangians consistent with the symmetries of the strong interaction 
and are solutions of a Lippman-Schwinger equation based on time-ordered 
perturbation theory. Although not all parameters and approximations used in 
the two systems are the same, this model should still be viewed as a benchmark 
calculation for pion production in A^A^ collisions within the distorted wave 
Born approximation. We will start this section with a description of the various 
ingredients of the model and then present some results. 

6.2.1 The NN interaction 

A typical example of a so-called reahstic model for A^A^ scattering is the Bonn 
potential [26] . The model used for the A" A" distortions in the final and initial 
states is based on this model, where a pseudo-potential is constructed based 
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Fig. 20. Comparison of the model predictions of Ref. [73] to the data. The references 
for the experimental data can be found in table 8. The solid lines show the results 
of the full model; the dashed line shows the results without the A contributions. 



on the t-channel exchanges of all established mesons below one GeV in mass 
between both nucleons and Delta isobars. 

The interaction amongst the various djTiamical fields in the model derived 
from the following Lagrange densities: 



Cnnit = ^^757^"^ ■ d'^Trip (73) 

Cnnp = 9NNph,T ■ p^^ + j^i'^.uT ■ id^p' - d^p^i^ (74) 

^NNu, = QNNui^lfi^^i' ( 75 ) 

i^NNa = gNNai^C^i' (76) 

^^NNao = QNNao^T ■ flo^ (77) 

CNi^.= ^-^i'f -d^TTi,^ + h.c. (78) 

Cnap = t^i'i.i.f ■ idf-p" - a>^)^, + h.c. . (79) 
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E,,JMeV] E^,[MeV] 

Fig. 21. The NN phase shifts for the model of Ref. [175] in the energy range relevant 
for pion production. The pion production threshold is at Elab = 286 MeV. The 
experimental data are from Refs. [177] (triangles) and [61] (circles). 

Note: The particle called 5 in the original Bonn publication [26] is nowadays 
called oq. The operator T as well as the fields are defined in Ref. [26]. Note 
that there is no tensor coupling for the ujNN vertex given for the fit to the 
elastic NN scattering data did not need any such coupling. 

There is a difference between the nucleon-nucleon model used [175] and the 
Bonn potential [26]. The original Bonn model has an energy dependent in- 
teraction, for it keeps the full meson retardation in the intermediate state. 
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This, however, leads to technical problems, once the model is to be evalu- 
ated above the pion production threshold due to the occurrence of three body 
singularities. In Ref. [24] those singularities were handled by solving the dy- 
namical equations in the complex plane. Unfortunately, this method is not 
useful for the application in a distorted wave Born approximation. Instead we 
used a model based on the so-called folded diagram formalism developed in 
Ref. [178]. This formalism, worked out to infinite order, is fully equivalent to 
time-ordered perturbation theory. When truncated at low order, however, it 
leads to energy independent potentials that can formally be evaluated even 
above the pion production threshold. In addition, the model of Ref. [175] is 
constructed as a coupled-channel model including the A'^A'' as well as the A^A 
and AA channels. This enables us to treat the A isobar on equal footing with 
the nucleons. 

The resulting phase shifts are shown in Fig. 21. Note that the model param- 
eters were adjusted to the phase shifts below the pion production threshold 
only, which is located at E^ab = 286 MeV. Fig. 21 thus clearly illustrates 
that using this model for the A'^A'^ interaction is indeed justified and we may 
conclude that — at least up to laboratory energies of 600 MeV-the A^A^ phe- 
nomenology is well understood. 

6.2.2 The nN model 

The ttN interaction that enters in the pion rescattering diagrams can be taken 
from a meson exchange model as well [179] . This allows a consistent treatment 
of the meson and nucleon dynamics. It should be stressed that this is the 
precondition for comparing the results of the phenomenological model to those 
of chiral perturbation theory, as we will do below. In addition, since we also 
want to include the rescattering diagram in partial waves higher than the s- 
wave, after a fit to the ttN data the pole contributions (nucleon and Delta) 
need to be removed from the amplitudes in order to avoid double counting 
with the direct production. This is possible only within a microscopic model. 

The main features of the nN model of Ref. [179] are that it is based on an ef- 
fective Lagrangian consistent with chiral symmetry to leading order and that 
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Fig. 22. Contributions to the potential of the model of Ref. [179]. 

the i-channel exchanges in the isovector (p) and isoscalar (a) channel are 
constructed from dispersion integrals . For details on how the i-channel ex- 
changes are included we refer to Ref. [179]. The diagrams that enter the poten- 
tial are displayed in Fig. 22. This potential is then unitarized with a relativistic 
Lippmann-Schwinger equation — in complete analogy to the nucleon-nucleon 
interaction. 



Within this model, at tree level the isoscalar and isovector nN interaction are 
give by the corresponding i-channel exchanges. In the latter case the unitariza- 
tion does not have a big influence in the near-threshold regime, and thus also 
the isovector scattering length is governed by the tree level p-exchange. The 
famous KFSR relation [181,182], that relates the couplings of the p-meson to 
pions and nucleons to the coupling strength of the Weinberg-Tomozawa term, 
is a consequence of this. On the other hand, for the isoscalar nN interaction 
the effects of the unitarization are large and lead to an almost complete can- 

There are ambiguities in how to extrapolate the results of the dispersion integrals 
to off-shell kinematics. This issue is discussed in detail in Ref. [180]. 
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Fig. 23. The ttN phase shifts for the model of Ref. [179] in the energy range rele- 
vant for pion production. The experimental data are from Refs. [183, 184] ■ Note the 
different scales of the various panels. 



cellation of the isoscalar potential with the iterated p-exchange in the near 
threshold regime. As one moves away from the threshold value this cancella- 
tion gets weaker leading to the strong variation of the off-shell isoscalar vrA^ 
T-matrix mentioned at the end of sec. 6.1. 

In Figure 23 the results for the model of Ref. [179] are compared to the data 
of Refs. [183,184]. As one can see the model describes the data well, especially 
in the most relevant partial waves: Su, S^i and P33. 
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6.2.3 Additional short range contributions and model parameters 

As was stressed above, a large class of additional mechanisms was suggested in 
the literature to contribute significantly to pion production in nucleon-nucleon 
collisions. Since they are all of rather short range and mainly influence the 
production of s-wave pions, in this work only a single diagram was included 
(heavy meson exchange through the uj — Fig. 2c) to parameterize these various 
effects. Consequently, the strength of this contribution was adjusted to repro- 
duce the total cross section of the reaction pp — > ppn^ close to the production 
threshold. The short-range contributions turn out to contribute about 20% to 
the amplitude. After this is done, all parameters of the model are fixed. 

6.2.4 Results 

The results of the model presented have already appeared several times in 
this report (Figs. 13,14,15,17,20,24) — mainly for illustrative purposes — and 
they are discussed in detail in Refs. [73,71]. 

Overall the model is rather successful in describing the data, given that only 
one parameter was adjusted to the total cross section for low energy neutral 
pion production (see sec. 6.2.3). One important finding is that the sign of the 
s-wave neutral pion production seems to be in accord with experiment, as is 
illustrated in Fig. 15 (c.f. corresponding discussion in sec. 4.3), in contrast to 
the early calculations using chiral perturbation theory. As we will see in sec. 
6.3, today we know that those early calculations using effective field theory 
were incomplete. 

The most striking differences appear, however, for double polarization observ- 
ables in the neutral pion channel, as shown in Figs. 13 and 14. As a general 
pattern the amplitudes seem to be of the right order of magnitude, but show 
a wrong interference pattern. To actually allow a detailed comparison of the 
model and data a partial wave decomposition of both is necessary. Work in 
this direction is under way. 

It is striking that for charged pion production the pattern is very different, for 
here almost all observables are described satisfactorily. In Fig. 24 the results of 
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the model for charged and neutral pion production are compared to the data 
for a few observables. In contrast to the neutral channel, the charged pion 
production is completely dominated by two transitions, namely ^Pi — >^ SiS, 
which is dominated by the isovector pion rescattering, and — >^ Sip, which 
governs the cross section especially in the regime of the Delta resonance. The 
prominence of the Delta resonance is a consequence of the strong transition 
^D2{NN) — >^ 5*2 (A^A) that even shows up as a bump in the NN phase shifts 
(see Fig. 21). This effect was first observed long ago and is well known (see, 
e.g., discussion in Ref. [185]). As a consequence, the NN — > A'"A transition 
potential should be rather well constrained by the NN scattering data and is 
not the case for all the many transitions relevant in case of the neutral pion 
production, where the Sp final state is not allowed due to selection rules (see 
sec. 1.6). One might therefore hope to learn more about the AN interaction 
from the pion production data. 

One can also ask how well we know the production operator. Fortunately, at 
the pion production threshold it is still possible to analyze meson production 
in A^A^ collisions within effective field theory. This analysis will give deeper 
insight into the production dynamics, as will be explained in the following 
section. The two approaches are then compared in sec. 6.5. 

6.3 Chiral perturbation theory 

The phenomenological approaches, such as the one described in the previous 
section, lack a systematic expansion. Thus it is neither possible to estimate 
the associated model uncertainties nor to systematically improve the models. 
On the other hand, for various meson production reactions the phenomenolog- 
ical approaches proved to be quite successful. One might therefore hope that 
effective field theories will give insights into why the phenomenology works, 
as was first stressed by Weinberg [186]. 

A first attempt to construct — model independently — the transition amplitude 
NN — > A'"A'"7r was carried out almost 40 years ago [187-189]. The authors 
tried to relate what is known about nucleon-nucleon scattering to the pro- 
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Fig. 24. Comparison of the model predictions to the data taken from Ref. [47] 
(pp — > ppTT°), Ref. [108] (pp dn^ ) and Ref. [94] (pp prnr^ ). 



duction amplitude via the low energy theorem of Adler and Dothan [190], 
which is a generalization of the famous Low theorem for conserved currents 
[191] to partially conserved currents via the PCAC relation (see, e.g., Ref. 
[192]). However, it soon was realized that the extrapolation from the chiral 
limit (mTT = 0) to the physical point can change the hierarchy of different 
diagrams. The reason for the non-applicability of soft radiation theorems to 
meson production close to the threshold is easy to see: a necessary condition 
for the soft radiation theorem to be applicable is that the energy emitted is 
significantly smaller than the typical energy scale, that characterizes variations 
of the nuclear wave function. In the near-threshold regime, however, the scale 
of variation is set by the inverse of the NN scattering lengths — thus a soft 
radiation theorem of the type of Low or that of Adler and Dothan could only 
be applicable to meson production if m^r ^ l/(M7va^). In reality, however, 
the pion mass exceeds the energy scale introduced by inverse NN scattering 
length by more than two orders of magnitude. The range of applicability of 
soft radiation theorems is discussed in Ref. [193] in a different context. 
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More recent analyses, however, show that in the case of pion production it 
is indeed possible to define a convergent effective field theory that allows a 
systematic study of the structure of the production operator. As we will show, 
most of the diagrams included in the Koltun and Reitan model [31] are indeed 
the leading operators in pion production. In addition this study will show 

• that the use of the distorted wave born approximation is justified 

• why neutral pion production is the more problematic case 

• the importance of loop contributions 

• that there is a close connection between pion production in nucleon-nucleon 
collisions and the three nucleon problem. 

The problem with strong interaction phenomena is their non-perturbative 
nature with respect to the coupling constants. To construct a controlled ex- 
pansion it is necessary, to identify a small expansion parameter. In general 
this is only possible for a limited energy range. The conditio sine qua non 
for constructing an effective field theory for any system is the separation of 
scales characteristic for the system. Once the — in this context light — scales 
are identified, one treats them dynamically, while all the dynamics that are 
controlled by the heavy scales are absorbed in contact interactions. As long 
as the relevant external momenta and energies are such that structures of the 
size of the inverse of the heavy scale can not be resolved, this procedure should 
always work. Weinberg [194] as well as Gasser and Leutwyler [195] have shown 
that this general idea works even when loops need to be included. 

In case of low energy pion physics it is the chiral symmetry that provides 
both preconditions for the construction of an effective field theory, in that 
it forces not only the mass of the pion m^, as the Goldstone boson of the 
chiral symmetry breaking, to be low, but also the interactions to be weak, 
for the pion needs to be free of interactions in the chiral limit for vanishing 
momenta. The corresponding effective field theory is called chiral perturbation 
theory (xPT) and was successfully applied to meson-meson [196] scattering. 
Treating baryons as heavy allows straightforward extension of the scheme to 
meson-baryon [197] as well as baryon-baryon [198-201] systems. In all these 
references the expansion parameter used wasp/A^ ~ p/ {4:71 f^r) ~ p/Mjy, where 
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denotes the chiral symmetry breaking scale, /jr the pion decay constant and 
Mjv the nucleon mass. Recently it was shown that also the Delta isobar can 
be included consistently in the effective field theory [202] . The authors treated 
the new scale, namely the Delta nucleon mass sphtting A = Ma — Mjv, to be 
of the order of p that is taken to be of the order of m^. 

An additional new scale occurs for meson production in nucleon-nucleon colli- 
sions, namely the initial momentum pi ~ y/rr^J^. Note that, although larger 
than the pion mass, this momentum is still smaller than the chiral symme- 
try breaking scale and thus the expansion should still converge, but slowly. A 
priori there are now two options to construct an effective field theory for pion 
production. The first option, called Weinberg scheme in what follows, treats 
all hght scales to be of order of m^. Thus, in this case, there is one expansion 
parameter, namely xw — rriT^/Mj^. The other option is to expand in two scales 
simultaneously, namely rriT^ and pj. In this case the expansion parameter is 



This scheme was advocated in Refs. [163,203] and apphed in Ref. [204]. The 
two additional scales, namely A and ttItt, are identified with 

^ Pi A Pi 2 /Qn^ 

where the former assignment was made due to the numerical similarity of the 
two numbers (A = 2.1mT, and Pi = 2.6m^). Only explicit calculations can 
reveal which one is the more appropriate approach. 

Within the Weinberg counting scheme, tree level calculations were performed 
for s-wave pion production in the reactions pp ppn^ [162,163,173] as well as 
pp pniT^ [164,165]. In addition, complete calculations to next-to-next-to- 
leading order (NNLO), where in the Weinberg scheme for the first time loops 
appear, are available for pp pp-K^ [205,206]. The authors found that some 
of the NNLO contributions exceeded significantly the next-to-leading (NLO) 
terms leading them to the conclusion that the chiral expansion converges only 
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Note that A stays finite in the chiral limit, whereas both pi as well as ttZtt vanish. 
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slowly, if at all. This point was further stressed in Refs. [207,155], however, it 
was shown recently that as soon as the scale induced by the initial momentum 
is taken into account properly (expansion in x iiot in Xw), the series 
indeed converges [203,204]. For illustration, in this section we compare the 
order assignment of the Weinberg scheme to that of the modified scheme. 

In appendix E the relevant counting rules of the new scheme are presented as 
well as justified via application to a representative example. Especially, it is 
not clear a priori what scale to assign to the zeroth component Iq of the four 
dimensional integration volume dH as it occurs in covariant loops. After all, 
the typical energy of the system is given by uIt^, but the momentum by pi. As 
is shown in the appendix by matching a covariant analysis to one carried out 
in time-ordered perturbation theory, in loops Iq ^ Pi- This assignment was 
also confirmed in explicit calculations [204] . 

The starting point is an appropriate Lagrangian density, constructed to be 
consistent with the symmetries of the underlying more fundamental theory 
(in this case QCD) and ordered according to a particular counting scheme. 
Omitting terms that do not contribute to the order we will be considering 
here, we therefore have for the leading order Lagrangian [208,199,197] 



2 M 2 2/2 



+A^^[iao 



1 



4/^ 



(tt X ■k)]N 



+^nW ■ a ■ ( Vtt + ^7r(7r • Vtt) ) N 

+^i[ido - A]*A + ^[N\T ■ S ■ V7r)*A + h.c] + 



(81) 



The expressions for interactions with more than two pions depend on the 
interpolating field used. The choice made here was the so called sigma gauge — 
c.f. Appendix A of Ref. [197], where also the corresponding vertex functions 
are given explicitly . 

For the next-to-leading order Lagrangian we get 



2^ As usual, all observables are independent of the choice of the pion field. 
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2mjv 

+-T^^(iA^V • (tt X Vtt) • VA^ + /i.e.) 

+^N^[ic2 + C3 - -^)^^ - C3(V7r)2 - 2cim^7r^ 

~o('^4 + )£ijk£abcCrkrcdi7radj7rb]N 

2 4mAr 

[iA^V • 7r<T • VA^ + h.c] 
[iN^'T-ivS-V-^A + h.c.] 



~N^{T-a-V7r)NN^N 

Jn 



d2 



SijkSabcdi'KaN^CyjnN N'^GkTcN ^ ■ ■ ■ , (82) 



where /^r denotes the pion decay constant in the chiral hmit, qa is the axial- 
vector couphng of the nucleon, Ha is the ANn couphng, and 5* and T are the 
transition spin and isospin matrices, normahzed such that 

^i^J ^ 3 - i^ijkCTk) , (83) 
TiT] = ^{2dij-ieijkTk) . (84) 

These definitions are in fine with the ones introduced in sec. 6.2.1. The dots 
symbohze that what is shown are only those terms that are relevant for the cal- 
culations presented. As demanded by the heavy baryon formalism, the baryon 
fields A^ and are the velocity-projected pieces of the relativistic fields ap- 
pearing in the interactions discussed in section 6.2.1; e.g. A'" = 1/2(1-|- i/jip, 
where v'^ denotes the nucleon 4-velocity. 

The terms in the Lagrangians given are ordered according to the conventional 
counting {p ~ m^). A reordering on the basis of the new scheme does not seem 
appropriate, for what order is to be assigned to the energies and momenta 
occurring depends on the topology of a particular diagram (see also Appendix 
E). 

The constants q can be extracted from a fit to elastic nN scattering. This was 
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Fig. 25. Illustration of resonance saturation. 



Table 9 

The various low energy constants Ci in units of GeV^^. The first two columns give 
the values from the original references (left column: tree level calculation of Ref. 
[171]; right column: one loop calculation of Ref. [172]), whereas the last two columns 
give those values reduced by the Delta contribution as described in the text. It is those 
numbers that where used in the calculations presented here. 



i 


„tree 
H 


loop 






1 


-0.64 


-Q.93 


-0.64 


-0.93 


2 


1.78 


3.34 


0.92 


0.64 


3 


-3.90 


-5.29 


-1.20 


-2.59 


4 


2.25 


3.63 


0.9 


2.28 



done in a series of papers with successively improved methods [171,172,209- 
212]. However, here we will focus on the values extracted in Refs. [171,172] for 
it is those that were used in the calculations for pion production in nucleon- 
nucleon collisions. In the former work the q were extracted at tree level and 
in the latter to one loop. The corresponding values are given in table 9. In 
both papers the Delta isobar was not considered as exphcit degree of freedom. 
In an effective field theory the low energy constants appearing depend on the 
dynamical content of the theory. Thus, in a theory without explicit Deltas, 
their effect is absorbed in the low energy constants [172]. This is illustrated 
graphically in Fig. 25. Thus we need to subtract the Delta contribution from 
the values given in the first columns of table 9. Analytical results for those 
contributions are given in Refs. [213,172] : 



Co — Co — 2c/i 



/)2 



9 (Ma - M, 



2.7 GeV-^ 



N 



(85) 



To match the results of the two references the large Nc value has to be used for 
Ha = S9A/V2 ~ 2.7. 
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a) 



b) 



Fig. 26. Illustration of the role of the ANtt contact term in NN NNn and three 
nucleon scattering. Solid lines denote nucleons, dashed lines denote pions. 

Thus, the only undetermined parameters in the interaction Lagrangian are di 
and ^2- Since they are the strength parameters of 4-nucleon contact interac- 
tions that do not contain any derivatives of the nucleon fields, they only con- 
tribute to those amphtudes that have an NN S-wme to NN S-w&ve plus pion 
p-wave transition. This automatically excludes a transition from an isospin 
triplet to an isospin triplet state, for this would demand to go from ^5*0 to ^5*0 
accompanied by a p-wave pion, which is forbidden by conservation of total 
angular momentum. Thus, only two transitions are possible: T — ^ T — 1, 
as it can be studied in pn — > ppTr~ , and T=l— >T = 0, asit can be studied 
in pp — > pmr'^. Even more importantly, in both channels the di appear with 
the same linear combination d which is completely fixed by the corresponding 
isospin factors: 



where we have introduced the dimensionless parameter 5 . In order for the 
counting scheme to work, S — 0{1) has to hold. As we will show below, this 
order of magnitude is indeed consistent with the data. Please note, the four- 
nucleon-pion contact term with p-wave discussed here was to some extend 
also investigated in Ref. [214], however, within a different expansion scheme. 
Thus direct comparison is not possible. 

The parameter d is very interesting, for it is at the same time the leading short- 
range-long-range contribution^^ to the three nucleon force, as illustrated in 

^^In this context, pion exchanges are called long-ranged, whereas any exchange of 
heavier mesons — absorbed in the contact terms — are called short-ranged. 
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Fig. 26. Naturally, it can also be fixed from pd scattering data directly. This 
was done in Ref. [215]. We come back to this point below. 

The effect of the large scale on the vertices is best illustrated with an example. 
In £(0) the so^alled Weinberg-Tomozawa (WT) term 1/(4/2) ATtr • (tt ^ 
appears. The corresponding recoil correction l/{%Mi^fl){iN^T ■ (tt x Vtt) • 
VA^ + h.c.) appears in Therefore the vertex function derived from the 
WT term is proportional to {qq + ko)/ f^, where qo (ko) denote the zero-th 
component of the 4-vectors for the outgoing (incoming) pion. If the WT term 
appears as the ttN vertex in the rescattering contribution (c.f. Fig. 2b), in 
threshold kinematics ko = m^/2 and go = ^tt- Also in threshold kinematics, 
where the 3-momentum transfer equals the initial momentum, we find for the 
contribution from the recoil term / {2Mn f"^). Obviously, since = MNnin, 
the contribution from the WT term and from its recoil term are of the same 
order. This changes if the WT term appears inside a loop, for then the scale 
for ko is also set by p — in this case the recoil term is suppressed by one chiral 
order compared to the WT term itself. The assignments made are confirmed 
by explicit calculation [204] as well as by a toy model study [174]. 

Since we know now the interactions of pions and nucleons we can investigate 
the relevance of sublcading loops, where we mean loops that can be constructed 
from a low order diagram by inserting an additional pion line. Obviously this 
procedure introduces at least one nNN vertex ~ p/ fn, a pion propagator 
~ l/p"^, an integral measure p'^/(47r)2 and either an additional ttNN vertex 
together with two nucleon propagators, or one additional nucleon propagator 
together with a factor 1/ fn- To get the leading piece of the loops in the lat- 
ter case, the integration over the energy variable has to pick the large scale 
and thus we are to count the nucleon propagator as 1/p, leading to an overall 
suppression factor for this additional loop of p^/M^. In the former case, how- 
ever, a topology is possible that contains a two-nucleon cut. This unitarity cut 
leads to an enhancement of that intermediate state, for it pulls a large scale 
(Mjv) into the numerator. Based on this observation Weinberg established a 
counting scheme for NN scattering that strongly differs from that in tttt as 
well as ttN scattering [186]. The corresponding factor introduced by this loop 
is p/Mj^ (and in addition typically comes with a factor of tt). As was stressed 
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Table 10 

Comparison of the corresponding chiral order in the Weinberg scheme (p ~ m.,^) 
and the new counting scheme (p ~ \Jm-,^MN ) for several nucleonic contributions for 
p-wave pion production. Subleading vertices are marked as 0. Here q denotes the 
external pion momentum. For simplicity we assume the outgoing nucleons at rest. 

p-wave diagrams 



(nucleons only) 




b) 




g) 



h) 



scale 


p ~ m^r 


p ~ 




LO 


LO 


qniT^/p^ 


LO 


N^LO 


q/Mn 


NLO 


N^LO 




N^LO 


N^LO 



by Weinberg, this suppression is compensated by the size of the corresponding 
low energy constants of the two-nucleon interaction and therefore all diagrams 
that can be cut by crossing a two nucleon line only are called reducible and 
the initial as well as final state NN interaction is summed to all orders. This 
is what is called distorted wave Born approximation. There is one exception 
to this rule: namely when looking at the two-nucleon intermediate state close 
to the pion production vertex in diagram a of Fig. 2. Either the incoming or 
the outgoing nucleon needs to be off-shell and thus this intermediate state 
does not allow for a two nucleon unitarity cut. Therefore those two-nucleon 
intermediate states are classified as irreducible. 
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There is a special class of pion contributions not yet discussed, namely those 
that contain radiative pions (on-shell pions in intermediate states). We re- 
strict ourselves to a kinematic regime close to the production threshold. As 
soon as an intermediate pion goes on-shell, the typical momentum in the cor- 
responding loop automatically needs to be of order of the outgoing momenta. 
This leads to an effective suppression of radiative pions; e.g., for s-wave pions 
the effects of pion retardation become relevant at N^LO. 

Note that each loop necessarily contributes at many orders simultaneously. 
The reason for this is that the two scales inherent to the pion production prob- 
lem can be combined to a dimensionless number smaller than one: m^/pj = x 
. However, what can always be done on very general grounds is to assign the 
minimal order at which a diagram can start to contribute and this is sufficient 
for an efficient use of the effective field theory. 

The next step is the consistent inclusion of the nuclear wave functions. How- 
ever, the NN potentials constructed consistently with chiral perturbation the- 
ory [198-201] are not applicable at the pion production threshold. Therefore 
we use the so-called hybrid approach originally introduced by Weinberg [216], 
where we convolute the production operator, constructed within chiral per- 
turbation theory, with a phenomenological NN — NA wavefunction. We use 
the CCF model described in the previous section [175]. 

Let us start with a closer look at the production of p-wave pions, for those 
turned out easier to handle than s-wave pions. The reason for this pattern lies 
in the nature of pions as Goldstone bosons of the chiral symmetry: since in the 
chiral limit for vanishing momenta the interaction of pions with matter has 
to vanish, the coupling of pions naturally occurs in the company of either a 
derivative or an even power of the pion mass . As a consequence, the leading 
piece of the nNN vertex is of p-wave type, whereas the corresponding s- 
wave piece is suppressed as uj^^/Mn, where uji^ denotes the pion energy. Note 

Only even powers of the pion mass are allowed to occur in the interaction, since, 
due to the Gell-Mann-Oakes-Renner relation [217], oc niq, where ruq is the 
current quark mass and in the interaction no terms non-analytic in the quark masses 
are allowed. 
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Fig. 27. Comparison of the predictions from effective field theory at LO (dashed line) 
and NLO (with the ttN parameters from a NLO (dotted line) and an NNLO (solid 
line) analysis) with the '^ai cross section for pp ppir^ [4?]. 

that also in the case of neutral pion photoproduction the s-wave amplitude is 
dominated by pion loops, whereas the p-wave amplitude is dominated by tree 
level diagrams [218,197]. 

The corresponding diagrams are shown in table 10 up to N^LO. So far in 
the literature calculations have been carried out only up to N^LO [203]. An 
important test of the approach is to show its convergence. For that we need 
an observable to which s-wave pions do not contribute. Such an observable is 
given by the spin cross section ^cxi recently measured at lUCF for the reaction 
pp ppiT^ [47] (c.f. table 7 in section 4.2). The parameter-free prediction 
of chiral perturbation theory compared to the data is shown in Fig. 27. As 
one can see, the total amplitude is clearly dominated by the leading order 
suggesting a convenient rate of convergence for the series. In addition, the 
prediction agrees with the data. 

Thus we are now prepared to extract the parameter d from data on the reaction 
pp pniT'^ [219]. As it was argued above, only the amplitude corresponding 
to the transition ^Sq — >^ Sip, called oq, is influenced by the corresponding 
contact interactions. The results of the chiral perturbation theory calculations 
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Fig. 28. oo ofpp npn^ in chiral perturbation theory. The different lines correspond 
to values of the parameter related to the three-nucleon force: 8 = 1 (long dashed line). 
(5 = (dot-dashed line), 5 = —0.2 (solid line), and 6 = —1 (short dashed line). Data 
are from Ref. [219]. 

are shown in Fig. 28. The figure shows four curves for different values of the 
parameter 5 defined in Eq. (86), namely the result for 5 = (dot-dashed line), 
for 5 = —0.2 — the authors of Ref. [220] claim this value to yield an impor- 
tant contribution to Ay in Nd scattering at energies of a few MeV — as well 
as the results we get when 6 is varied within its natural range 6 = +1 and 
6 = —1 shown as the long-dashed and the short-dashed curve, respectively. 
Thus we find that the results for are indeed rather sensitive to the strength 
of the contact interaction. This might be surprising at first glance, since we 
are talking about a subleading operator, however it turned out that the lead- 
ing (diagrams a) and b) in table 10) and subleading contributions (the same 
diagrams with a A intermediate state) largely cancel. Thus, to draw solid con- 
clusions the p-wave calculations should be improved by one chiral order, the 
corresponding diagrams of which are shown in table 10g)-j). 

Please note that the rescattering contribution involving C4 that occurs at the 
same order as the d contribution turned out to be sensitive to the regulator 
used for the evaluation of the convolution integral with the nuclear wave- 
function. This cutoff dependence can be absorbed in d, which in turn is now 
cutoff-dependent. Therefore, in order to compare the results for d from the 
pion production reaction to those extracted from the three-nucleon system 

The calculation of Ref. [220] suffers from numerical problems. 
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[215] a consistent calculation that is not possible at present has to be per- 
formed. Note that the d parameter as fixed in Ref. [215] (there it is called E) 
also turned out to be sensitive to the regulator. On the long run, however, a 
consistent description of pion production and three-nucleon scattering should 
be a rather stringent test of chiral effective field theories at low and intermedi- 
ate energies. As should be clear from the discussion above, to yield values for 
the low energy constant 5 that are compatible both calculations pd scatter- 
ing as well as 7r-production have to be performed using the same dynamical 
fields — at present the A-isobar is not considered as explicit degree of freedom 
in Ref. [215], but plays a numerically important role in the extraction of S 
from NN NNn. 

In section 4.3 it was shown that the differential cross section as well as the 
analysing power for the reaction pn ppT[~ is sensitive to an interference 
term of the s-wave pion production amplitude of the An amplitude (^Pq 
Sqs) and the p amplitudes of Aqi. ^Si SqP and ^Di S^p. Obviously, 
the 4-nucleon contact interaction contributes to the former. Thus, once a 
proper chiral perturbation theory calculation is available for the s-wave pion 
production — whose status will be discussed in the subsequent paragraphs — 
the reaction pn — > pp7r~ close to the production threshold might well be the 
most sensitive reaction to extract the parameter d. 

Let us now turn to the s-wave contributions. The leading diagrams containing 
nucleons only are shown in table 11; those that contain the A are shown in 
table 12. Again, the Weinberg scheme and the new scheme are compared. The 
list is complete up to NNLO in both schemes. Please note, however, that one 
class of diagrams {k and I) that is of NLO in the Weinberg counting in the 
new scheme is pushed to N^LO! 

Thus in the new counting scheme — contrary to the Weinberg scheme — the 
leading pieces of some loops appear one order lower than the tree level isoscalar 
rescattering amplitudes. As can also be read from the table, the corresponding 
order is m,r/-^Ar- If we consider, in addition, that due to the odd parity of the 
pion the initial state has to be a p-wave, the loops themselves have to scale 
as ■sfm^. Since no counter term non-analytic in rriT^ is allowed, the loops have 
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Table 11 

Comparison of the corresponding chiral order in the Weinberg scheme (p ~ m.,^) 
and the new counting scheme (p ~ ^/nhJd^ ) for several nucleonic contributions for 
s-wave pion production. Subleading vertices are marked as 0. Not shown explicitly 
are the recoil corrections for low-order diagrams. For example, recoil corrections to 
diagram b) appear at order pm-,^ / M'j^ . 

s-wscve diagrams 
(nucleons only) 



d) 




b) 



g) 



h) 




scale 


p ~ 


p ~ 




LO 


LO 


p/M 


NLO 


LO 




N^LO 


NLO 


pm^/M% 


N^LO 


N^LO 


ml/{pMN) 


NLO 


N^LO 


ml/pMl 


N^LO 


N^LO 
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Table 12 

Comparison of the corresponding chiral order in the Weinberg scheme (p ~ rriT^) 
and the new counting scheme (p ^ \/m.,^Miy ) for the leading and next-to-leading 
A contributions for s-wave pion production. Subleading vertices are marked as 0. 
Not shown explicitly are the recoil corrections for low-order diagrams. 

s-wave diagrams scale A ~ tut^ A ~ 

(Deltas only) \/rrh^^Mj^ 



pm^/AMAT NLO NLO 



p^/{AM^) N^LO NLO 



to be finite or to cancel exactly. This requirement is an important consistency 
check of the new counting scheme. 




The details of the loops calculations in threshold kinematics can be found in 
Ref. [204]. Here we only give the results. Note, we only evaluate the leading 
order pieces of the integrals corresponding to the diagrams of table (11). For 
example, in the integrals we drop terms of order compared to /q- 

After these simplifications, straightforward evaluation gives for the production 
amplitude from the loops with nucleons only 

A ^ 3 / - ^ W ^Nm^ , 



i«=i 

where n denotes the diagram (labels as in the figure) . For the isospin functions 
we find 
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Those can be easily evaluated in the different isospin channels. We find for 
i{TT3) — (TTsltrf + Le + if\ll) (it is sufficient here to look at the pp initial 
state only) 

^(11) = -1-^ + ^=0' ^(00) = -1 + + ^ = ■ (89) 

The first observation is that the NLO contributions are of the order of mag- 
nitude expected by the power counting, since 

where we used = 93 MeV. The power counting proposed in Refs. [163,203] 
thus is indeed capable of treating properly the large scale inherent to the the 
NN NNn reaction. 

We checked that our results for the individual diagrams agree with the leading 
non-vanishing pieces from the calculations of Rcf. [205] . In addition, for 
almost all diagrams given in table 11, Ref. [205] gives explicit numbers for 
the amplitudes in threshold kinematics. It is intriguing to compare those to 
what one expects from the different counting schemes. This is done in table 
13, where the first line specifies the particular diagram according to table 11 
and the second gives the result of the analytical calculation of Ref. [205] 
(normalized to the first column) . In the following two lines those numbers are 
compared to the expectations based on the counting schemes — first showing 
those for the Weinberg scheme and then those for the new counting scheme. 
As can clearly be seen, the latter does an impressive job of predicting properly 
the hierarchy of diagrams. Thus, at least when ISI and FSI are neglected, the 

In this reference the same choice for the pion field is made and thus a comparison 
of individual diagrams makes sense. Note that there is a sign error in the formula 
for diagram / in Ref. [205]. We are grateful to F. Myhrer for helping us to resolve 
this discrepancy. 

In Ref. [205] the full result for the particular loops are given. Thus, any loop 
contains higher order contributions. 
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Table 13 

Comparison of the results of the analytic calculation of Ref. [205] with the expec- 
tations based on the two counting schemes as discussed in the text. The diagrams 
are labeled as in Fig. 11 (the label shows that here the recoil term of diagram j 
is calculated; it appears at NNLO in the Weinberg scheme as well as in the new 
counting scheme). 



Diagram 


d 


e 


g 


3R 


k 


1 


Ref. [205] 


1.0 


-1.0 


0.1 


0.4 


0.03 


0.02 




1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


p ~ ^/tUt^Mn 


1.0 


1.0 


0.4 


0.4 


0.06 


0.06 



counting scheme proposed is capable of dealing with these large momentum 
transfer reactions. 

It is striking that the sum of loops in the case of the neutral pion production 
vanishes = in Eq. (89)). In addition, for neutral pion production 

there is no meson exchange current at leading order and the nucleonic current 
(diagrams b) and c) in table 11) gets suppressed by the poor overlap of the 
initial and final state wave functions (see discussion in sec. 1.4) — an effect 
not captured by the counting — and interferes destructively with the direct 
production off the Delta (diagrams a) and b) in table 12). Thus the first 
significant contributions to the neutral pion production appear at NNLO. 
This is the reason why many authors found many different mechanisms, all 
of similar importance and capable of removing the discrepancy between the 
Koltun and Reitan result and the data, simply because there is a large number 
of diagrams at NNLO. The situation is very different for the charged pions. 
Here there is a meson exchange current at leading order and there are non- 
vanishing loop contributions. We therefore expect charged pion production to 
be significantly better under control than neutral pion production and this is 
indeed what we found in the phenomenological model described in the previous 
section. 

Next let us have a look at the loops that contain Delta isobars (c.f. Fig. 
12). In Ref. [204] it was shown that the individual loops diverge already at 
leading order, because the Delta-Nucleon mass difference introduces a new 
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scale. Therefore, as was argued above, the sum of the diagrams has to cancel, 
for at NLO there is no counter term. Is is an important check of the counting 
scheme that this does indeed happen. We take this as a strong indication that 
expanding in x is consistent with the chiral expansion. 

The observation that we have now established a counting scheme for reactions 
of the type NN NNtt also has implications for the understanding of other 
reactions. One example is the analysis elastic 7rd scattering, which is commonly 
used to extract the isoscalar nN scattering length that is difficult to get at 
otherwise [167] . As mentioned above, due to chiral symmetry constraints the 
isoscalar scattering length does not get a contribution at leading order and 
thus, for an accurate extraction of this important quantity from data on ird 
scattering or on Trd bound states, a controlled calculation of the few body- 
corrections is compulsory. One of these corrections are the so-called dispersive 
corrections (see Ref. [221] and references therein): loop contributions to the 
elastic scattering that have intermediate two-nucleon states. Obviously, the 
imaginary part of those loops gives the essential contribution to the imaginary 
part of the nd scattering length , however there is also a real part to these 
loops that needs to be calculated within a scheme consistent with that used for 
the calculation of the other contributions. Within chiral perturbation theory 
that has not been done up to now. Given the progress reported here, such a 
calculation is now feasible. The same technique can then also be used to cal- 
culate the corresponding corrections for n^He scattering, recently calculated 
in chiral perturbation theory for the first time [223]. 

6.4 On the significance of off-shell effects 

A few years a ago there existed a strong program at several hadron facilities 
to measure bremsstrahlung in NN collisions with the goal to identify the 
true off-shell behavior of the A^A^ interaction. Indeed it was found that the 
predictions for several highly differential observables are significantly different 
when different NN potentials are employed. 

About one third of the imaginary part was found to be related to the reaction 
TTd^jNN [222]. 
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Fig. 29. Visualisation of Eq. (90). 



On the other hand, a change in the off-shell behavior of any T-matrix can be 
realized on the level of the effective interaction by a field redefinition and it 
is known since long that S matrix elements do not change under those trans- 
formations (as long as one works within a well defined field theory) [224,225]. 
Thus, for any given set of fields the off-shell amphtudes might well have a 
significant impact on the values of various observables, but it is an intrin- 
sic feature of quantum field theory that they can not be separated from the 
short range interactions constructed within the same model space. In a very 
pedagogical way those results were presented later in Refs. [226,227]. 

How can one understand this seemingly contradictory situation: on the one 
hand off-shell amplitudes enter the evaluation of matrix elements and in some 
cases influence significantly the result (c.f. discussion in previous sections), 
on the other hand they are claimed not to have any physical significance? 
Following Ref. [226] to start the discussion, let us consider some general half- 
off shell nN T-matrix. In a covariant form the list of its arguments contains 
the standard Mandelstam variables s, t as well as — the four momentum 
squared of the off-shell particle(s). For our discussion let this be the incoming 
pion and for simplicity omit all spin indices. We now define Tr(s, t, (f) via 

T{s, t, q^) = f (s, t, q^) + {q^ - ml)Tn{s, t, q^) , 

where T(s, t, q^) is arbitrary up to the condition that it has to agree with the 
on-shell amphtude for (f = m^. Obviously, this is always possible and we can 
assume Tr to be smooth around the on-shell point. Then, when introduced 
into the NN — > A^'A^tt transition amplitude, we get 

A = T ^ W^NN = f ^ W^NN + TrW^nn , (90) 
q^ — m^ q^ — 

where Wt^nn denotes the nNN vertex function. This example highlights two 
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important facts: i) a change in the off-shell dependence of a particular ampli- 
tude can always be compensated via an appropriate additional contact term. 
The only quantities that are physically accessible are S'-matrix elements and 
those are by definition on-shell; and ii) the contribution stemming from the 
off-shell 'kN T-matrix is a short range contribution. Therefore the distinction 
between short range and off-shell rescattering is artificial. 

6.5 Lessons and outlook 

The previous section especially made clear that it is rather difficult to con- 
struct a model that gives quantitatively satisfactory results for the reaction 
pp ppT[^ . On the other hand, it turned out that for the reaction pp pmr^ 
current models as well as the effective field theory approach do well. In the 
previous sections this difference was traced back to a suppression of meson 
exchange currents in the neutral pion production. What is the lesson to be 
learned from this? First of all, any model for meson production close to the 
threshold should contain the most prominent meson exchange currents; as long 
as these are not too strongly suppressed, one can expect them to dominate the 
total production cross section close to the threshold. However if there is no 
dominant meson exchange current, then there is a large number of sub-leading 
operators that compete with each other and make a quantitative understand- 
ing of the cross section difficult. The optimistic conclusion to be drawn from 
these observations is that (also for heavier mesons) the leading meson ex- 
change currents should give a reasonable description of the data, while the 
contributions from irreducible loops largely cancel. Obviously, in heavy meson 
production there is no reason anymore to consider pions only as the exchange 
particles. In this sense it is tt"^ production that is the typical case, whereas 
the 7r° is exceptional due to the particular constraints from chiral symmetry 
in that channel. Note that also in the case of pion photoproduction close to 
the threshold the 7r° plays a special role, in that the s-wave amplitude is 
dominated by loops [218]. 

Let us look in somewhat more detail at the production operator for neutral 
pion production. Within the model described above, the most prominent di- 
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Fig. 30. The one sigma exchange as it is perturbatively built up in the effective field 
theory, starting from the left with the lowest order diagram (NLO). The chiral order 
increases by one power in x between each diagram from the left to the right. 

agram for neutral pion production close to the production threshold is pion 
rescattering via the isoscalar pion-nucleon T-matrix that, for the kinematics 
given, is dominated by a one-sigma exchange (diagram b) of Fig. 2, where 
the T-matrix is replaced by the isoscalar potential given by diagram e) of 
Fig. 22). Within the effective field theory the isoscalar potential is built up 
perturbatively. This is illustrated in Fig. 30. As was shown above, the leading 
piece of the one-sigma exchange gets canceled by other loops that cannot be 
interpreted as a rescattering diagram (the sum of diagrams d), e) and f) of 
table 11 vanishes) and are therefore not included in the phenomenological ap- 
proach. This is an indication that in order to improve the phenomenological 
approach, at least in case of neutral pion production, pion loops should be 
considered as well . 

There is one more important conclusion to be drawn from the insights re- 
ported in the two previous sections: We can now identify what was missing 
in the straightforward extension of the Bonn potential reported in Ref. [24]. 
This is especially relevant, since we will see that the failure to describe the low 
energy pion production data in that approach does not point at a missing de- 
gree of freedom in the nucleon-nucleon phenomenology, but at an incomplete 
treatment of the cut structure. 

Diagrammatically, the dressed ttNN vertex function is given by 



Note that within the effective field theory approach the convergence of the series 
shown in Fig. 30 should also be checked, as pointed out in a different context in 
Ref. [228]. 
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Fig. 31. Leading cut structure of the one pion exchange potential in nucleon-nucleon 
scattering. 




where the solid dot indicates the so-called bare vertex (in this picture it is 
to be understood as a parameterization of both the extended structure of the 
nucleon as well as meson dynamics not included in the non-pole ttN T-matrix 
T^^, such as np [229] or na [230] correlations) as well as the dressing due to the 
ttN interactions parameterized in J"^^. The latter structure, however, contains 
a ttN cut not considered in Ref . [24] . It is straightforward to map the different 
three-body cuts in the one pion exchange potential to those that must occur 
in the pion production operator in order to allow for a consistent description of 
scattering and production. This is indicated in Fig. 31: the additional cut due 
to the pion dynamics in the vertex function is related to the pion rescattering 
diagram in the production operator. Therefore, at least close to threshold in 
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Ref . [24] , the most relevant mechanism that leads to inelasticities was missed. 

As we saw, however, as we move away from the near-threshold region the 
contribution from the Delta isobar becomes more and more significant. Con- 
sequently, the results of Ref. [24] look significantly better at higher energies. 

7 Production of heavier mesons 

We start this section with some general remarks and then discuss some special 
examples. The phenomenology will be discussed only very briefly; for details 
on this we refer to the original references or the recent reviews [1,2]. The list of 
reactions discussed in detail is by no means complete. For example, we will not 
discuss here 77' production, for this was already discussed in great detail in Ref. 
[2]. Neither will the production of vector mesons be discussed in detail. Here 
we refer to the recent contributions to the hterature (Ref. [59] and references 
therein) . 

As can be seen from the headlines already, special emphasis will now be put 
on the physics that can be extracted from the particular reactions. 

7.1 Generalities 

As the mass of the produced meson increases, the initial energy needs to in- 
crease as well and consequently also the typical momentum transfer at thresh- 
old. This has various implications: 

• the NN interaction needed for a proper evaluation of the initial state in- 
teraction is less under control theoretically; for the T = channel there is 
not even a partial wave analysis available above the 77-production threshold, 
due to a lack of data 

• a larger momentum transfer makes it more difficult to construct the pro- 
duction operator; at least one should expect that the exchanges of heavier 
mesons become even more relevant compared to the pion case 
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• it is not known, how to formulate a convergent effective field theory ^ 

• less is known about the interactions of the subsystems 

• the treatment of three-body singularities due to the exchange of light par- 
ticles requires more care. 

In the upcoming section we will not discuss in detail the role of the three- 
body singularities. Only recently was a method developed that in the future 
will allow inclusion of those singularities in calculations within the distorted 
wave born approximation [231]. This method was already applied in a toy 
model calculation [232] , where its usefulness was demonstrated. 

In calculations these singularities can occur only if the initial state interaction 
is included. For the production of mesons heavier than the 77, however, no 
reliable model is available at present. Up to now no model for the production 
of heavy mesons takes three-body singularities into account, and thus their 
role is yet unclear. 

1.2 Remarks on the production operator for heavy meson production 

As in sec. 6 we will only look at analyses, that work within the distorted wave 
born approximation. For the two-baryon states nothing changes, other than 
that the nucleon-nucleon interaction needed for the initial state interaction 
gets less reliable from the phenomenological point of view as we go up in en- 
ergy. The construction of the production operator, however, is now even more 
demanding, for with increasing momentum transfer heavier exchange mesons 
can play a significant role. Therefore, in order to get anything useful out of 
a model calculation, as many reaction channels as possible should be studied 
simultaneously. Only in this way can the phenomenological model parameters 
can be fixed and useful information be extracted. In this context also the si- 
multaneous analysis of pp and pn induced reactions plays an important role. 

In principle one might expect an approach like that presented in sec. 6.3 to work 
for the production of all Goldstone bosons. However, the next threshold after pion 
production is that for eta production and it already corresponds to an initial mo- 
mentum of pi=770 MeV. Thus the expansion parameter would be ^Jm^jM^ = 0.8. 



105 



for different isospin structures in the production operator will lead to very dif- 
ferent relative strengths of the two channels (c.f. discussion in sec. 1.5). Then, 
once a basic model is constructed that is consistent with most of the data, 
deviations in particular reaction channels can be studied. 

In Ref. [34] the model described in sec. 6.2 was extended to rj production 
in nucleon-nucleon collisions, where the single channel ttN T-matrix used in 
the pion production calculations was replaced by the multi channel meson- 
baryon model of Ref. [233] to properly account for the exchange of heavy 
mesons. So far only s-waves were considered in the calculation. Results for 
the various total cross sections using two different models for the final state 
A'^A^ interaction are shown in Fig. 32. The studies of Ref. [34] indicate, that 
a complete calculation for A^A^ NNrj will put additional constraints on 
the relative phases of the various meson-baryon— > riN transition amplitudes. 
Unfortunately, up to now the rj is the heaviest meson that can be investigated 
using this kind of microscopic approach, for there is no reliable model for the 
A'A^ interaction for energies significantly above the rj production threshold. 

The strategy that therefore needs to be followed is to study many reaction 
channels consistently. This was done, for example, in a series of papers by 
Kaiser and others for tt, rj, rj' [155], uj [234] as well as strangeness produc- 
tion [235] and by Nakayama and coworkers for cu [59,236,237], [237,238], 
rj [239,240] and rj' [241] production. Both groups construct the production 
operator as relativistic meson exchange currents, where the parameters are 
constrained by other data such as decay ratios. The striking difference be- 
tween the two approaches is that the former studies s-waves only, does not 
consider effects of the initial state interaction and treats the final state interac- 
tion in an approximate fashion. The latter group includes the ISI through the 
procedure of Ref. [40] (c.f. sec. 3), treats the FSI microscopically and includes 
higher partial waves as well. It is important to stress that, where compatible, 
the two approaches give quahtatively similar results, as stressed in Ref. [234]. 

Recently it was observed that the onset of higher partial waves can strongly 
constrain the production operator [240]. In the same reference it was demon- 
strated that to unambiguously disentangle effects from FSI or higher partial 
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waves, polarization observables are necessary. Thus one should expect that 
once a large amount of polarized data is available on the production of heavy 
mesons, the production operators and thus the relevant short-range physics 
for the production of a particular meson, can be largely fixed. As an example, 
this particular issue is discussed in detail in the next section. 

1.3 The reaction NN NNrj or properties of the 5'ii(1535) 

The T) meson is a close relative of the pion, since it is also a member of the 
nonet of the hghtest pseudo scalar mesons. It is an isoscalar with a mass of 
547.3 MeV — the difference in mass from the pion can be understood from its 

content of hidden strangeness through the Gell-Mann-Okubo mass formula 
[242,243]. 

As for the pion, the rj couples strongly to a resonance, but a resonance with 
different quantum numbers: is it the positive parity A(1232) for the pion 
that in the A rest frame leads to p wave pion production, the rj production 
is dominated by the negative parity S'ii(1535) , which leads to s wave rj 
production in its rest frame. As a consequence near-threshold rj production 
is completely dominated by the 1511(1535) resonance. This statement can well 
be reversed: studying rj production in various reaction channels close to the 
production threshold allows selective study of the properties of the 5*11(1535) 
resonance in various environments. In this context it is interesting to note that 
the nature of this particular resonance is under discussion: within the so-called 
chiral unitary approach the resonance turns out to be dynamically generated 
[244,245], while on the other hand, detailed studies within the meson exchange 
approach [246], as described in sec. 6, as well as quark models (see Ref. [247] 
and references therein) call for a genuine quark resonance. One expects that 
a molecule behaves differently in the presence of other baryons than a three- 
quark state due to the naturally enhanced affinity to meson baryon states of 
the former. In this sense data on rj production in few-baryon systems should 

The quantum numbers are chosen in accordance with the partial wave in which 
the resonance would appear in ttN scattering. Therefore positive parity resonances 
appear in odd partial waves. 
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Fig. 32. Results for the total cross sections of the reactions pp — > ppr], pn pur], 
and pn drj from Ref. [34] employing different NN models for the final state 
interaction. The solid lines represent the results with the CCF NN model [175] 
whereas the dashed-dotted lines were obtained for the Bonn B model [25]. Data are 
from Refs. [33,109-112,120]. 

provide valuable information about this lowest negative parity excitation of 
the nucleon. 

The reaction NN NNrj was studied intensively both theoretically 
[34,155,239,240,248-259] as well as experimentally (c.f. table 8) in recent years. 
The outstanding feature of rj production in NN collisions is the strong effect of 
the T]N FSI, most visible in the reaction pn drj, where in a range of 10 MeV 
the amplitude grows by about an order of magnitude [120]. This phenomenon 
was analyzed by means of Faddeev calculations, showing this enhancement to 
be consistent with a real part of the rjN scattering length of 0.4 fm [248]. It is 
interesting that the same value of the rjN scattering length was recently ex- 
tracted from an analyses of the reaction 'yd rjpn [260] and that this value is 
consistent with that stemming from the microscopic model for meson-baryon 
scattering of Ref. [233]. First three body calculations for the reaction with 
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a two-nucleon pair in the continuum show only a minor impact of the rjN 
interaction on the invariant mass spectra [261]. 

Most model calculations for the production operator share the property that 
the production operator was calculated within the meson exchange picture . 
This was then convoluted with the NN FSI, treated either microscopically or 
in some approximate fashion. In most analyses it was found that the domi- 
nant production mechanism is via the 5'ii(1535), with the exception of Ref. 
[249] , where the reaction pp — > ppr] turned out to be dominated by isoscalar 
meson exchanges, in analogy with the heavy meson exchanges in the reaction 
pp — > pp7r°. The largest qualitative differences between the various models 
is the relative importance of the exchanged mesons. In Refs. [155,252-255] 
the p-exchange turned out to be the dominant process, whereas it played a 
minor role in Refs. [34,256,249]. In Ref. [259] the p-exchange was not even 
considered, but 77 exchange turned out to be significant there. It remains to 
be seen, however, if models that are dominated by isoscalar exchanges, Refs. 
[259,249], are capable of accounting for the large ratio of pn- to pp-induced 
eta production, as discussed in sec. 1.5. 

As stressed above, in this section we do not want to focus on the details of 
the dynamics of the production operator, but instead on the physical aspect 
of what we can learn about the y^ii resonance from studying eta meson pro- 
duction in A'^A'^ collisions. Thus, in the remaining part of this section we will 
restrict ourselves to a model-independent analysis of a particular set of data 
based on the amplitude method introduced in sec. 4.3. Here we closely follow 
the reasoning of Ref. [240] . 

To be definite we will now concentrate on the measurement of angular distri- 
butions and invariant mass spectra for pp ppr] at Q=15 MeV [113]. The 
experiment shows that the angular distribution of both the two-proton pair 
and of the rj are fiat, suggesting that only s-waves are present. On the other 
hand, the invariant mass distribution da/drupp deviates significantly from what 

The only exception to the list given above is Ref. [257] , where the ratio of pp- to 
pn-induced 77 production is explained by the instanton force. The relation between 
this result and the hadronic approach is unclear [257] . 
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Fig. 33. Invariant mass spectrum for the two-proton system for the reaction 
pp — > pprj at Q=15 MeV. The dashed line shows the Ss contribution, the dot-dashed 
line the Ps distribution and the solid line the incoherent sum of both. The data is 
taken from Ref. [113]. 

is predicted based on the presence of the strong pp FSI only (c.f. discussion 
of sec. 2). This is illustrated in Fig. 33, where the distribution for pp S-wave, 
7] s-wave {Ss) is shown as the dashed line. Since the rjN interaction is known 
to be strong, its presence appears as the natural explanation for the devia- 
tion of the dashed line from the data. The structure, however, is even more 
pronounced at Q = 42 MeV (c.f. Ref. [240]) in contrast to what should be 
expected for the rjN interaction. 

On the other hand, the discrepancy between the data and the dashed curve 
can be well accounted for by a Ps configuration (given by phase-space times 
a factor of p'^, c.f. sec. 2.1). This is also illustrated in Fig. 33, where the pure 
Ps configuration is shown as the dashed-dotted line and the total result — the 
incoherent sum of Ps and Ss cross section — as the solid line. 

How is this compatible with the angular distribution of the two nucleon system 
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being flat? This question can be very easily addressed in the amplitude method 
described in sec. 4.3. As long as we resrict ourselves to rj s-waves and at most 
P waves in the NN system, the only non-vanishing terms in Eq. (48) are Q, 

— * 

containing the amplitude for the Ss final state, and A, containing the two 
possible amplitudes for the Ps final state. The explicit expressions for the 
amplitudes were given previously in Eqs. (55). 

In this particular case we find for example (c.f. Eqs (49)- (54)) 

^o=i(igr + |AT) 

Aoiao = 

A,,ao^\{\Q\'-\A\') 
where 

\Q\' = \ai\\ 

\A\'^p''(\a, - {l/3)as\' + x^asl' - {2/3)Re{ala,)}) , 

and p'x = {p' ■ p). Here the amplitudes Oi, 02 and 03 correspond to the tran- 
sitions ^Pq —^^ Sqs, ^Sq -^^ PqS, and — P2S, respectively. From these 
equations we directly read that 

• if 03 is negligible, the pp differential cross section da/dx is flat 

• the observable ao{l — A^x) directly measures the A^A^ P-wave admixture 

• any non-zero value for the analysing power is an indication for higher partial 
waves for the 77. 

It is easy to see that for 03 = the differential cross section has to be fiat, since 
then the only partial wave that contributes to the A^A^ P-waves is ^Sq Pqs 
and a J = initial state does not contain any information about the beam 
direction. 

Certainly, the 77A" final state interaction has to be present in this reaction 
channel as well and it is important to understand its role in combination with 
the two-nucleon continuum state. It should be clear, however, that in order 
to understand quantitatively the role of the 77A" interaction, it is important to 
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pin down the NN P-wave contribution first. This is why a measurement of 
A^x for pp — > ppr) is so important. 

What does it imply if the conjecture of a presence of NN P-waves at rather 
low excess energies were true? In Ref. [240] it is demonstrated that the need 
to populate predominantly the ^-S'o — >^ Pqs instead of the — >^ P2S very 
strongly constrains the NN — > -S'nA^ transition potential. This indicates that 
a detailed study of NN — > NNrj should reveal information about the Sn in 
a baryonic environment that might prove valuable in addressing the question 
of the nature of the lowest negative-parity nucleon resonance. 

It should be mentioned that in Ref. [262] an alternative explanation for the 
shape of the invariant mass spectrum was given, namely an energy depen- 
dence was introduced to the production operator. Based on the arguments 
given in the previous chapters, we do not believe that this is a natural ex- 
planation. However, as we just outlined, with A^x an observable exists that 
allows to unambiguously distinguish between the two possible explanations. 
The experiment is possible at COSY [263]. 

7.4 Associated strangeness production or the hyperon-nucleon interaction 
from production reactions 

In sec. 5 the small breaking of SU (2) isospin symmetry present in the strong 
interaction was discussed. If we include strange particles in the analysis we can 
also study the breaking of fiavor SU (3). It is well-known that the light mesons 
and baryons can be arranged according to the irreducible representations of the 
group SU{3). The mass splittings within the multiplet can be well accounted 
for by the number of strange quarks in some baryon or meson. However, not 
much is known about the dynamics of systems that contain strangeness. Many 
phenomenological models for, e.g., hyperon-nucleon scattering [268-271] use 
the fiavor SU{3) to fix the meson baryon-meson couplings. The remaining 
parameters, like the cut-off parameters, are then fit to the data. As we will 
discuss below, so far the existing data base for hyperon-nucleon scattering is 
insufficient to judge, if this procedure is appropriate. 
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Fig. 34. Comparison of the quality of available data for the reactions pA elastic 
scattering (data are from [264-266]) and pp K+Ap at TLab=2.3 GeV [267]'^'^ . 
In both panels the red curve corresponds to a best fit to the data. In the left panel 
the dashed lines represent the spread in the energy behavior allowed by the data, 
according to the analysis of Ref. [266]; analogous curves in the right panel would lie 
almost on top of the solid line and are thus not shown explicitly. 

As was stressed previously, effective field theories provide the bridge be- 
tween the hadronic world and QCD. In connection with systems that contain 
strangeness there are still many open questions. Up to now it is not clear if the 
kaon is more appropriately treated as heavy or as light particle. In addition, 
in order to establish the counting rules it is important to know the value of 
the SU{3) chiral condensate. For a review this very active field of research as 
well as the relevant references we refer to Ref. [272]. 

To further improve our understanding of the dynamics of systems that contain 
strangeness, better data is needed. The insights to be gained are relevant not 
only for few-body physics, but also for the formation of hypernuclei [273], 
and might even be relevant to the structure of neutron stars (for a recent 
discussion on the role of hyperons in the evolution of neutron stars see Ref. 
[274]). Naturally, the hyperon-nucleon scattering lengths are the quantities of 
interest in this context. 
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Fig. 35. Values allowed for the spin singlet and spin triplet scattering length by the 
AN elastic scattering data according to Ref. [266]. The dark shaded area denotes 
the la range for the parameters and the light shaded area the 2a range. The cross 
shows the best fit value (og = 1.8 fm and at = 1.6 fm). 

In the left panel of Fig. 34 we show the world data set for elastic AA^ scatter- 
ing. In Ref. [266] a Likelihood analysis based on the elastic scattering data was 
performed in order to extract the low energy AA^ scattering parameters. The 
resulting contour levels are shown in Fig. 35, clearly demonstrating that the 
available elastic hyperon-nucleon scattering data do not significantly constrain 
the scattering lengths: the data allows for values of (-1, 2.3) as well as (6,1) 
(all in fermi) for (as,aj) respectively. Later models were used to extrapolate 
the data. However, also in this way the scattering lengths could not be pinned 
down accurately. For example, in Ref. [269] one can find six different models 
that equally well describe the available data but whose (S-wave) scattering 
lengths range from 0.7 to 2.6 fm in the singlet channel and from 1.7 to 2.15 fm 
in the triplet channel. Production reactions are therefore a promising alterna- 
tive. In the literature the reactions K^d — > '-jh.n [9], — > K^An (Ref. [53] 
and references therein) and pp pK^A [122] were suggested. Therefore the 
method described in sec. 2.2, that applies to the latter two reactions, is an im- 
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portant step towards a model-independent extraction of the hyperon-nucleon 
scattering lengths. 

A natural question that arises is the quality of data needed e.g. for the reac- 
tion pp — > pK^h. in order to significantly improve our knowledge about the 
hyperon-nucleon scattering lengths. In Ref. [32] it is demonstrated that data 
of the quality of the Saclay experiment for pp — > K^X [267], shown in right 
panel of Fig. 34 that had a mass resolution of 2 MeV, allows for an extrac- 
tion of a scattering length with an experimental uncertainty of only 0.2 fm. 
Note, however, that the actual value of the scattering length extracted with 
Eq. (22) from those data is not meaningful, since the analysis presented can 
be applied only if just a single partial wave contributes to the invariant mass 
spectrum. The data set shown, however, represents the incoherent sum of the 
^Si and the ^-S'o hyperon-nucleon final state. 

The two spin states can be separated using polarization measurements. In 
sec. 4.3.1 as well as in Ref. [32] it was shown, that for the spin singlet final 
state the angular distributions of various polarization observables are largely 
constrained. This is sufficient to extract the spin dependence of the A-nucleon 
interaction from the reaction pp — > pK'^A. One more issue is important to 
stress: to make sure that the structure seen in the invariant mass spectrum 
truly stems from the final state interaction of interest, a Dalitz plot analysis 
is necessary, for resonances can well distort the spectra. This is discussed in 
detail in sec. 4.1 as well as at the end of sec. 2.2. 

7. 5 Production of scalar mesons or properties of the lightest scalar 

The lightest scalar resonances ao(980) and /o(980) are two well-established 
states seen in various reactions [146], but their internal structure is still under 
discussion. Analyses can be found in the literature identifying these structures 
with conventional qq states (see Ref. [275] and references therein), compact 
qq-qq states [276,277] or loosely bound KK molecules [278,279]. In Ref. [280] 

Shown is only the low tail of the data taken at TLab = 2.3 GeV and an angle 
of 10" 
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A = — •' — 



Fig. 36. Graphical illustration of the leading contribution to the fo—ao mixing matrix 
element A defined in Eq. (91). 

a close connection between a possible molecule character of the light scalar 
mesons and chiral symmetry was stressed. It has even been suggested that at 
masses below 1.0 GeV a complete nonet of 4-quark states might exist [281]. 

Resolution of the nature of the hght scalar resonances is one of the most 
pressing questions of current hadron physics. First of all we need to understand 
the multiplet structure of the light scalars in order to identify possible glueball 
candidates. In addition, it was pointed out recently that ao(980), /o(980) as 
well as the newly discovered Dg [282] might well be close relatives [283]. Thus, 
resolving the nature of the light scalar mesons allows one simultaneously to 
draw conclusions about the charmed sector and will also shed light on both 
the confining mechanism in light-light as well as light-heavy systems. 

Although predicted long ago to be large [284] , the phenomenon of Oq — /o mix- 
ing has not yet been estabhshed experimentally. In Ref. [284] it was demon- 
strated that the leading piece of the /o — ao mixing amplitude can be written 



as^^ 



A = (/o|r|ao) = igf,KK9aoKKV~s {pro - Pk+) + O , (91) 



where pk denotes the modulus of the relative momentum of the kaon pair 
and the effective coupling constants are defined through F^.^^ = 91kkPk- 
Obviously, this leading contribution is just that of the unitarity cut of the 
diagrams shown in Fig. 36 and is therefore model-independent. Note that in 
Eq. (91) electromagnetic effects were neglected, because they are expected to 
be small [284]. 



Here we deviate from the original notation of Achasov et al. in order to introduce 
dimensionless coupling constants in line with the standard Flatte parameterization. 
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Fig. 37. Modulus of the leading piece of the mixing amplitude A defined in Eq. (91). 
The two kinks occur at the K+K' (at 987.35 MeV) and the K°K° (995.34 MeV) 
threshold respectively. 

The contribution shown in Eq. (91) is unusually enhanced between the K'^K~ 
and the K^K^ thresholds, a regime of only 8 MeV width. Here it scales as 
(this formula is for illustration only — the Coulomb interaction contributes with 
similar strength to the kaon mass difference [285]) 



where m^, ma and rus denote the current quark mass of the up, down and 
strange quark, respectively, and m — (m^ + md)/2. This is in contrast to 
common CSB effects ^° which scale as (m^ — mij)/{m + m^), since they have 
to be analytic in the quark masses. It is easy to see that away from the kaon 
thresholds A returns to a value of natural size. This ^/s dependence of A is 
depicted in Fig. 37. 

Within a microscopic model for tttt and KK scattering the mixing of ao and /o 
was studied in Ref . [286] . Within this model both resonances are of dynamical 
origin. The only mixing mechanism considered was the meson mass differences. 
Within this model the predictions of Ref. [284] were confirmed. 



Here we denote as common CSB effects those that occur at the Lagrangian level. 
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As was demonstrated by Weinberg for the case of the deuteron [287], the 
effective couphngs of resonances to the continuum states contain valuable in- 
formation about the nature of the particles. In the case of the deuteron this 
analysis demonstrated that the effective coupling of the deuteron to the pn 
continuum, as can be derived from the scattering length and the effective 
range, shows that the deuteron is a purely composite system. Recently it was 
demonstrated that, under certain conditions that apply in the case of Qq and 
/o, the analysis can be extended to unstable scalar states as well [288]. Ac- 
curate data on the effective couplings of the scalars to kaons should therefore 
provide valuable information about their nature. As was argued in sec. 5, the 
occurrence of scalar mixing dominates the CSB effects in production reactions. 
Therefore quantifying the mixing matrix element might be one of the cleanest 
ways to measure the effective decay constants of and /q. In sec. 5 it was 
demonstrated that from studies of Tir] production in NN and dd collisions one 
should be able to extract the Cq — /o mixing amplitude. Those arguments were 
supported in sec. 4.3.2, where it was shown, that the reaction pp — > dK^K^ is 
indeed dominated by KK s-waves. We should therefore expect a significant 
signal for the mixing as well. In the years to come we can thus expect the 
experimental information about the scalar mesons to be drastically improved. 

8 Meson production on light nuclei 

In this section we wish to make a few comments concerning meson production 
on light nuclei. Note, dd induced meson production was discussed to some 
extent in section 5. A more extensive discussion can be found in Ref. [2]. 

8.1 Generalities 

Almost all general statements made for meson production in NN coUisions 
apply equally well for meson production involving light nuclei. Naturally, now 
the selection rules are different (see discussion at the end of sec. 1.6) and the 
possible breakup of the nuclei introduces additional thresholds that must be 
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Fig. 38. Possible diagrams to contribute to the reaction Nd — ^ 3Nx. In diagram a) 
and b) the three nucleons in the final state from a bound state, whereas diagram 
c) shows the dominant diagram in quasi-free kinematics. Here it is assumed that 
only a particular nucleon pair (to be identified through proper choice of kinematics) 
interacts in the final state. 

considered in theoretical analyses. 

A comparison of meson production in two-nucleon collisions and in few- 
nucleon systems should improve our understanding of few-nucleon dynamics. 
As a result of a systematic study of, for example, both A^A^ and pd induced 
reactions, a deeper insight into the importance of three-body forces should be 
gained. For example, a recent microscopic calculation using purely two-body 
input to calculate the reaction Tr+^i/e — > ppp, found that the data [289,290] 
call for three-body correlations [291], confirming earher studies [292]. It should 
be stressed, however, that this field is still in its infancy and a large amount 
theoretical effort is urgently called for. 

8.2 2 — >• 2 reactions 

The reaction channel studied best up to date is pd — >^ Hex. Data exist mainly 
from SATURNE for x = tt [293], x ^ ou [294], x ^ rj [295,296], as well as x = 77' 
and (f) [297]. Most of the experiments were done using a polarized deuteron 
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beam. In addition the reactions dd — > an° [140] and dd — > ar] [298] were 
measured. 

The reaction pd He7r° could be understood quantitatively from diagram 
a) of Fig. 38 [299]. However, it turned out that diagram a) alone largely under 
predicted the data [300,301] for the reaction pd -^^ Rer}, whereas diagram b) 
contributes sufficiently strongly to allow description of the data [300,302,303]. 
In all these approaches the individual amplitudes were taken from data di- 
rectly; so far no microscopic calculation exists. The relevance of the pion ex- 
change mechanism was explained in Ref. [304] as what the authors called a 
kinematic miracle: in the near-threshold regime the kinematics for eta pro- 
duction almost exactly matches that for pp — > dir'^ followed by 7r+n r]p, 
with the intermediate pion on-shell. 

The striking feature of reactions with an 77-nucleus final state is the pro- 
nounced energy dependence that was already described in sec. 7.3 for the re- 
action pn — > dr) — often interpreted as a signal of an existing (quasi)bound state 
in the 77-nucleus system. Indeed, since the quark structure of the rj uu + dd 
is proportional to the number operator, one should expect the 77-nucleus in- 
teraction to get stronger with increasing number of light quark flavors present 
in the interaction region. For a more detailed discussion of this issue we refer 
to Ref. [2]. 

8.3 Quasi-free Production 

A dcutcron is a loosely bound state of a proton and a neutron. For properly 
chosen kinematics, the deuteron can therefore be viewed as an effective neu- 
tron beam/target as alternative to neutron beams (see Ref. [2] and references 
therein). The corresponding diagram is shown in Fig. 38c). 

Obviously, the existence of a three-nucleon bound states already shows that 
there must be a kinematic regime, where the interaction of all three nucle- 
ons in the final state is significant, namely, when all three have small relative 

Note that this reaction is charge symmetry breaking and was discussed in sec. 5. 
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Fig. 39. Momentum distribution of proton-spectators in the reaction pd^STr'^2Tr~ p. 
The data are taken from Ref. [305] while the solid line is is based on a calculation 
that includes the quasi-free production as well as meson rescattering. The figure 
is taken from [306], where also details about the corresponding calculation can be 
found. 

momenta. On the other hand, if the spectator nucleon escapes completely un- 
affected, its momentum distribution should be given by just half the deuteron 
momentum, convoluted with the deuteron wavef unction times the phase-space 
factors. Thus, we should expect a momentum distribution for the spectator 
that shows in addition to a pronounced peak from the quasi-free production 
a long tail stemming from rescattering in the final state. 

Experimental data exist for the proton-spectator momentum distribution for 
the reaction pd— >37r+27r~p [305]. In Fig. 39 this data is shown together with 
the results of a calculation [306] that considers both the quasi-free piece as 
well as a rescattering piece. The data clearly shows the quasi-elastic peak. 

Thus, when considering spectator momenta that are of about 100 MeV or 
less, to a good approximation, the reaction should be quasi-free. 

Experimentally this conjecture can be checked by comparing data for pp in- 
duced reactions with those stemming from a pd initial state with a neutron 
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spectator. Those comparisons were carried out at TRIUMF for pion produc- 
tion [74,75] as well as at CELSIUS for rj production [33,120], showing that the 
quasi-elastic assumption is a valid approximation. 

In sec. 5 it was argued that a forward-backward asymmetry in the reaction 
pn — > d{7r^'r])s-wave is a good system from which to extract the /o — oq mixing 
matrix element. If investigated at COSY, this reaction can only be studied 
with a deuteron as effective neutron target. Since we are after a CSB effect, 
the expected asymmetry is of the order of a few percent. Thus we need to en- 
sure a priori that the spectator does not introduce an asymmetry of this order 
through its strong interaction with, say, the deuteron. Theoretically it would 
be very demanding to control such a small effect in a four-particle system. 
Fortunately, we can test experimentally to what extent the spectator intro- 
duces a forward-backward asymmetry by studying, for example, the reaction 
pd — > TT'^dn. If the reaction were purely quasi-elastic, the angular distribution 
of the TT+ii system in its rest frame would be forward-backward symmetric, for 
it would be stemming from a pp initial state. However, any interaction of the 
spectator with either the pion or the deuteron, should immediately introduce 
some forward-backward asymmetry. Experimental investigations of this very 
sensitive test of the spectator approach are currently under way [307]. 

9 Summciry 

The physics of meson production in micleon-nucleon collisions is very rich. 
The various observables that are nowadays accessible experimentally are, for 
example, influenced by baryon resonances and final state interactions as well 
as their interference. It is therefore inevitable to use polarized observables to 
disentangle the many different physics aspects. 

From the authors personal point of view, the most important issues for the 
field are: 

• that NN and dd induced reactions are very well suited for studies of 
charge-symmetry breaking [5]. Especially, investigation of the reactions 
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pn — > d{Tr'^r))s-wave and dd — > a{7r'^r))s-wave should shed hght on the na- 
ture of the scalar mesons. For the experiments planned in this context we 
refer to Ref. [308]; 

• that sufficiently strong final state interactions can be extracted from produc- 
tion reactions with large momentum transfer, as, for example, pp — > pK'^A. 
The condition for an accurate extraction is a measurement with high reso- 
lution, as should be possible with the HIRES experiment at COSY [309] ; 

• that an effective field theory was developed for large momentum transfer 
reactions such as A^A'^ — > NNn. Once pushed to sufficiently high orders, 
those studies will not only provide us with a better understanding of the 
phenomenology of meson production in A'^A'^ collisions but also allow us 
to identify the charge symmetry breaking operators that lead to the cross 
sections reported in Refs. [6,140]; 

• that resonances can be studied in A^A^ induced reactions. Those studies are 
complementary to photon induced reactions, since the resonances can be 
excited by additional mechanisms, which can be selected using, for example, 
spin particles as spin filter [4]. In addition, also the properties of resonances 
in the presence of another baryon can be systematically investigated. 

There are very exciting times to come in the near future, for not only does the 
improvement of experimental apparatus in recent years permit the acquisition 
of a great deal of high-precision data but also because or theoretical under- 
standing has now also improved to a point that important physics questions 
can be addressed systematically. 
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Appendices 

A Kinematical variables 

In this report we only deal with reactions, in which the energies of the final 
states are so low that a non-relativistic treatment of the baryons is justified. 
This greatly simplifies the kinematics. Thus, in the center of mass system we 
may write for the reaction NN — > B1B2X, 



where p denotes the momentum of the initial nucleons, p' the relative momen- 
tum of the final nucleons and q' the center of mass momentum of the meson 
of mass rrix- The reduced mass of the outgoing two baryon system is denoted 
by/xi2 = MiM2/(Mi + M2). 

The kinematical variable traditionally used for the total energy of a meson 
production reaction is 77, the maximum meson momentum in units of the 
meson mass. Then one gets from Eq. (A.l): 



Another often used variable for the energy of a meson production reaction is 




(A.l) 




(A.2) 



g = - (Ml + M2 + m^) . 



(A.3) 



It is straightforward to express 77 in terms of Q: 




(A.4) 
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Table A.l 

The values for T^ab, Q and r] for the different channels for pion production. 

where Mf — Mi + M2 and the reduced mass of the full system is given by 
II — MfTJix/ {Mf + mj:). For the latter approximation we used that in the close 
to threshold regime Q <^ {Mf + rrix). 

Traditionally 77 is used for the pion production only, whereas Q is used for 
the production of all heavier mesons. To simplify the comparison of results for 
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those different reaction channels, we present in table A.l the various values 
for r), Q as well as T^ab for the different pion production channels. 



For completeness we also give the relation to the Laboratory variables 



Lab 



2M 



Plab-Tlb + '^MTLab. 



(A.5) 
(A.6) 



B Collection of useful formulas 



B. 1 Definition of Coordinate system 



In this appendix wc give the explicit expressions for the vectors relevant for the 
description of a general 2—^3 reaction. We work in the coordinate system, 
where the z-axis is given by the beam momentum p. These formulas are 
particularly useful for sec. 4.3. Thus we have 



p = p 







p' = p' 



\ cos(^p/) / 
Prom these one easily derives 



q' = q' 



^ Sm{6q') COs{(f)q') ^ 

sm{9q/) sin(0g/) 

\ COs{9q') J 



i{p X p') — ipp' 



i{p' X q')^ip'q' 



^ — sin(^p/) sin(0p/) ^ 
sin(^p/) cos(0p/) 




\ 



cos[ti qi ) sm[ti pi ) sm[(pp' ) — cos(^t/p/ ) sm[Hq' ) sm^^^g- 
- cos(^^q/) sin(^^p/) cos((;6p/) + cos(^^p/) sin(^^y/) cos(0q') 
— sin(^p/) sm{9qi) sm{(f)q' — 0^/) 



(B.l) 



as well as the analogous expression for i{px q'). 
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B.2 Spin traces 



In this appendix some relations are given, that are useful to evaluate expression 
that arise from in the amplitude method appearing described in sec. 4.3. 



The sum over the spins of the external particles leads to traces in spin space, 
such as 



tr(aO = 0, 
tr((7i(7j) = 25ij , 
tr((7j(7;3(7j) = 2ieii3j , 



(B.3) 
(B.4) 
(B.5) 
(B.6) 



To identify dependent structures the following reduction formula is useful 



a' 6- (cx d)^ 

{a -b) a ■ {cx d) + {a ■ d) a ■ {b X c) + {a ■ c)a ■ (dxb) . (B.7) 

To recall the sign or, better, the order of the vectors appearing, observe that 
on the right hand side the vectors other than a are rotated in cyclic order. To 
prove Eq. (B.7) we use 



^jma {^ka^li ^ki^la) ^il^jmk ^ki^jml ■ 



C Partial wave expansion 



In this appendix we give the explicit relations between the partial wave ampli- 
tudes for reactions of the type NN — > NNx and the spherical tensors defined 
in Eq. (34), where x is a scalar particle. The relations between the spherical 
tensors and the various observables is given in tables 4 and 5. 



127 



In terms of the partial wave amplitudes, we can write for two spin-^ particles 
in the initial state 



rpk:iq:i,k4,q4, 



1 



J2{S'M's,p',q'\M\SMs,p} {SMs,m^\S'M's,p\q') 



(/l)U/2)t|0/, 



(2L + 1)(2L + 1) 
\ (2J + 1)(2J+1) 

X {S'M's, L'M[|/m;)(/m;, l'm[\JMj){SMs, LO\JMj) 
X {S'M's, L'M'S'M'SQ'M'., l'fn'i\JMj){SMs, LO\JMj) 



X {SMs\Ti%jS^JSMs){S'M's\r;^^^^^^ 
X Y,^'^ {q ')Yl'm', {P ')Yv^'^ {q ')*1i^m^ {p T 
X M"(s,e)M"(s,e)^ . 



(/l)t^(/2)t| 



(C.l) 



In order to proceed the following identities are useful [310]: 



^/imi(p)^Z2m2(p) = II 



Im 



\ 



(2/1 + 1) (2/2 + 1) 

{21 + l)47r 
X {hmi, I2m2\lm){h0, l20\lO)Yim{p) 



,1 



-yV2k + l{-a,k 



(C.2) 
(C.3) 



The latter, for instance, allows evaluation of the matrix element of the spin 
operators: 
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{SMs I r£,r£, | SMs) 

= E(^^5|r£jmim2)(mim2|rg)j^M5) 

x(l(M5 - m2)|r(f,Jmi)(m2|r(J^J^(M5 - mi)) 
= ^(-)'^i+'/y(2fci + l)(2fc2 + l) 
x(>5M5|^(M5 - ma), ^m2)(-SM5|^mi, ^{Ms - mi)) 

x{^{Ms-m2),ki (-gi)|imi)(^m2,A;2 (-g2)|^(M5 - mi)) .(C.4) 

It is convenient to couple the remaining spherical harmonics to a common 
angular momentum and to define 

^>ZMjp)nm,(^) ECLMl, lrhi\XQ)B9.^^{q,p) (C.5) 

A 

where we used the fact that the sum of the projections turns out to be equal 
to Qi + q2 — Q; B is then 

Bfi,i^,P) = E ^{Lf^LMmYi^^mjQ) ■ (C.6) 
and normalized such that 

j dnpdngB9.^^{q,p) = SxoSioS'Jq^^ . (C.7) 

Some properties of B are derived in the next section. 

After putting together the individual pieces we arrive at the final result: 

LIX 

where p = {hqi, A;2?2, hqs, hq^}, and 

A'fr , = E C^ffM'^iM'') t (C.8) 

OL^Ot 
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with 



Ll,X 4^ ' 

X {S'M's, UM'j\3'M\)(j'M], l'm[\JMj) {SMs, LO\JMj) 

X {S'M's, I'MlIIMj) Q'M'j, Tm[\JMj) {SMs, LO\JMj) 

X {L'M'l, L' - M'L\LM){lm'i, T - m[\lm){lM, lm\\Q) 

X (L'0,Z'0|LO)(/'Oj'0|/~0) 

X {SMs\]^{Ms - m2), ^m2){SMs\^mi, ^{Ms - mi)) 

X (^(mi + qi), ki (-gi)|^mi)(^m2, k2 (-g2)|^(Ms - "^i)) 

1 11 1 - 

X (2K1 - ^3),A;3 q3\^m[){-m2,k4 g4|-(M^ -m'J) , (C.9) 

where the STim nms over {Ms- M'^. Af|. Af^. nii. v)'{\ and 



K — 



{2U + 1)(2P + 1)(2L + 1)(2L + 1){2L' + 1)(2L^ + l)(2A;i + 1)(2A;2 + 1) 
\ (2L + 1)(2[+1)(2J + 1)(2J+1) 



D On the non— factorization of a strong final state interaction 

In this appendix we will demonstrate the need for a consistent treatment 
of both the NN scattering and production amplitudes in order to obtain 
quantitative predictions of meson-production reactions. 

Let us assume a separable NN potential 

V{p',k) = ag{p')g{k) , (D.l) 

where a is a coupling constant and g{p) an arbitrary real function of p. With 
this potential the T-matrix scattering equation can be readily solved to yield 

T(v' k) = ^(P^'^) m 2) 

' l-R{p')+iK{p')V{p',p') ' ^ > 

with 

R{p') ^ mP dk' ^^p^ (D.3) 

and k{p) — pTT/i denotes the phase space density here expressed in terms of 
the reduced mass of the outgoing two nucleon pair = miv/2. 
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Note that for an arbitrary function g{k), such as g{k) = 1 as discussed be- 
low, R{p') may be divergent. In this case R is to be understood as prop- 
erly regularized. The principal value integral R(p') given above is therefore a 
model-dependent quantity, for it depends on the regularization scheme used. 
The condition that the on-shell NN scattering amphtude should satisfy Eq. 
6 relates this to the on-shell potential, V{p',p'): 

Rip') = 1 + Kip') cotiSip'))Vip',p') , (D.4) 

where it is assumed that r]ip') = 1. This shows that, for a given potential, the 
regularization should be such that Eq. (D.4) be satisfied in order to reproduce 
Eq. 6. Indeed, in conventional calculations based on meson exchange models, 
where one introduces form factors to regularize the principal value integral, 
the cutoff parameters in these form factors arc adjusted to reproduce the A^A^ 
scattering phase shifts through Eq. 6. Conversely, for a given regularization 
scheme, the NN potential should be adjusted such as to obey Eq. (D.4). 
This is the procedure used in effective field theories [311], where the coupling 
constants in the A^'A'" potential are dependent on the regularization. 

We also assume that the production amplitude M is given by a separable 
form, 

MiE,k)=Pgik)hip) , (D.5) 
where is a coupling constant and hip) an arbitrary function of the relative 
momentum p of the two nucleons in the initial state. With this we can express 
the total transition amplitude as 

Equation D.6 is the desired formula for our discussion. It allows us to study 
the relationship between the A'"A'" potential and the production amplitude 
M{E, p') exphcitly as different regularization schemes are used. For this pur- 
pose let us study the simplest case of a contact A'"A'" potential (setting the 
function g — 1) in the limit p' — > 0. If we regularize the integrals by means of 
the power divergent subtraction (PDS) scheme [311] we get 

1 — ap, 
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where denotes the regularization scale. Substituting this result into Eq. 
(D.4), we obtain 



for the NN coupling strength. Note that for = the PDS scheme reduces 
to that of minimal subtraction [311]. Since the total production amplitude A 
should not depend on the regularization scale we immediately read off Eq. D.6 
that 



Therefore the model clearly exhibits the point made in section 2: Namely, the 
necessity of calculating both the production amplitude and the FSI consis- 
tently in order to allow for quantitative predictions. 

E Chiral Counting for Pedestrians 

In this appendix we demonstrate how to estimate the size of a particular loop 
integral. This is a necessary step in identifying the chiral order of a diagram. It 
should be clear, however, that the same methods can be used to estimate the 
size of any integral. However, the importance of the chiral symmetry is that 
it ensures the existence of an ordering scheme that suppresses higher loops. 

The necessary input are the expressions for the vertices and propagators at 
any given order. For the chiral perturbation theory those can be found in Ref. 
[197]. In addition we need an estimate for the measure of the integral. 

Once each piece of a diagram is expressed in terms of the typical mo- 
menta/energies, one gets an estimate of the value of the particular diagram. 
The procedure works within both time-ordered perturbation theory and co- 
variant theory. Obviously, for each irreducible diagram both methods have to 

give the same answer. If a diagram has a pure two-nucleon intermediate state, 
as is the case for the direct production, the covariant counting can only give 
the leading order piece of the counting within TOPT. 

In this appendix we study only diagrams that are three-particle irreducible; 
i.e., the topology of the diagram does not allow an intermediate two-nucleon 




/? cx (1 — a/j,) 



-1 
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Fig. E.l. ^ typical loop that contributes to pion production in nucleon-nucleon col- 
lisions. Solid lines denote nucleons and dashed lines pions. The solid dots show 
the points of interactions. The horizontal dashed lines indicate equal time slices, as 
needed for the evaluation of the diagram in time-ordered perturbation theory. 

state to go on-shell. The reducible diagrams require a different treatment and 
are discussed in detail in the main text. There is another group of diagrams 
mentioned in the main text that is not covered by the counting rules presented, 
namely radiation pions. Those occur if a pion in an intermediate state goes 
on-shell. It is argued in the main text that these are suppressed, because the 
pion, in order to go on-shell, is only allowed to carry momenta of the order 
of the external momenta, and thus the momentum scale within the loop is of 
order of the pion mass and not of the order of the initial momentum. Therefore 
we do not consider radiation pions any further. 



E.l Counting within TOPT 



As mentioned above, if we want to assign a chiral order to a diagram, all we 
need to do is to replace each piece in the complete expression for the evaluation 
of the matrix element by its value when all momenta are of their typical 
size. In case of meson production in nucleon-nucleon collisions this typical 
momentum is given by the initial momentum pj. Time-ordered perturbation 
theory contains only three-dimensional integrals and thus we do not need to 
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fix the energy scale in the integral. 
The counting rules are 

• the energy of virtual pions is interpreted as 0{pi) and thus 

■ every time slice that contains a virtual pion is interpreted as 1/pi (see Fig. 
E.l) 

■ for each virtual pion line put an additional 1/pi (from the vertex factors), 

• interpret the momenta in the vertices as Pi 

• every time slice that contains no virtual pion is interpreted as l/mTr ; most 
of these diagrams, however, are reducible (c.f. main text) 

• the integral measure is taken as pf/(47r)^. 

Here we used that pf/Mjv — itLt^ ^ pi, in accordance with Eqs. (80), and thus 
nucleons can be treated as static in the propagators if there is an additional 
pion present. However if there is a time slice that contains two nucleons only, 
the corresponding propagator needs to be identified with the inverse of the 
typical nucleon energy 1/m^ and the static approximation is very bad [155]. 

In the diagram of Fig. E.l three ttA^A^ vertices, each ~ pj//^, appear as well 
as the 7171 NN Wcinbcrg-Tomozawa vertex pi/ f^. In addition the three time 
slices give a factor l/p^j and we also need to include a factor 1/p^, since there 
are two virtual pions. We therefore find 

~ [fj fl [Pij W ~ p. \mJ ■ 

Here we used ^Tif.,^ — Mj^ and Pi ~ ^/Wj^m^. 
E.2 Counting within the covariant scheme 

Naturally, as TOPT and the covariant scheme are equivalent, the chiral order 
that is to be assigned to some diagram needs to be the same in both schemes. 
The reason why we demonstrate both is that here we are faced with a problem 

There is one exception to this rule: if a time derivative acts on a pion on a vertex, 
where all other particles are on-shell, then energy conservation fixes the energy. 
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in which the typical energy scale 171.,^ and the typical momentum scale pi are 
different. In the covariant approach a four-dimensional integral measure en- 
ters and naturally the question arises whether or pj is appropriate for the 
zeroth component of this measure po- For example, in Ref. [163] it was argued, 
that one should choose m^r, although this choice is by no means obvious from 
the structure of the integrals. However, given the experience we now have in 
dealing with loops in time-ordered perturbation theory, where these ambigui- 
ties do not occur, the answer is simple: we just have to assign that scale to po 
that will reproduce the same order for any diagram as in the counting within 
TOPT [204]. Once the choice for po is fixed, the much easier to use covariant 
counting can be used to estimate the size of any loop integral. Thus have the 
following rules: 

• the energy of virtual pions is interpreted as 0{pi) 

• each pion propagator is taken as 0{l/pf) 

• each nucleon propagator that cannot occur in a two-nucleon cut is taken as 
O{l/po) (the leading contribution of a nucleon propagator that can occur 
in a two-nucleon cut is 0(l/m^); most of these diagrams, however, are 
reducible (c.f. main text)) 

• interpret the momenta in the vertices as pi 

• the integral measure is taken as Popf/(47r)^ (when the diagram allows for a 
two-nucleon cut the measure reads (m^p^)/(47r)^). 

Thus we have for the diagram of Fig. E.l 

" v/J Uv U; (4^)' n \mj [pJ ■ 

Thus, we need to assign po ~ Pi in order to get the same result in both schemes. 
As a side result we also showed, that the nucleons are indeed static in leading 
order inside loops that do not have a two-nucleon cut, as pointed out in Ref. 
[206]. 



^^With the same exception as in the time ordered situation (c.f. corresponding 
footnote) . 
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List of Figures 



1 Comarison of two different reactions that can be used to 

extract parameters of a strong final state interaction. Shown 
is a typical low momentum transfer reaction (diagram a), 
e.g., K~d '-/An) and a large momentum transfer reaction 
(diagram h), e.g., pp pK~^A). 

2 Some possible contributions to the production operator for 
neutral pion production in pp collisions. Solid lines are 
nucleons, dashed lines are pions. 

3 Sketch of the production reaction showing the initial state 
interaction amongst the two nucleons (IS I), the final state 
interaction of all outgoing particles (FSI) as well as the 
production operator A as defined in the text. The left diagram 
indicates the complexity possible, whereas the right diagram 
shows the first and potentially leading term for the FSI only. 

4 Illustration of the momentum mismatch. The picture shows 
why one should a priori expect meson exchange currents to be 
very important close to the threshold. 

5 Diagrammatic presentation of the distorted wave Born 
approximation. 

6 Energy dependence of two different spin cross sections for the 
reaction pp — > pp7r°. The data for ^ao (^(Ji), shown as solid 

( opaque ) dots, are compared to a fit of the form y — Arf with 
r = 6 (t — d)). The data are from table V of Ref. [47] with an 
additional error of 10 % on the total cross section. 
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7 Illustration of the sources of a discontinuity in the amplitude 
A from two baryon intermediate states. The dashed horizontal 
lines indicate the presence of a two baryon unitarity cut 
(Diagrams with initial state interaction are not shown explicitly 

(c.f Fig. 5)). 23 

8 The integration contours in the complex ml ^ plane to be used 
in the original treatment by Geshkenbein [55,56] (a) and in the 
derivation of Eq. (18) (h). The thick lines indicate the branch 

cut singularities. 26 

9 Illustration of the choice of variables in the over all 

center-of-mass system. 34 

10 Sketch of a Dalitz plot for a reaction with a three-body final 
state (particles are labeled as 1,2 and 3). Regions of possibly 
strong final state interactions in the (ij) system are labeled by 
F(ij-); a resonance in the {12) -system would show up like the 
band labeled as R. The region of a possible interference of the 

two is labeled as I. 35 

11 Angular distributions for the reaction pp dK^K^ measured 
at Q = A& MeV [66]. The solid line shows the result of the 
overall fit including both KK s-waves as well as p-waves. To 
obtain the dashed (dotted) line, the parameters for the p-wave 
(s-wave) were set to zero (see teoct). The small error shows 
the statistical uncertainty only, whereas the large ones contain 
both the systematic as well as the statistical uncertainty added 
linearly. 38 

12 Various mass distributions for the reaction pp — > dK^K'^ 
(Line code as in Fig. 11). The small error bars show the 
statistical uncertainty only, whereas the large ones contain both 
the systematic as well as the statistical uncertainty (c.f. Ref. 

[66]). 40 
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13 Some polarization observables reported in Ref. [47] for the 
reaction pp — > pp7r° at rj — 0.83 as a function of the pion angle 
compared to predictions of the model of Ref . [71]. The solid 
lines show the results for the full model whereas contributions 
from the Delta where omitted for the dashed lines. 48 

14 Some polarization observables reported in Ref. [47] for the 
reaction pp — ppn^ at rj = 0.83 as a function of cos{9p) 
compared to predictions of the model of Ref. [71]. The solid 
lines show the results for the full model whereas contributions 
from the Delta were omitted for the dashed lines. 49 



15 Sensitivity of the analyzing power as well as the differential 
cross section for the reaction pn — > ppn~ to the sign of the 
^Pq -^^ SqS amplitude ai. The lines corresponds to the model of 
Ref. [73]: the solid line is the model prediction whereas for the 

dashed line the sign of ai was reversed. The experimental data 

are from Ref [74] and [75] at TLab = 353 MeV (rj = 0.65^. 53 

16 Illustration of the coordinate system used in the analysis for 

the reaction pp — > dK^K^ . 59 

17 Demonstration of the selectivity of the spin cross sections. 
Shown are the spectra of the reaction pp prnx^ as a function 
of e for an excess energy of 25 MeV. In each panel the solid 
line shows the full result for the corresponding cross section, 
the dot-dashed, long-dashed, dashed, and dotted lines show the 
Bp, Ss, Sd, and Pp contribution respectively. The curves are 
from the model of Ref. [71]. 62 

18 The lowest meson production thresholds for single meson 
production in proton-proton collisions together with the 
corresponding energy ranges of the modern cooler synchrotrons. 63 



138 



19 Illustration of different sources of charge symmetry breaking: 
diagram a) shows CSB in the production operator through 

TT — T] mixing and diagram b ) shows CSB in the propagation of 
the scalars. Thin solid lines denote nucleons, thick solid lines 
scalar mesons and dashed ones pseudoscalar mesons. The X 
indicates the occurrence of a CSB matrix element. 68 

20 Comparison of the model predictions of Ref. [73] to the data. 
The references for the experimental data can be found in table 
8. The solid lines show the results of the full model; the dashed 

line shows the results without the A contributions. 75 

21 The NN phase shifts for the model of Ref . [175] in the energy 
range relevant for pion production. The pion production 
threshold is at E^ab — 286 Me V. The experimental data are 

from Refs. [177] (triangles) and [61] (circles). 76 

22 Contributions to the potential of the model of Ref. [179]. 78 

23 The ttN phase shifts for the model of Ref . [179] in the energy 
range relevant for pion production. The experimental data are 
from Refs. [183,184]. Note the different scales of the various 
panels. 79 

24 Comparison of the model predictions to the data taken from 



26 Illustration of the role of the ANn contact term in NN — > NNn 
and three nucleon scattering. Solid lines denote nucleons, 



Ref [47] (pp ppTi"), Ref [108] (pp d7r+) and Ref [94] 
(pp priTi^ ). 
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25 Illustration of resonance saturation. 
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dashed lines denote pions. 
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27 Comparison of the predictions from effective field theory at LO 
(dashed line) and NLO (with the ttN parameters from a NLO 
(dotted line) and an NNLO (solid line) analysis) with the ^ai 
cross section for pp ppir^ [4 7]. 92 

28 of pp npiT^ in chiral perturbation theory. The different 
lines correspond to values of the parameter related to the 
three-nucleon force: 6=1 (long dashed line). 6 = (dot-dashed 
line), 6 = —0.2 (solid line), and 6 = —1 (short dashed line). 

Data are from Ref [219]. 93 

29 Visualisation of Eq. (90). 100 

30 The one sigma exchange as it is perturhatively built up in the 
effective field theory, starting from the left with the lowest 
order diagram (NLO). The chiral order increases by one power 

in X between each diagram from the left to the right. 102 

31 Leading cut structure of the one pion exchange potential in 
nucleon-nucleon scattering. 103 

32 Results for the total cross sections of the reactions pp pprj, 
pn — * pnrj, and pn drj from Ref. [34] employing different 
NN models for the final state interaction. The solid lines 
represent the results with the CCF NN model [175] whereas 
the dashed-dotted lines were obtained for the Bonn B model 

[25]. Data are from Refs. [33,109-112,120]. 108 

33 Invariant mass spectrum for the two-proton system for the 
reaction pp pprj at Q=15 MeV. The dashed line shows the 
Ss contribution, the dot-dashed line the Ps distribution and 
the solid line the incoherent sum of both. The data is taken 

from Ref [113]. 110 
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34 Comparison of the quality of available data for the reactions pA 
elastic scattering (data are from [264-266]) and pp — > K^Kp 
at TLab^2.3 GeV [267] . In both panels the red curve 
corresponds to a best fit to the data. In the left panel the dashed 
lines represent the spread in the energy behavior allowed by the 
data, according to the analysis of Ref. [266]; analogous curves 
in the right panel would lie almost on top of the solid line and 

are thus not shown explicitly. 113 

35 Values allowed for the spin singlet and spin triplet scattering 
length by the AN elastic scattering data according to Ref. 
[266]. The dark shaded area denotes the la range for the 
parameters and the light shaded area the 2a range. The cross 
shows the best fit value (as — 1.8 fm and at — 1.6 fm). 114 

36 Graphical illustration of the leading contribution to the /o — Qq 
mixing matrix element A defined in Eq. (91). 116 

37 Modulus of the leading piece of the mixing amplitude A defined 
in Eq. (91). The two kinks occur at the K^K'^ (at 987.35 

MeV) and the K°K° (995.34 MeV) threshold respectively. 117 

38 Possible diagrams to contribute to the reaction Nd 3Nx. 
In diagram a) and b) the three nucleons in the final state 
from a bound state, whereas diagram c) shows the dominant 
diagram in quasi-free kinematics. Here it is assumed that only 
a particular nucleon pair ( to be identified through proper choice 

of kinematics) interacts in the final state. 119 



^ The data was taken inclusively (only the kaon was measured in the final state), 
however, for the small invariant masses shown only the AK channel is open. 
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39 Momentum distribution of proton- spectators in the reaction 
pd— >37r"'"27r~p. The data are taken from Ref. [305] while 
the solid line is is based on a calculation that includes the 
quasi-free production as well as meson rescattering. The figure 
is taken from [306], where also details about the corresponding 
calculation can be found. 121 

E.l A typical loop that contributes to pion production in nucleon- 
nucleon collisions. Solid lines denote nucleons and dashed 
lines pions. The solid dots show the points of interactions. The 
horizontal dashed lines indicate equal time slices, as needed 
for the evaluation of the diagram in time-ordered perturbation 
theory. 133 
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